Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,122,676	T→A	8.8%	M102L (ATG→TTG)	cpsG ←	phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,122,676	0	T	A	8.8%	153.5 / 6.0	64	M102L (ATG→TTG)	cpsG	phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base T (31/27);  new base A (2/4);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

GTTGATTTGCTGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTT  >  NC_000913/2122530‑2122824                                                                                                                                                   |                                                                                                                                                     gTTGATTTGCTGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCg                                                                                                                                                    >  2:470153/1‑149 (MQ=255) gTTGATTTGCTGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCg                                                                                                                                                    <  1:402272/149‑1 (MQ=255) gTTGATTTGCTGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCg                                                                                                                                                    <  2:198666/149‑1 (MQ=255)    gATTTGCTGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGAt                                                                                                                                                 <  2:722155/149‑1 (MQ=255)      tttGCTGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTa                                                                                                                                               <  1:650823/149‑1 (MQ=255)        tGCTGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGCTCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTAtg                                                                                                                                             >  1:1146596/1‑149 (MQ=255)        tGCTGATAGCGACCGCGTTTGGGTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCAGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTAtg                                                                                                                                             <  1:286631/149‑1 (MQ=255)         gCTGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGAGCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTAtgg                                                                                                                                            <  1:71811/149‑1 (MQ=255)          cTGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCGATGGATTAtaac                                                                                                                                           <  1:401481/149‑4 (MQ=255)           tGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCGATGGATTAtaacg                                                                                                                                          <  1:170199/149‑5 (MQ=255)           ggATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCGATGGATTAtaacg                                                                                                                                          <  2:441366/148‑5 (MQ=255)                   aCCGCGTTTGGTTTCTTCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGTAAGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTGATGCCGTTCTAATCCATCGGATTATGGCTGGCGGTaa                                                                                                                                  <  1:857122/149‑1 (MQ=255)                   aCCGCGTTTGGTTTCATCGCCGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTaa                                                                                                                                  >  1:780340/1‑149 (MQ=255)                    ccGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGACCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAAc                                                                                                                                 >  1:117437/1‑149 (MQ=255)                        gTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGGAACTTCa                                                                                                                             >  1:260898/1‑149 (MQ=255)                              tttCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGCGCCCCCTCGCGAACCAGCAACATGACGTTATAATCCATCGGATTATGTCTGGCGGTAAATTCAAttcac                                                                                                                       >  2:1078961/1‑145 (MQ=255)                              tttCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCGATGGATTAtaa                                                                                                                                            >  2:1124383/1‑126 (MQ=255)                              tttCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCGATGGATTAtaa                                                                                                                                            <  1:1124383/128‑3 (MQ=255)                              tttCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATgccg                                                                                                                       >  2:4679/1‑149 (MQ=255)                                 cATCGACGGGAGGAAAGTCTTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCGGATCGGGCGAGCCCCCTCGCGAACCAGCTGCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCa                                                                                                                    <  1:143459/149‑1 (MQ=255)                                 cATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCa                                                                                                                    <  1:512624/149‑1 (MQ=255)                                   tCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATc                                                                                                                  >  1:1147498/1‑149 (MQ=255)                                   tCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATc                                                                                                                  >  1:823100/1‑149 (MQ=255)                                           aggaAAGTCGTTGGCTTCAGCCAGACGCTGGCCGACGCGCCTTCCGGTTTAACCGAGGAGCGGGCGAGCCCCGGCGCGACCCGCCTAAATGCCTTTATAAGCCAACGGATTAGGGGTGGCGGTACCTGAATTGCGGCCTGCCACGcggc                                                                                                          >  2:159572/1‑146 (MQ=255)                                           aggaAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCgc                                                                                                          >  1:643205/1‑148 (MQ=255)                                           aggaAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCga                                                                                                          <  2:279213/149‑1 (MQ=255)                                              aaaGTCGTTGGCTTCAGCCAGACGCTGGACGTAGCGCAGTCCGGTATCTCCGCTGATCGGTCGAGCCACCTCGCGAACCAGCTTCATGCCCTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCAACCACGCCGAGAt                                                                                                       >  2:1142415/1‑149 (MQ=255)                                                  tCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGaa                                                                                                   <  2:605006/149‑1 (MQ=255)                                                    gTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGCTCGGGCGAGCCCCCTCGCGAGCCAGCGTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAAAGCCGCCATCCACGCCGAGATGGAACg                                                                                                 >  1:666559/1‑149 (MQ=255)                                                      tGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTg                                                                                               <  2:441282/149‑1 (MQ=255)                                                           tCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGaa                                                                                          <  2:1147498/149‑1 (MQ=255)                                                              gCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATa                                                                                       <  1:604529/149‑1 (MQ=255)                                                               ccAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAg                                                                                      >  1:542118/1‑149 (MQ=255)                                                                        ggACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCg                                                                             >  2:718063/1‑149 (MQ=255)                                                                        ggACGTCGCGCAGTCCGGTATCTCCACTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTCGCGAAATAGATCTCTTCg                                                                             >  1:876759/1‑149 (MQ=255)                                                                          aCGTCGCGAAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGt                                                                           <  2:687301/149‑1 (MQ=255)                                                                           cGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCTTCGGATTATGGCTGGCGGTTCCTGCAATGCCGGCATCCACGCGGAGATGGCGCGTGGCGCAATTGATCTCTTAggag                                                                          >  2:500350/1‑147 (MQ=255)                                                                              cgcgCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCg                                                                       <  1:470153/149‑1 (MQ=255)                                                                               gcgcAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCgg                                                                      >  2:202409/1‑149 (MQ=255)                                                                                 tcAGTCCGGTATCTCCTCTGATCGGGCGAGCCCCCTCGCGAACCATCTTCATGCCGTTTTAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGAc                                                                    <  2:831879/148‑1 (MQ=255)                                                                                 gcAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGAc                                                                    <  2:542118/149‑1 (MQ=255)                                                                                     tCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATAc                                                                <  2:780340/149‑1 (MQ=255)                                                                                       cGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCa                                                              <  2:341086/149‑1 (MQ=255)                                                                                         gTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAat                                                            >  2:817386/1‑149 (MQ=255)                                                                                             ctcCGCTGATCGGGCGAGCCCCCTCGCGCACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCCTCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGACGGACATACCAATATcc                                                        >  1:794194/1‑149 (MQ=255)                                                                                               ccGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAg                                                      <  1:302023/149‑1 (MQ=255)                                                                                                  cTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCAc                                                   <  1:4679/149‑1 (MQ=255)                                                                                                      tCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCCTCCACGCCGAGATGGAACGGGGCGAAATAGATCTCTTCGGTGCCGGACATACCAAGGTGCAGCACGTcc                                               >  1:96591/1‑148 (MQ=255)                                                                                                        gggCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAAc                                             >  1:614523/1‑149 (MQ=255)                                                                                                         ggCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACg                                            >  1:565090/1‑149 (MQ=255)                                                                                                                 cccTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATc                                    <  1:470063/149‑1 (MQ=255)                                                                                                                 cacTCGCGAACGAGCTTCATGCCGTTATAATACATCTGATTATGGCTGGCGGTAACTTCAATGCCGCAATCAACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGTCATACCAATATCCAGCACGTCAACGCCCGCATc                                    <  1:1078961/147‑1 (MQ=255)                                                                                                                  ccTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGTTGCCGGACATACCAATATCCAGCACTTCAACGCCCGCATcc                                   >  2:1093155/1‑149 (MQ=255)                                                                                                                  ccTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATcc                                   >  2:548026/1‑149 (MQ=255)                                                                                                                    tCGCGAACCAGCCTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGAAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGATATAGATCTCTACGGTTCCGGACATCCCAATATCCAGCACGTCAACGCCCGCATccag                                 >  1:1149957/1‑147 (MQ=255)                                                                                                                        taaCCTGCTTAATGCCGTATTAATCCATCGTATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGTAACGTGGCGAAATATTTCTCTTCGGTGCCGGACATACCAATTTCCCGCACGTCAACGCCCGCATCCTGTaaa                             <  2:876759/148‑1 (MQ=255)                                                                                                                           ccAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCt                          >  1:628206/1‑149 (MQ=255)                                                                                                                                       cGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCATTGCCGCCATCCACGCCGAGATGGAACGTGCCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCACCGCCCGCATCCTGTAAACATTAcgccagcgcc              >  1:1180559/1‑149 (MQ=255)                                                                                                                                             aaTCCATCGGATTATGGCTGGCGGT‑ACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTa       <  1:832519/149‑1 (MQ=255)                                                                                                                                              aTCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTa       <  2:1146596/149‑1 (MQ=255)                                                                                                                                                 cATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAgt    >  2:1136118/1‑148 (MQ=255)                                                                                                                                                 cATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAgg    >  2:448275/1‑149 (MQ=255)                                                                                                                                                  aTCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAggg   >  1:374508/1‑148 (MQ=255)                                                                                                                                                  aTCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGt   >  2:622125/1‑149 (MQ=255)                                                                                                                                                   tCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGtt  >  1:171148/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GTTGATTTGCTGATAGCGACCGCGTTTGGTTTCATCGACGGGAGGAAAGTCGTTGGCTTCAGCCAGACGCTGGACGTCGCGCAGTCCGGTATCTCCGCTGATCGGGCGAGCCCCCTCGCGAACCAGCTTCATGCCGTTATAATCCATCGGATTATGGCTGGCGGTAACTTCAATGCCGCCATCCACGCCGAGATGGAACGTGGCGAAATAGATCTCTTCGGTGCCGGACATACCAATATCCAGCACGTCAACGCCCGCATCCTGTAAACCTTTCGCCAGCGCCAGTTTTAAGGTT  >  NC_000913/2122530‑2122824

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑