Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,826,462	0	G	A	45.1%	60.6 / 38.3	71	G24E (GGA→GAA)	srlE	glucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base G (27/12);  new base A (0/32);  total (27/44)
Fisher's exact test for biased strand distribution p-value = 2.19e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGG  >  NC_000913/2826315‑2826609                                                                                                                                                      |                                                                                                                                                     gATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCTTTCAACTTGAATAAAAAGTTCACCTGGCTGGAGAAACATGATGACGCATATTCGGATCGTAATAGGAAT‑GGGCGGCTGGGGCGGCTGGTTTGCTCTG‑TAAGCCACGACGGGaaa                                                                                                                                                   >  1:160442/1‑151 (MQ=255) gATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGaaa                                                                                                                                                   >  2:673975/1‑151 (MQ=255) gATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGaaa                                                                                                                                                   >  2:191271/1‑151 (MQ=255) gATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGaaa                                                                                                                                                   >  1:318523/1‑151 (MQ=255) gATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGaaa                                                                                                                                                   >  1:15177/1‑151 (MQ=255)   tCTGACAACCGCCATGTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGaaaaa                                                                                                                                                 <  1:85825/151‑1 (MQ=255)           ccGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGt                                                                                                                                         >  2:672586/1‑150 (MQ=255)                aTTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTatat                                                                                                                                   >  2:1270110/1‑151 (MQ=255)                aTTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTatat                                                                                                                                   <  1:402330/151‑1 (MQ=255)                aTTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCGG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTatat                                                                                                                                   <  1:1122726/151‑1 (MQ=255)                 tttttGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATc                                                                                                                                  <  2:131095/151‑1 (MQ=255)                   tttGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCAc                                                                                                                                <  1:814070/151‑1 (MQ=255)                     tGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCAccg                                                                                                                              <  2:1118402/151‑1 (MQ=255)                     tGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCAccg                                                                                                                              <  1:1185594/151‑1 (MQ=255)                           aaaGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTa                                                                                                                        >  2:54349/1‑151 (MQ=255)                              gATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTAcc                                                                                                                      >  1:632862/1‑150 (MQ=255)                                ggcaatgccaagataaacaaacaGTTCGCCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCgg                                                                                                                   >  1:1366375/16‑151 (MQ=255)                                    cATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGc                                                                                                                  >  1:190689/1‑148 (MQ=255)                                        cAACTTGAAAAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCGGGGGCGGCCCGCTGGGGCTA‑AAAGCCCCGCGGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGAt                                                                                                           <  1:1283845/151‑1 (MQ=255)                                        aaaCTTGAAAAAAAGGTTCACGTGGGGGGAGAAACATCAGGACGCATATTCGGATCAAAAAGGGAGG‑GGGGGGCGGGGGGGGCCCGCTTGAGCTG‑AAAGCCACGCCGGAAAAAAAAA‑TCGTCTAAATCACCGCCGGTACCCGGCCTGCGAt                                                                                                           <  1:1405634/150‑1 (MQ=255)                                         aaCTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGAtt                                                                                                          >  2:1168825/1‑151 (MQ=255)                                                aaaaaaaGTTCTCCGGGAGGGAGCAACATCATGCCG‑GTTTTAGGGGGAAAAAAGGAAA‑GGGGGGCGGGGGGCGCCCCGTTGGGTTGAAAAGCCCCGCGGGAAAAAAAAATTCGTCTATACCACCGCCGGTACCCGGCCTGCGATTGTTGAc                                                                                                    <  1:904037/149‑1 (MQ=255)                                                   aaaaGTCCCCCTGGCGGGAGCAACACCATGACCCATATTGGAATGGAAAAAGGACC‑GGGGGGCTGGGGCGGCCCGCTTGAGCTA‑AAAGCCCCGCCGGAAAAAAAAA‑TCGTTTATACCACCGCCGGTACCCGGCCTGCGATTGTTGACAAAc                                                                                                <  2:338767/151‑1 (MQ=255)                                                         tccccgggCGGGAGAAACATCAGGACCCATTTCCGGACGAAAAAAGGAAG‑GGGTGGCGGGGGCGCCCCGCTGGAGCGA‑AAGGCCCCGCGGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGcac                                                                                          <  1:456874/145‑1 (MQ=255)                                                         tccccgggCGGGAACAACATCAGGCCGCATTTTCGAACCAAAAAAGAAGC‑GGGTGGCGGGGGGGGCCCGCTTGAGCTG‑AAACCCACGCCGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGcac                                                                                          <  2:560679/145‑1 (MQ=255)                                                         tccccgggCAGGACCAACATCAGGACCCATTTTCGAATCAAAAAAGGAAC‑GGGTGGCTGGGGCGCCCCCCTTGAGCTA‑AAACCCCCCCCGGAAAAAAAAA‑CCGTCTATACCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGcac                                                                                          <  2:593709/145‑1 (MQ=255)                                                          cACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGcaca                                                                                         >  1:459039/1‑151 (MQ=255)                                                              tGGCGGGAGCAACACCATGCCGCATATTCGGATCGAAAAAGGACG‑GGGGGGCGGGGGCGGCCCGCTTGGGCTA‑AAAGCCACGCGGGAAAAAAAAA‑TCGTCTATATCCCCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCtt                                                                                     <  2:630003/151‑1 (MQ=255)                                                               ggCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTa                                                                                    >  2:562902/1‑151 (MQ=255)                                                                    gtgCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCAACGCCGGGAAAAAAAA‑TCGCCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGAAAAACTGGCACAGCTTActgg                                                                                >  1:1155135/3‑150 (MQ=255)                                                                    gagaacaatcatGACGAATATTGGGAACAAAAAAGAAAC‑GGGTGGCGGGGGCGGCCCGCTTGAGCTA‑AAACCCCCGCGGGAAAAAAAAA‑CCGTCTATACCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTActgg                                                                                <  1:644867/143‑1 (MQ=255)                                                                    gAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAATCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTActgg                                                                                >  1:1015569/1‑150 (MQ=255)                                                                    gAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTActgg                                                                                >  1:281899/1‑150 (MQ=255)                                                                     ggCAACATCAGGACGAATATGGGGATCAAAAAAGGAAG‑GGGTGGCGGGGGCGGCCCGCTTGAGCTA‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTActggc                                                                               <  2:990247/149‑1 (MQ=255)                                                                     aGCAACATCATGACGCATATTCGGATCGAAAAAGGAGC‑GGTGGGCTGGGGGGCCCCGCTGGAGCTA‑AAAGCCGCGCCGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTActggc                                                                               <  1:1059468/150‑1 (MQ=255)                                                                     aGCAACATCAGGACGCATATTCGGACCAAAAAGGGAAG‑GGGTGGCTGGGGCGGCCCGTTTGAGCTA‑AAAGCCCCCCCGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTActggc                                                                               <  1:912156/150‑1 (MQ=255)                                                                      gaaaCATCAGGACGAATATTCGGATCAAAAAAGGAAC‑GGGGGGCGGGGGCGCCCCGCTTGAGCGA‑AAAGCCACGCCGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTActggct                                                                              <  2:15191/148‑1 (MQ=255)                                                                        aaCATCAGGCCGAATATTCGGAGCAAAAAAGGAAG‑GGGGGGCGGGGGCGGCCCGCTTGGGCTG‑AAACCCACGCGGGAAAAAAAAA‑TCGTTTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACtggctggc                                                                           <  2:1020298/151‑1 (MQ=255)                                                                         aCATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCa                                                                          >  1:734970/1‑151 (MQ=255)                                                                         aCATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCCATATCAACTCCGGTACCCGGCCTGCGCTGGTTGACAACCTGGCACAGCTTACTGGCTGTCa                                                                          >  1:681137/1‑151 (MQ=255)                                                                          attcatGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCGG‑AAAGCCACGCGGGAAAAAAAAA‑TCGTCTATATCCCCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAg                                                                         <  1:94073/149‑1 (MQ=255)                                                                           atcatGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAgg                                                                        >  1:659117/1‑151 (MQ=255)                                                                                  cgaattttCGGATCAAAAAAGAACC‑GGGTGGCTGGGGCGCCCCGCTTGAGCTA‑AAACCCCCGCCGGAAAAAAAAA‑TCGTCTATACCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGa                                                                 <  2:878407/145‑1 (MQ=255)                                                                                    cATATTGGGATCAAAAAAGAAAG‑GGGTGGCGGGGGCGGCCCGTTTGAGCTG‑AAAGCCCCGCCGGAAAAAAAAA‑ACGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGAc                                                                <  2:47559/150‑1 (MQ=255)                                                                                    aatatTCGGACCAAAAAAGGAGG‑GGGGGGCGGGGGCGCCCCGCTGGAGCTA‑AACGCCCCGCCGGGAAAAAAAA‑CCGTCTATAGCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGAc                                                                <  2:817166/149‑1 (MQ=255)                                                                                      tatcgggacagAAAAAGGAGG‑GGGTGGTGGGGGCGGCCCGTTTGGGTAA‑AAAGCCACGCGGGAAAAAAAAA‑TCGTCTAAATCTCCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGa                                                             <  2:1223769/141‑1 (MQ=255)                                                                                       ttttgggatgAAAAAGGGACGGGGGGGGGGGGGCGGCCCGTTGGAGCGA‑AAAGCCACGCGGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGa                                                             <  2:1136805/142‑1 (MQ=255)                                                                                             aTCAAAAAAGGAAG‑GGGGGGGGGGGGCGGCCCGCTTGGGCGA‑AAACCCACGCCGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAaag                                                      <  2:281026/151‑1 (MQ=255)                                                                                             aTCAAAAAAGAAGC‑GGGTGGCTGGGGGGCCCCGTTGGGGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCCCCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAaag                                                      <  1:1296185/151‑1 (MQ=255)                                                                                              tCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAaaga                                                     >  1:1112539/1‑151 (MQ=255)                                                                                              ccaaaaaaGGAGG‑GGGGGGCGGGGGCGCCCCGCTGGGGCTA‑AAGCCCCCGCGGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAaaga                                                     <  2:31710/148‑1 (MQ=255)                                                                                               caaaaaaGGAGC‑GGGGGGCGGGGGCGCCCCCGTTGAGCTG‑AACGCCACGCGGGAAAAAAAAA‑TCTTTTAAACCCCCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAaag                                                      <  1:1481878/147‑1 (MQ=255)                                                                                               cGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAaagaa                                                    >  2:207469/1‑151 (MQ=255)                                                                                                  aaagGGAAG‑GGGGGGCGGGGGCGGCCCGTTTGAGCTG‑AAAGCCCCGCGGGAAAAAAAAA‑TCTTCTATACCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGt                                                 <  2:640220/151‑1 (MQ=255)                                                                                                        aaC‑GGGTGGCGGGGGCGGCCCGCTTGAGCTA‑AAAGCCACGCCGGAAAAAAAAA‑TCGTCAATATCTCCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAAcc                                            <  1:867944/150‑1 (MQ=255)                                                                                                         agggggTGGCGGGGGCGGCCCGTTGGAGCTA‑AAGGCCACGCCGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAAcc                                            <  2:31082/147‑1 (MQ=255)                                                                                                          c‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCcg                                          >  1:946773/1‑150 (MQ=255)                                                                                                             ggTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCgg                                       >  2:272471/1‑151 (MQ=255)                                                                                                             ggTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGGCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCgg                                       >  2:1141153/1‑151 (MQ=255)                                                                                                                ggCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAgg                                    >  2:368552/1‑151 (MQ=255)                                                                                                                ggCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAgg                                    >  1:766803/1‑151 (MQ=255)                                                                                                                  cgggggCGGCCCGCTTGGGCTA‑AACGCCCCGCGGGAAAAAAAAA‑TCGTTTATATCCCCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGc                                   <  1:1033735/148‑1 (MQ=255)                                                                                                                        cggcCCGTTTGAGCTG‑AAACCCACGCGGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAAt                             <  2:296782/150‑1 (MQ=255)                                                                                                                        cgccccGCTTGGGCTA‑AAAGCCACGCGGGAAAAAAAAA‑TCGTCAATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAAt                             <  2:248584/147‑1 (MQ=255)                                                                                                                                           acccccgcgGGAAAAAAAAA‑CCTTCTATATCCCCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGAc          <  2:342778/146‑1 (MQ=255)                                                                                                                                           acccccGCCGGAAAAAAAAA‑CCGTCTATACCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGAc          <  2:935040/146‑1 (MQ=255)                                                                                                                                           acccACGCGGAGAAAAAAAA‑GCGTATTTACCCCCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGAc          <  1:901651/149‑1 (MQ=255)                                                                                                                                             ccACGCGGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGAc          <  2:830327/149‑1 (MQ=255)                                                                                                                                             ccACGCGGGAAAAAAAAA‑CCGTCTATATCCCCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGAc          <  2:885845/149‑1 (MQ=255)                                                                                                                                                   cGGAAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTggcgg  <  2:755466/151‑1 (MQ=255)                                                                                                                                                   cGGAAAAAAAAA‑CCGTCTAAATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTggcgg  <  1:949162/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAAC‑GGGTGGCTGGGGCGGCCCGCTTGAGCTG‑AAAGCCACGCCGGGAAAAAAAA‑TCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGACGGATTTAAAGAAGGTGAACCCGCGGAGGCGGAAATTGGTGTCGCGGTAATCGACTGTGGCGG  >  NC_000913/2826315‑2826609

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