Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	40.5%	126.0 / 18.7	85	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (11/39);  new base G (34/0);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 2.79e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508511‑3508806                                                                                                                                                     |                                                                                                                                                         tCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGC‑GC‑GGaaaa                                                                                                                                                       >  1:837758/1‑150 (MQ=255)  ccACGTTGTTGTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGaaaaa                                                                                                                                                      <  2:846999/150‑1 (MQ=255)     cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGgaaaaa                                                                                                                                                     >  1:643127/1‑147 (MQ=255)     cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGgaaaaa                                                                                                                                                     >  2:400698/1‑147 (MQ=255)        tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATtac‑a                                                                                                                                               >  2:780530/1‑150 (MQ=255)          tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATcac‑aca                                                                                                                                             >  1:95155/1‑150 (MQ=255)          tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGC‑GG‑GGAAAAAATcac                                                                                                                                                 >  2:854359/1‑147 (MQ=255)           ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCtc‑tcac                                                                                                                                            <  2:809922/150‑1 (MQ=255)           ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGC‑GC‑GGGAAAAATtac‑aca                                                                                                                                             >  1:355248/1‑149 (MQ=255)               aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTa                                                                                                                                        <  1:1079462/150‑1 (MQ=255)               aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTa                                                                                                                                        <  2:1256283/150‑1 (MQ=255)               aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGC‑GC‑GGGAAAAATAAC‑ACACCGt                                                                                                                                         >  2:298727/1‑149 (MQ=255)                 aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGcc                                                                                                                                     <  1:149180/151‑1 (MQ=255)                 aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCTGCGCGGC‑GGgaaaaacacaaaccgccg                                                                                                                                         >  2:325564/1‑136 (MQ=255)                  gCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGccc                                                                                                                                    <  1:136379/151‑1 (MQ=255)                   ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GCGGGAAAAATTaac‑aacaggagtccc                                                                                                                                    >  2:121251/1‑140 (MQ=255)                   ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGC‑GC‑GGGAAAAAACAC‑AAACAGTCGgccc                                                                                                                                   >  2:1031769/1‑151 (MQ=255)                     ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCca                                                                                                                                 <  1:448068/151‑1 (MQ=255)                     ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGC‑GC‑GGCAAAAATCAT‑ACACgggataacgc                                                                                                                                  >  1:628775/1‑142 (MQ=255)                       agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCca                                                                                                                                 <  1:820456/149‑1 (MQ=255)                         agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATTAC‑ACACCGGAGccccccc                                                                                                                                >  1:467506/1‑146 (MQ=255)                          gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGG‑GC‑CGAAAAAATTAC‑ACCCCGGAGCCCCcac                                                                                                                                >  1:1108984/1‑147 (MQ=255)                             cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATCAC‑ACACCGGAACCCCCCCAACa                                                                                                                            >  1:82960/1‑148 (MQ=255)                             aaaTATGAAGATGCTCACGGGCGAGGGAGTAAGCCGTCGGATCAAAACTCATAGGCTGTCTCCTGCAAACCGGAGGGAGGAAAGACAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑CGAAAAATTtac‑acaaaataccccc                                                                                                                                   >  2:958150/2‑131 (MQ=255)                                  tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATTAC‑ACACGGTAGgccccccccccgaaagacg                                                                                                                    >  1:1205403/1‑137 (MQ=255)                                   gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATCAA‑ACACCGGAGCCCCCACATCACA‑AAACggg                                                                                                                  >  1:328438/1‑149 (MQ=255)                                   gCAGATGCTCATGGGCGAGGCTGTATACCGTCGCATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCACTTTGCGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACAGGGTAGCCCCCACATCAGA‑AAGCgg                                                                                                                   <  2:444751/150‑2 (MQ=255)                                    cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGC‑CC‑GGGAAAATTTAC‑CCACCGGAGGCCCCCCAACAGG‑AAGCgtgg                                                                                                                 >  1:41968/1‑151 (MQ=255)                                       aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGC‑GC‑GGGAAAAATTAA‑AAACCGGAGCCCCCACACCAAA‑AAACgggg                                                                                                                 >  2:1056320/1‑148 (MQ=255)                                        tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGG‑GCGGGAAAAATTAAA‑AAAC‑GTAGCCCCCACCTCAGA‑AAACggggtggg                                                                                                             >  2:198311/1‑149 (MQ=255)                                                   gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GG‑GGGAAAAATTAA‑ACACAGTAGCCCCCCCATCAGA‑AAAGGTGGTGCGCACCCCGAGc                                                                                                   >  1:632141/1‑150 (MQ=255)                                                   gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATTAC‑ACACGGGAGACCCCCCAACACA‑AAACGGGGGGGGcac                                                                                                          >  1:868421/1‑143 (MQ=255)                                                   gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATCAC‑ACACCGTAGCCCCCCCCTCAGA‑AAGCGGGGGGGGCACAcc                                                                                                       >  1:555901/1‑146 (MQ=255)                                                   gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATCAA‑AAACGGGAGGCCCCCCCACACA‑AAAAGTGGGGGGGACACCCAGc                                                                                                   >  2:531267/1‑150 (MQ=255)                                                    aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGAAAAAATAACAACCCGGGCG‑CCCCACAACACA‑AAGCGTGGTGTGCGCAcc                                                                                                       >  2:250340/1‑145 (MQ=255)                                                            aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGG‑GG‑GGGAAAAATCAC‑ACCCGGGAGGCCCCCCATCACA‑AAAACCGGGGCGGACACCGAGCCGCCGAGCg                                                                                          >  2:381131/1‑150 (MQ=255)                                                             cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                                                                        <  1:1068594/151‑1 (MQ=255)                                                             cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                                                                        <  1:217856/151‑1 (MQ=255)                                                                  cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATCAC‑ACACCGTAGCCCCCCCAACAGA‑AAACGTGGTGGGCACACCGAGCCGCGTAGCGATCA‑Tcc                                                                                  >  1:130672/1‑151 (MQ=255)                                                                       cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGc                                                                             <  2:144316/151‑1 (MQ=255)                                                                        aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATCAC‑ACACGGGAGCCCCCCCATCAGA‑AAACGTGGGGCGCACACCGAGCTGCGTAGCGATCC‑Tccc                                                                                 >  1:308019/1‑146 (MQ=255)                                                                         aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATTAC‑ACACCGTAGCCCCCCCACCAGA‑AAACGGGGTGGGCACACCCAGCCGCGTAGCGAACA‑Tcccct                                                                               >  1:840392/1‑145 (MQ=255)                                                                         aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTg                                                                           <  2:869596/151‑1 (MQ=255)                                                                          aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGG‑GC‑GGGAAAAATCAC‑ACAAAGGAACCCCCCCCTCCGA‑AAACGAGGTGGGGACACCGAGCTGCGTAGCGAGCA‑ACCCTTGc                                                                             >  2:1276045/1‑148 (MQ=255)                                                                              cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATaa                                                                       <  1:1227065/150‑1 (MQ=255)                                                                              cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATaa                                                                       <  2:996842/150‑1 (MQ=255)                                                                              cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGG‑GC‑CGGAAAAATCAC‑ACACGGGAGGCCCCCCCCCAGA‑AAAAGGGGGGCGCACACCGAAGGGCGGAGGGGGCA‑TCCCCTGCTGGTaa                                                                       >  1:47200/1‑150 (MQ=255)                                                                                tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGC‑GC‑GGGGAAAATCAC‑ACACCGGAGCCCCCCCATCAGA‑AAACGGGGGGGGCACACCGAGCTGCGGAGCGCTCA‑TCCCCTGCTGATAAAcc                                                                    >  1:921133/1‑151 (MQ=255)                                                                                 aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATCAC‑ACACCGGAGCCCCCCCATCAGA‑AAGCGGGGGGCGCCCACCCAGCTGCGTAGCGGACC‑TCCCCTGCTGATAAAcc                                                                    >  1:238447/1‑150 (MQ=255)                                                                                 aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATCAC‑ACACCGGAGCCCCCCCATCAGA‑AAACGGGGGGCGCACACCGAGCCGCGGAGCGCTCA‑TCCCCTTCTGATaaaa                                                                     >  1:140052/1‑148 (MQ=255)                                                                                 aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGC‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCa                                                                   <  2:410452/151‑1 (MQ=255)                                                                                 aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCa                                                                   <  1:239086/151‑1 (MQ=255)                                                                                 aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGGGC‑GG‑GGAAAAAATCAC‑ACACGGGAGCCCCCCCATACGAGAAGAG‑GGGGGGGACACCGCGCTGCGGAGGGGTCA‑TCCCGTGCTGaaaa                                                                       >  2:506447/1‑147 (MQ=255)                                                                                  gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATCAC‑AAAACGTAGCCCCCACAACAGA‑AAAGGGGGGGGGCCCACCCAGCTGCGGAGCGATCA‑TCCCCTGCTGATaaa                                                                      >  1:838256/1‑147 (MQ=255)                                                                                   ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATTAA‑AAACGGTAGCCCCCCCCTCAGA‑AAACGGGGGGGGGACACCCAGCTGCGTAGCGATCA‑TCCCCTTCTGTTAAACCa                                                                   >  1:1220896/1‑149 (MQ=255)                                                                                      gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCAc                                                               <  1:854390/150‑1 (MQ=255)                                                                                      gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGC‑GG‑GGGAAAATTTTG‑ACAAGGGAACCCCCACCTCACA‑AAAAGCGGTGCGGAGCCCGAGCAGCGGGGCCAACA‑TCCCTTGCTTGCAAACCAGc                                                                 >  1:109397/1‑148 (MQ=255)                                                                                        ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATAAC‑ACACAGGAACCCCCACCACAGA‑AAAAGGGGGGGGCAACCCCAGCCGCGGAGCGATAA‑TCCCCTGCTTATAAAACAGCACttg                                                            >  2:579677/1‑149 (MQ=255)                                                                                                aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACtt                                                    <  2:578439/151‑1 (MQ=255)                                                                                                 aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                   <  1:1082601/151‑1 (MQ=255)                                                                                                    ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                                <  1:687657/151‑1 (MQ=255)                                                                                                      agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATCAA‑ACACCGTAACCCCCCCATCAGA‑AAGCGGGGGGGGCACACCCAGCTGCGGAGCGATCA‑TCCCCTTCTTATAAAACAACACCGGCACCACCTCCCggg                                              >  1:1165285/1‑151 (MQ=255)                                                                                                      agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                              <  2:26859/151‑1 (MQ=255)                                                                                                      agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                              <  2:625876/151‑1 (MQ=255)                                                                                                      agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGC‑GG‑GGGAAAAATAAC‑ACACCGGAACCCCCACCTCAGA‑AAACGGGGGGCGCCCACCGAGCTGCGGAGAGGTCA‑TCCCTTTCTGGTAAAAAACCACCGGCACCACTTCGCggg                                              >  2:1034310/1‑151 (MQ=255)                                                                                                       gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                              <  2:446997/150‑1 (MQ=255)                                                                                                       gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                             <  1:601951/151‑1 (MQ=255)                                                                                                            gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggt                                        <  1:530125/151‑1 (MQ=255)                                                                                                             aaagaaaagaGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAAACAC‑ACACGGTAGCCCCCCCATCCGA‑AAACGGGGGGGGCACACCCAGCCGGGGGAGGAGCA‑TCCCCTTCTTGTAAAACAACACTGGCACCCCTTCCCGGGGAtgggg                                       >  2:247897/1‑149 (MQ=255)                                                                                                              aagaaaagaGCGACCATTTTGTGGCGGGGGCGC‑GC‑GGGAAAAATAAA‑ACACCGGAGCCCCCACATCAGA‑AAGCGGGGGGCGCACCCCGAGCTGCGTAGCGATCA‑TCCCTTTCTTATAAAACAGCACCGGGACCACCTCCCGGGGGtggtgg                                      >  2:686552/1‑151 (MQ=255)                                                                                                               agaaaagaGCGACCATTTTGTGGCGGAGGCGC‑GC‑GGAAAAATTAAA‑AAACGGGTACCCCCCCAACAGC‑AAGAGCGGGGGGCACACCGCGCCCCGGAGGGGTCATTCCCTTCTTGTAAAACCA‑CACCGGGACCACCTCGGGGGtgttg                                         >  2:57126/1‑143 (MQ=255)                                                                                                                gaaaagaGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATAAC‑AAACGGCAGCCCCCCCATCAGA‑AAACGGGGGGGGCACACCCCGCCGCGGAGCGAGCA‑TCCCTTGCTGGTAAAACAACACCGGCACCCCATCCCGGGGAtggtgg                                      >  2:405569/1‑149 (MQ=255)                                                                                                                    agagCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGa                                <  1:780558/151‑1 (MQ=255)                                                                                                                      agCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                              <  1:1132264/151‑1 (MQ=255)                                                                                                                       gCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                             <  1:1028995/151‑1 (MQ=255)                                                                                                                        cGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAATTCAC‑ACACGGTAGCCCCCCCATCAGA‑AAACGGGGGGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGATGGTGGGTGTGAcc                              >  1:29383/1‑149 (MQ=255)                                                                                                                         gACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAAt                           <  2:767892/151‑1 (MQ=255)                                                                                                                          aCCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                          <  1:198319/151‑1 (MQ=255)                                                                                                                             aTTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                       <  1:1082972/151‑1 (MQ=255)                                                                                                                               tttGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                      <  1:1159495/150‑1 (MQ=255)                                                                                                                               tttGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                      <  2:921095/150‑1 (MQ=255)                                                                                                                                 tgtgGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGa                    <  2:411425/150‑1 (MQ=255)                                                                                                                                     gcggtggcgC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt               <  1:1135350/151‑1 (MQ=255)                                                                                                                                     gcggtggcgC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt               <  1:1019408/151‑1 (MQ=255)                                                                                                                                      cggtggcgC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGCAGCGTCCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                 >  1:411088/1‑148 (MQ=255)                                                                                                                                          ggcgC‑GC‑GGAAAACTTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTCGCTGTGACCAATGGTCGGATCGTTGCCGt          >  2:429092/1‑151 (MQ=255)                                                                                                                                          ggcgC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg           <  1:29894/150‑1 (MQ=255)                                                                                                                                          ggcgC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt          <  1:1225080/151‑1 (MQ=255)                                                                                                                                           gcgc‑gc‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg         <  2:856242/151‑1 (MQ=255)                                                                                                                                            cgc‑gc‑gGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg         >  2:354857/1‑150 (MQ=255)                                                                                                                                                     aaaaaTTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc  >  2:208705/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                         TCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCAC‑ACACGGTAGCCCCCACATCAGA‑AAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCA‑TCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508511‑3508806

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