Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,689	0	G	T	38.7%	49.8 / 20.3	63	intergenic (+43/+22)	yegJ/yegK	uncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (25/13);  new base T (10/14);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 7.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-02

GTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCC  >  NC_000913/2151541‑2151834                                                                                                                                                         |                                                                                                                                                   gTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTc                                                                                                                                                  <  1:23121/151‑1 (MQ=255)        cgcgGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTCAAcc                                                                                                                                           <  1:555354/151‑1 (MQ=255)                            cggTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTTTTTTCAACCCCCATGAGAAaatcac                                                                                                                           >  2:1305732/1‑147 (MQ=255)                               ttatccttttGTTCTTTTATGGGTCCCCTATCAAAAAGAAGAAGCCGACGATCTAAGAGAAAATTTATTTCGTTATTCAGCCTTTATTTTTTTTTTTTAAC‑CTGGTTT‑TTAAT‑TT‑GTTGTCTTAAACCGCCATGAGAAACTCActc‑t                                                                                                                       <  2:826782/139‑1 (MQ=255)                                tataCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTTTCTTCAACCGCCATGAGAAACTCACTC‑TGcc                                                                                                                    >  2:600334/1‑149 (MQ=255)                                  taCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑ATGAA‑TT‑GTTTTGTTCTACCGCCCAGAGAAAGTCACCC‑TTccccc                                                                                                                 >  1:608930/1‑148 (MQ=255)                                   aCCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑GTGAA‑AT‑GTTGTTCTCAACCCCCATGAGAAACTAACTC‑TAccc                                                                                                                   >  1:281679/1‑147 (MQ=255)                                    cttttGTTCCTTTTTTGGTCACCTACCAAAAAAAAGAGGCCGACGAGCTAAGAGAAAAATTATTTCGCTAATCGCCCGTTATTTTTA‑TGTTTAAACCGGGCGG‑TTAAA‑TT‑TTTGTCTAAAACCGCCATGAGAAACTCACTC‑TGCCca                                                                                                                  <  2:107841/144‑1 (MQ=255)                                    ccATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTTAA‑TT‑GTTGTCTTCAACCCCCCTGAGAAACTCACTC‑TGccccc                                                                                                                 >  1:823411/1‑146 (MQ=255)                                       ttGATCCTTTTTGTATCACCTATCCAAAAGAAGAACCCGAGGATCAAAGAAAAAAATTAGTTCGCTAATCAGCCTTTATTTTTTTTGTTAAAC‑CGGGGGG‑TGTTA‑TT‑TTTTTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAgc                                                                                                              <  1:1010049/148‑1 (MQ=255)                                        tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTTAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGcc                                                                                                           >  1:45296/1‑150 (MQ=255)                                         gATCTTTTATTGGTCACTTACCAAAAAGAAGAAGCCGACAAGCAAAGAAAAAAATTATTACGAAAATCACCTTTTATATTTATTTTTTAAC‑CTGGCGG‑TTAAT‑TT‑GTTTTCTCAAGCCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCa                                                                                                          <  2:1032707/150‑1 (MQ=255)                                            ccTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTTTCTTCAACCGCCATGAGAAAATTACTC‑TGCCCCCAGCGCCAGTg                                                                                                       >  2:1168700/1‑150 (MQ=255)                                                 aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTCAACCGCCATGAGAAAATCACTC‑TGCCCACAGCGCCAGTGCCaa                                                                                                   >  1:1107555/1‑149 (MQ=255)                                                  ttGGTCCCCTATCAAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTTTTTCGCTAATCACCCTTTATTTTTATTGTTAAAC‑CGGGCGG‑TTGAA‑TT‑GTTGTCTCCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAAc                                                                                                  <  2:210232/149‑1 (MQ=255)                                                    ggTCCCCTATCAAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTAGTTCGATAATCAGCCTTTATATTTTTTTTTAAAA‑CTGGCGT‑TTAAA‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACg                                                                                                 <  1:775760/148‑1 (MQ=255)                                                       ccccTATCAAAAAAAAAAACTCAACGATCTAGAAGAAAAATATTTTCTCTACGCAGCTTCTATATTT‑TTTTTTGAA‑CGTGGGG‑TTTAA‑TTGGTGGTCTTCAGCCGCCATGAGAAGCTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTat                                                                                           <  1:1087495/149‑1 (MQ=255)                                                       cACCTAACAAAAAAAAAAAGCCGACGAGAAAAGACAAAATTTATTTCCCTACCCACCCTTTATATTTATTGTTAAAA‑CGGGCGT‑TTAAA‑TT‑GTTTTCTCCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGtt                                                                                               <  2:682139/147‑1 (MQ=255)                                                           tATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTTAA‑TT‑GTTTTCTTCAAACCCCCTGTGAAACTCACTC‑TGCCCCCACCGCCACTGCCAACGTTTTat                                                                                           >  1:515229/1‑147 (MQ=255)                                                              ccAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑GTGAA‑TT‑TTTGTTTTCAAACCCCATGAGAAAATCACTC‑TGCCCACAGCGCCCGTGCCAACGTTTTatcatca                                                                                      >  1:1021886/1‑149 (MQ=255)                                                                aaaaaaAGAAGCCGACGAGCTAAGAACAAATTTATTTCTCTAATCAGCCTTTATATTTTTTGTTTAAC‑CGGGGGG‑TTGAT‑TT‑TTTTTCTTAAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATcc                                                                                   <  1:1290287/150‑1 (MQ=255)                                                                 aaaGAAAAAGCCGAGAAACAAAGACAAAAATTAGTCCTCTAATCAGCCTTTTTTTTTTTTGTTAAAC‑CTGGCGG‑TTGAA‑TT‑GTTGTCTCCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGt                                                                                 <  2:560548/151‑1 (MQ=255)                                                                  aagaagaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGTA‑TT‑GTTTGCTTCAAACGCCATGAGAAAATCCCCC‑TGCCCACAGCGCCAGTGCCAACGTTTTATGATCATTCGt                                                                                 >  2:409391/1‑150 (MQ=255)                                                                   agaagacaCAGACGAAATAAGAGAAAAGTTAGTCCAATGAGCAGCCTTTATTTTTTTTGTAAAAA‑CGGGGGT‑TTTAA‑TT‑GTTTTTTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAc                                                                               <  2:114533/151‑1 (MQ=255)                                                                     aagaagCCGTGGATATAAGAAAAAAATTATTTCCTTAATCACCTGTTATTTTTATTTTTAAAA‑CTGGGGG‑TTAAT‑TT‑TTTTTCTACACCCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTAc                                                                               <  1:822101/149‑1 (MQ=255)                                                                             gacgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTTTCTTCAACCGCCATGAGAAAATCACTC‑TGCCCACAGCGCCAGTGCCCACGCTTTATCATCCTCCCTTCGCCCGtttt                                                                      >  2:81047/1‑149 (MQ=255)                                                                             gacgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTTAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGGGCCAACGTTTTATCATCATCCGTACGCTCGTTTa                                                                      >  1:1168101/1‑150 (MQ=255)                                                                                     aagaaaaaaaTTATTCCCCTAACCGGCTTTTATTTTTTTTTTTAAAC‑CGGGGGG‑TTAAT‑TT‑GTTTTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAccgcc                                                                 <  2:251528/141‑1 (MQ=255)                                                                                     aaGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCggg                                                              >  1:675328/1‑150 (MQ=255)                                                                                        agCAAAATTAGTTCGCTCATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑CTGGA‑TT‑GTTGTTTTCCATCGCCATGAGAAAATCCCTC‑TTCCCACAGGGCCACCGGCAACGGTTTTTAATCATCCCCTCGCTCGTTTACCCCCGGGGTgg                                                          >  1:297745/1‑151 (MQ=255)                                                                                           aaatttttctCGCAAATCAGCTTTTAT‑TTTTTTTTTTAAA‑CCGGGGGTTTTAT‑TT‑TTTTTCTAAAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACa                                                       <  1:1048207/142‑1 (MQ=255)                                                                                           aaaaTTAGTCCGATAACCCCCTTTTAT‑TTTTTTTTTTAAA‑CCGGGGGTTTGAT‑TT‑TTTTTCTCCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACa                                                       <  2:999562/151‑1 (MQ=255)                                                                                            aaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTTAA‑TT‑TTTGTCTTTAACCGCCATTAGAAAATCACTC‑TGCCCACAGCGCCAGGGCCCACCTTTTTTCATAATCCGTACGCTCGTTTTCCCCCCGGc                                                             >  2:213115/1‑144 (MQ=255)                                                                                             aattggtccgtaaaTCAGCCTTTATATTTATTG‑TTAAA‑CCGGGGG‑TTTAATTT‑GTTTTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaa                                                     <  1:171300/139‑1 (MQ=255)                                                                                                 attccataaaTCAGCTTTTTTTTTTTTTGTTAAAC‑CTGGTGT‑TTAAT‑TT‑GTTTTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaaaa                                                   <  2:3851/141‑1 (MQ=255)                                                                                                   ttCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑GTTGA‑TT‑CTTGGCTTCACACCCCATTAAAAAATTTCTCGTGCCC‑CAGCCCCCCGGCGAACGTTTTTTAATTATCCCTACGCCCCTCTACCTCCGGGGGGGACcaaaaa                                                  >  2:866157/1‑147 (MQ=255)                                                                                                    tCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGGCTTCCAATTCCCTGAGAAACTCCCCC‑TGCCCACAGCGCCCGTGCCAAAGCTTTTTTAAGATCCGTACGCCCGTTTAACGCCCGGCTGGACaaaaagtcct                                              >  1:850615/1‑146 (MQ=255)                                                                                                       cTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑CTGAA‑TT‑GTTGTCTTCAACCGCCATGAGAAAATCACTC‑TGCCCACAGCGCCCGTGCCCAAGTTTTATCATTATCCGTACGCTCGGTTACCGCCGGGCTGGGCAAAAAATGCttt                                            >  1:445071/1‑150 (MQ=255)                                                                                                          accAGCCTTTAT‑TTTTTTTTTTAAA‑GTGGGGG‑TTTAATTT‑GTTTTTTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACa                                        <  1:684466/149‑1 (MQ=255)                                                                                                          aTCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTCAACCGCCATTAGAAAACCACTC‑TGCCCACAGCGCCAGCGGCAACGTTTTATCATTATCCGTACGCTCGTTTTCCCCCCGGCCTGACAAAAAATGCTTTa                                           >  1:218150/1‑148 (MQ=255)                                                                                                                 gTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑AT‑TTTGTCTTCAACCCCCCTTAGAAAATCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCCTACGCCCGTTTACCGCCGGGCCGGACAAAAACTGTTTTAACAAtt                                     >  2:906118/1‑147 (MQ=255)                                                                                                                    atatTTATTGTTTAAA‑CTGGCGG‑TTTGA‑TT‑GTTTTCTTTAACCCCCATGAGAAAATCACTC‑TGCCCACAACGCCAGGGCCCACGTTTTTTCATCATCCCTACGCTCGTTTTCCCCCCGGCTGGACAAAAAATGCTTTAAAAATTCAgg                                 >  2:60651/1‑148 (MQ=255)                                                                                                                    atatTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑TTTGTCTTCAACCGCCATGAGAAAATCACTC‑TGCCCACAGCGCCAGTGCCAACGGTTTATCATCATCCGTACGCTCGCTTACCGCCGGGCTGGACAAAAAATGCTTTAACAATTCAgg                                 >  1:504400/1‑148 (MQ=255)                                                                                                                     tatttattGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTCAACCGCCCTGAGAAACTTACTC‑TTCCCCCAGCGCCAGTGCCCACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAAATCAGGtaa                              >  1:332346/1‑150 (MQ=255)                                                                                                                       ttttttttttAAA‑CGGGGGG‑TTAAT‑TT‑GTTGTCTCCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGta                               <  1:148593/147‑1 (MQ=255)                                                                                                                       tttattGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTTTCTTCAAACGCCATGAGAAACTCAATC‑TGCCCACAGCGCCAGTGCCAAAGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCCGGACAAAAACTGCTTTAACAATTCAGGtaa                              >  1:158001/1‑148 (MQ=255)                                                                                                                       tttattGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTTAACCCCCCTGAGAAACTCACTC‑TGCCCACCGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGACAAAAACTTCTTTAACAATTCAGGtaata                            >  1:25857/1‑150 (MQ=255)                                                                                                                           ttGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTTTCTTCAAACCCCATGAGAAAACCACTC‑TGCCCCAAACGCCAGTGCCAACGGTTTTTTATCATCCCCACGCTCCTTTAACCAACGGCCGGACAAAAAACCCTTTAACAATTCTGGgaataa                           >  2:55549/1‑147 (MQ=255)                                                                                                                           ttGTTAAAA‑CGGGCGT‑TTAAT‑TT‑GTTGGCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATaaa                          <  1:733316/148‑1 (MQ=255)                                                                                                                              tttAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTCAACCCCCATGATAAACTCACTC‑CGCCCACCACGCCAGTGCCAACGTTTTATTATCATCCCTACCCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCCGGTAATAAAACGa                      >  2:205692/1‑149 (MQ=255)                                                                                                                              tttAAA‑CGGGGGT‑TTGAT‑TT‑TTTGTCTCAACCCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATc                        <  1:1265597/147‑1 (MQ=255)                                                                                                                              tttAAA‑CGGGCGT‑TTGAA‑TT‑TTTGTCTTCAGCCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAg                     <  2:497517/150‑1 (MQ=255)                                                                                                                               ttaaaaggggCGT‑TTAAT‑TT‑TTTTTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGa                      <  1:1323421/141‑1 (MQ=255)                                                                                                                                  acgcTGGCGG‑TAAAA‑TT‑ATTATATCCCACCGCCATGAGAAACTCACTC‑TGCCCGCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgtt                   <  2:1302036/146‑1 (MQ=255)                                                                                                                                  aacggggggg‑TTGAT‑TT‑TTTTTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgtt                   <  1:1067228/141‑1 (MQ=255)                                                                                                                                  aa‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgtt                >  1:925676/1‑151 (MQ=255)                                                                                                                                  aa‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgtt                >  1:364114/1‑151 (MQ=255)                                                                                                                                   aatctgggg‑TTGAATTT‑TTTGTCTCAAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgt                 <  2:263822/144‑1 (MQ=255)                                                                                                                                         gtggtttaattt‑gtttTCTCCACCCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTc               <  2:379998/138‑1 (MQ=255)                                                                                                                                          gg‑TTGAA‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGt         >  2:726344/1‑151 (MQ=255)                                                                                                                                           g‑TTGAA‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGtt        >  2:62387/1‑151 (MQ=255)                                                                                                                                             ttGAA‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtg       >  1:935244/1‑151 (MQ=255)                                                                                                                                                aa‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgctg    >  2:659595/1‑151 (MQ=255)                                                                                                                                                 a‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgctgc   >  1:492140/1‑151 (MQ=255)                                                                                                                                                    t‑gttttCTCCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTgcc  <  2:94325/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                   GTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAA‑CTGGCGG‑TTGAA‑TT‑GTTGTCTTCAACCGCCATGAGAAACTCACTC‑TGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCC  >  NC_000913/2151541‑2151834

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