Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	42.0%	140.1 / 14.9	82	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (37/10);  new base C (0/34);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 3.82e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGC  >  NC_000913/907456‑907750                                                                                                                                                   |                                                                                                                                                     aGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGc                                                                                                                                                  >  2:1234861/1‑151 (MQ=255)  gcgcGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCt                                                                                                                                                 >  3:1586516/1‑151 (MQ=255)   cgcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTc                                                                                                                                                <  2:1474797/151‑1 (MQ=255)      ggTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCAcc                                                                                                                                             <  2:400175/151‑1 (MQ=255)       gTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGGACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCAc                                                                                                                                              >  2:1010950/1‑149 (MQ=255)        tATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGc                                                                                                                                           >  3:1288613/1‑151 (MQ=255)         atcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGcc                                                                                                                                          <  3:1433490/151‑1 (MQ=255)         atcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGc                                                                                                                                           <  2:730830/150‑1 (MQ=255)           catcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCa                                                                                                                                         <  2:1230487/150‑1 (MQ=255)              ccGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATaa                                                                                                                                      >  3:1093028/1‑150 (MQ=255)               cGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGGGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATaaa                                                                                                                                     >  3:848379/1‑150 (MQ=255)                 ccAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATg                                                                                                                                  >  2:1457564/1‑151 (MQ=255)                      agATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCa                                                                                                                              >  3:126108/1‑150 (MQ=255)                      agATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCa                                                                                                                              >  3:417198/1‑150 (MQ=255)                        aTCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAAc                                                                                                                           >  2:529330/1‑151 (MQ=255)                              tttCACCCCCTAAACCACCGCCCGGGGGGGGGTGGGCGGTACGAAAGTAACCACATTGAAAACCCACCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTg                                                                                                                         <  2:1635960/147‑1 (MQ=255)                              ttccacccccTAAATCACCCCCCGGGGGGGGGGGGGGGGTACGAAAGAAATCAATTTGAGAACCCACACAAAAACCCGGGGGGTGGGGGGGGGGTGGGGAATAACCCGGCGGCTCTTCCCCCTCACCGCCATAAAATGTTCAAACTGtt                                                                                                                       <  2:1388338/143‑1 (MQ=255)                                 cacccccTAAACCACCGCCCGGGGGGGGGGGGGGGGTAGGGAAGAAATCAAATGGAAAACCCCCCCCAAAACCCAGGGGGGGGGGGGGGGGAGGGGAACCAACCAGGGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGa                                                                                                                  <  2:1321053/148‑1 (MQ=255)                                 cAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                                                   >  2:817406/1‑150 (MQ=255)                                    ccccaaaaTCACCGCCCGGGGGGGGGGGGGCGGTAGGAAAGTAATCACATGGAAAACCCCCCCCAAAACCCAGGGGGGGGGGGGGGGGATGGGGACCAACCAGGGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATaa                                                                                                               <  3:1391789/145‑1 (MQ=255)                                    ccccTGAATCAACGCCCGGGGGGTGGTGGGCGGTAGGAAAGTAATCACATGGAAAACCCACCCCAAAACCCGGGGGGTGGGGAGGGGGAGGGGGATCAACCAGGGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                                                   <  2:970160/147‑1 (MQ=255)                                     cccTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATaa                                                                                                               >  3:823750/1‑150 (MQ=255)                                     cccTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAAc                                                                                                              >  3:362150/1‑151 (MQ=255)                                       cTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                                            >  2:255742/1‑151 (MQ=255)                                        tGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                                            >  2:971361/1‑150 (MQ=255)                                          aaTCACCGCCCTGGCGGGGGGGGGGGGTACGGAAGTAATCACATGGAGAACCCACCCCAAAACCCGGCGGGGGGGGGGGGGGAGGGGGATCAACCGGCGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGa                                                                                                           <  2:1624158/149‑1 (MQ=255)                                          aaTCACCGCCCGGGGGGTGGGGGGGGGTAGGGAAGAAATCACATGGAAAACCCCCCCCAAAACCCGGGGGGTGGGGGGGGGGAGGGGGATAAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGa                                                                                                           <  2:1004154/149‑1 (MQ=255)                                          aaTCACCGCCCGGGGGGGGGGGGGCGGTACGAAAGAAATCACATTGAAAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGGATGGGGATCAACCGGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGc                                                                                                         <  2:1004617/151‑1 (MQ=255)                                          aaTCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGc                                                                                                         >  3:216404/1‑151 (MQ=255)                                              accccccGGGCGGGGGTGGGCGGAACGAAAGTAACCCTATTGAAAAACACCCCAAAAACCCGGGGGGGGGGGGGGGGGAGGGGGACAAACCAGCGCCGTTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGa                                                                                                           <  2:742061/141‑1 (MQ=255)                                                cGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                                   >  2:1293763/1‑151 (MQ=255)                                                cGCCCTGGCGGGGGGGGGCGGTACGGAAGTAATCACATTGAAAACCCACACCGAAACCCGGGGGGGGGGGGGGGGGATGGGGATCAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                                   <  2:758001/151‑1 (MQ=255)                                                 gCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTg                                                                                                  >  2:921487/1‑151 (MQ=255)                                                  cccTGGGGGTGGTGGGGGGTACGGAAGTAATAACATGGAGAACCCCCCCCGAAACCCGGGGGGGGGGGGGGGGGAGGGGAACAAACCGGGGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTg                                                                                                  <  3:1072287/150‑1 (MQ=255)                                                      gggggtgggGGGCGGTACGAAAGAAATCACATTGAAAACCCCCACCAAAACCCAGGGGGGGGGGAGGGGGAGGGGGATCACCCGGCGGCGCTTCCCGCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTc                                                                                               <  3:895527/146‑1 (MQ=255)                                                      gggggggggggggggTAGGGAAGTAATCACATAAAAAAACCCCCCCAAACCCCGGGGGGGGGGGGGGGGGATGGGAATAACCCGGCGGCGCTCCCCGCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                                             <  3:808777/142‑1 (MQ=255)                                                      ggggggggggggCGGTGGGGAAGTAATCAGAGGGAGAAACCACGCAGAAATCCGGGGGGGGGGGAGGGGGAGGGGAAAAAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                                              <  2:328309/141‑1 (MQ=255)                                                                ggCGGTAGGGAAGTAATCCCATTGAGAACCCCCCCCAAAACCCAGGGGGGGGGGGGGGGGATGGGGACCAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGa                                                                                   <  2:1234433/151‑1 (MQ=255)                                                                     tACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATg                                                                              >  3:1467077/1‑151 (MQ=255)                                                                      aCGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGc                                                                             >  3:493310/1‑151 (MQ=255)                                                                             tAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                      >  3:420110/1‑151 (MQ=255)                                                                             tAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                      >  2:1138883/1‑151 (MQ=255)                                                                              aaTCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCg                                                                     >  3:449871/1‑151 (MQ=255)                                                                              aaTCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCg                                                                     >  3:871086/1‑151 (MQ=255)                                                                              aaTCACATTGAAAACCCCCCCCGAAACCCAGGGGGGGGGGGGGGGGAGGGGGATCAACCGGGGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtgt                                                                        <  2:1325940/148‑1 (MQ=255)                                                                               accacaTGGAAAACCCACCCCGAAACCCAGGGGGGGGGGGGGGGGGTGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGc                                                                    <  2:367309/149‑1 (MQ=255)                                                                                tCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCt                                                                   >  3:169205/1‑151 (MQ=255)                                                                                  acaTTGAAAACCCACACCGAAACCCAGGGGGGGGGGGGGGGGAGGGGAATCAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTg                                                                  <  3:941103/150‑1 (MQ=255)                                                                                      tGAAAACCCACCCCAAAACCCGGGGGGTGGGGAGGGGGAGGGGAATCAACCAGCGGCGTTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGc                                                                <  2:1451658/148‑1 (MQ=255)                                                                                       gagaAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                             >  2:1088459/1‑150 (MQ=255)                                                                                          aaaCCCCCCCGAAACCGGGGGGGGGGGGGGGGGGAGGGGGACCACCCGGCGGCGCTTACCGCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatg                                                          <  3:1488446/150‑1 (MQ=255)                                                                                           acccACACCAAACCCCGGGGGGGGGGGGGGGGGATGGGAATAACCCGGCGCCGTTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatgat                                                        <  3:998074/149‑1 (MQ=255)                                                                                                 acCGAAATCCAGCGGGTGGAGCGGGGGTTGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGACCAACTCTGGTCCGGCGGCTTCGAAGATGCGGCGTTCGCTGGCCGGATGATCTaac                                                   >  3:633291/1‑150 (MQ=255)                                                                                                    gAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                                               >  3:1483746/1‑151 (MQ=255)                                                                                                     aaagccgggggggggggggggggggggggAACAACCGGCGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaaca                                                  <  3:848428/138‑1 (MQ=255)                                                                                                     aaacccgggggggggggggggggAGGGGGATAAACCGGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaac                                                   <  2:612131/137‑1 (MQ=255)                                                                                                     aaacccggggggTGGGGGGGGGGAGGGGGATCACCCAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                                               <  2:1371890/141‑1 (MQ=255)                                                                                                       aTCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                             >  2:1150395/1‑150 (MQ=255)                                                                                                       aTCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                             >  3:122048/1‑150 (MQ=255)                                                                                                         ccgggggggggggggggggAGGGGGATCACCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCaa                                            <  2:1255698/144‑1 (MQ=255)                                                                                                            gggggTGGGGAGGGGGAGGGGGATCACCCGGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCa                                         <  2:834457/147‑1 (MQ=255)                                                                                                             gggggggggggggggAGGGGGATCAACCGGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGAt                                      <  2:558085/150‑1 (MQ=255)                                                                                                              ggggggggggggggAGGGGGATCAACCGGGGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGa                                       <  2:1586677/149‑1 (MQ=255)                                                                                                              ggggggggAGGGGGATGGGGATCACCCGGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATa                                     <  2:1288733/151‑1 (MQ=255)                                                                                                              gggTGGGGGGGGGGATGGGGATAACCCGGCGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGa                                       <  2:46931/149‑1 (MQ=255)                                                                                                              gggTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGAt                                      >  3:1453433/1‑150 (MQ=255)                                                                                                               gggggggggggggggggGAATAACCCGGGGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAg                                    <  2:347607/148‑1 (MQ=255)                                                                                                                ggggggggggggAGGGGGATCACCCAGCGGCGTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGt                                   <  1:4370/149‑1 (MQ=255)                                                                                                                  ggggggggggAGGGGGATCACCCAGGGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgc                                 <  2:1622808/151‑1 (MQ=255)                                                                                                                  ggggggggggAGGGGGACCACCCAGGGGCTCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTg                                  <  2:772520/150‑1 (MQ=255)                                                                                                                    ggggggggATGGGGACCACCCGGCGGCCCTTCCCGCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcag                               <  3:188013/148‑1 (MQ=255)                                                                                                                    ggAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgc                                 >  3:1538886/1‑149 (MQ=255)                                                                                                                     gggggggagggggaTCACCCAGGGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagc                              <  2:361071/149‑1 (MQ=255)                                                                                                                       gggggAGGGGAATCAACCGGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCaa                            <  2:437189/151‑1 (MQ=255)                                                                                                                                ggatCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                     >  3:881929/1‑149 (MQ=255)                                                                                                                                       cggggggctCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAAt             <  3:72386/145‑1 (MQ=255)                                                                                                                                       ccAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATg            >  3:495187/1‑151 (MQ=255)                                                                                                                                       ccAGCGGCCTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATg            <  2:767436/151‑1 (MQ=255)                                                                                                                                        cAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATgg           >  2:1040734/1‑151 (MQ=255)                                                                                                                                        cAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATgg           >  3:1427509/1‑151 (MQ=255)                                                                                                                                          gcggccCTCGCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCTAGCGCAATgg           <  3:1452675/149‑1 (MQ=255)                                                                                                                                                 ttACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGc  >  3:513055/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                     AGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGC  >  NC_000913/907456‑907750

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