Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	45.8%	83.2 / 16.3	60	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (2/30);  new base G (27/0);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 8.39e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCG  >  NC_000913/3508527‑3508797                                                                                                                                    |                                                                                                                                              aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGGAAAAATAACA‑CACGGTAGcc                                                                                                                          >  2:75412/1‑151 (MQ=255)     ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCca                                                                                                                      <  1:1054313/151‑1 (MQ=255)      gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGTGCGCGGGAAAAATccgac                                                                                                                                    >  1:962487/1‑133 (MQ=255)      gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGGGCGGGAAAAATCACACCCCGGAAGcccccc                                                                                                                      >  1:289394/1‑150 (MQ=255)          gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCa                                                                                                                 <  1:578390/151‑1 (MQ=255)           aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAg                                                                                                                <  1:758198/151‑1 (MQ=255)            tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGCGGGGAAAATTACA‑CACGGTAGacccc                                                                                                                       >  1:514022/1‑143 (MQ=255)             cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGCGGGAAAAATCACA‑CACGGGAGACCCcacacccga                                                                                                               >  1:592485/1‑145 (MQ=255)             cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGaa                                                                                                              <  2:1129356/151‑1 (MQ=255)             cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGa                                                                                                               <  1:713032/150‑1 (MQ=255)                  tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGGAAAAATCACA‑CACCGGAGCCCCCCCACCAGaaa                                                                                                             >  1:512040/1‑147 (MQ=255)                   gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑TACGGTAGCCCCCACATCAGAAAGCgtg                                                                                                        <  1:865270/151‑1 (MQ=255)                     aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGGGGGAAAAATTACA‑CACGGTAGCCCCCCCATCCGAAAACgggg                                                                                                       >  2:1158807/1‑150 (MQ=255)                      gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGCGGGAAAAATCACA‑CACCGGAGCCCCCACATCAGAAAgagggg                                                                                                       >  1:1093191/1‑144 (MQ=255)                         gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGGAAAAATCACA‑AAAGGGAGCCCCCCCAACAGAAAGCgggggggg                                                                                                   >  1:637949/1‑146 (MQ=255)                              tGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGCGGGAAAAATCACA‑AACGGTAGCCCCCCCATCAGAAAAAGGGGGGCGCACAcc                                                                                             >  2:1120660/1‑151 (MQ=255)                               gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCg                                                                                            <  1:480734/151‑1 (MQ=255)                                 gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAg                                                                                          <  2:1190195/151‑1 (MQ=255)                                         gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                   <  1:80126/150‑1 (MQ=255)                                                  cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGGAAAATTCACA‑AACGGTAGCCCCCACATCAGAAAACGGGGGGGGCACACCCAGCCGCGGAGCGGTCATcc                                                                         >  2:864651/1‑151 (MQ=255)                                                  cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                                                         <  2:1039581/151‑1 (MQ=255)                                                   ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                                                        <  1:66319/151‑1 (MQ=255)                                                   ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                                                        <  1:281534/151‑1 (MQ=255)                                                       cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGCGGCAAAAATTAAA‑CACGGTAGGCCCCCCCACAGAAAGCGAGGGGGGCACACCCAGCTGCGGAGCGATCATcccct                                                                      >  2:890277/1‑147 (MQ=255)                                                        aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                                                                   <  1:628109/151‑1 (MQ=255)                                                           aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGGAAAAATCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCTggt                                                                >  2:329512/1‑149 (MQ=255)                                                            cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGCGGGAAAAATTAAA‑CACGGGAAGCCCCACAACAGGAAAAGGGGGGCGGCCACCCAGCCGCGGAGGGGTCATccccttc                                                                    >  1:703180/1‑142 (MQ=255)                                                               aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGGAAAAATTACA‑CACGGTAGCCCCCCCCTCAGAAAGCGGGGGGGGCACACCCCGCTGCCGAGCGATCAACCCCTGCTGATaaa                                                             >  1:40091/1‑150 (MQ=255)                                                                 aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                           <  2:430370/150‑1 (MQ=255)                                                                         tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCTGGGCGGCGGGAAAA‑TCACA‑AACAGCAGACCCCCCACCACAAAGAGGGGGGCGCCCACCCCGGCGCGTAGAGCGCATCCCTTGTTGTTAAAACaccccc                                                     >  1:1168961/1‑143 (MQ=255)                                                                          ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                                                    <  1:1119112/148‑1 (MQ=255)                                                                               aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGCGGGAAAAATCACA‑CAAGGTAGCCCCCCCATCAGAAAACGGGGGGCGCACAACCCGCTGCGTAGCGATCATTCCCTTCTGGTAAAACAGCACCgggccc                                               >  1:1166107/1‑144 (MQ=255)                                                                               aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACt                                            <  1:615306/151‑1 (MQ=255)                                                                                aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGCGGGAAAAATCACA‑AACGGTGCGCCCCCCATCAGAAAACGTGGGGGGCACACCCGGCTGCGGAGGGATCATCCCTTTCTGATAAAACCGCACTGGcacc                                               >  2:50217/1‑147 (MQ=255)                                                                                 aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGCGGGAAAAATCACA‑CACCGGAGACCCCACATCAGAAAGCGTGGGGGGCACACCGAGCTGCGGAGCGATCATCCCTTTCTGATAAACCAGCAACGGCACCACCTc                                          >  1:1175731/1‑151 (MQ=255)                                                                                 aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                          <  2:607672/151‑1 (MQ=255)                                                                                 aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                          <  2:117823/151‑1 (MQ=255)                                                                                 aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                          <  2:242666/151‑1 (MQ=255)                                                                                  ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGGAAAAATTACA‑AACGGTAGGCCCCCCAACAAAAAACGGGGGGCGCCCACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAAACCACACTGGCCCCACTTc                                          >  2:686760/1‑150 (MQ=255)                                                                                   cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgc                                        <  1:954018/151‑1 (MQ=255)                                                                                    ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGCGCGGGAAAAATTACA‑CACCGGAGCCCCCCCATCAAAAAACGGGGGGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGGTAAAACAACACCGGCACCACCTcgc                                        >  2:664780/1‑150 (MQ=255)                                                                                       gggaggAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                    <  2:26515/151‑1 (MQ=255)                                                                                        ggaggaAAGAAAAGAGCGACCATTTTGTGGCGG‑GGGCGCGGGGGAAAAATTAAA‑CACGGGAGACCCCACAACAGAAAGCGGGGGGCGCGCACCCAGCCGCGGAGCGAACATCCCCTGCTGGTAAAACAGCACCGGCACCACCTCCCggg                                     >  1:99380/1‑149 (MQ=255)                                                                                           ggaAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg                                <  2:456673/151‑1 (MQ=255)                                                                                           ggaAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg                                <  2:590126/151‑1 (MQ=255)                                                                                            gaaagaaaAGAGCGACCATTTTGTGGCGG‑TGGTGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggt                               <  2:1013870/151‑1 (MQ=255)                                                                                                gaaaagaGCGACCATTTTGTGGCGG‑TGGCGCGCGGGAAAAATTACA‑CACGGGAGCCCCCCCATCAGAAAACGTGGGGCGCACACCGAGCTGCGGAGCGATCATCCCTTGCTGATAAACCAACACTGGCACCACTTCGCGGGGAtggtgg                             >  2:46576/1‑149 (MQ=255)                                                                                                  aaagaGCGACCATTTTGTGGCGG‑TGGCGCGCGGGAAAAATCACA‑AACGGTAGCCCCCCCCTCCGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACCGGCACCACCTCGCGGGGAtggggg                             >  2:12079/1‑147 (MQ=255)                                                                                                  aaagaGCGACCATTTTGTGGCGG‑GGGCGCGCGGAAAAAATCACA‑CCCGGTAGGCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCTGCGTAGCGGTCATCCCCTTCTCATAAAACACGACCGGCACCCCCTCCCGGGGAtggtgg                             >  2:255198/1‑147 (MQ=255)                                                                                                     gagCGACCATTTTGTGGCGG‑TTGCGCGCGGGAAAAATCAAA‑AACAGGAAGCCCCACATCATAAAGCGCGGGGCGCACAACCAGGTGCGCGGCGATCATACCTTGCTGATAAAAAAGCACAGGCCCCCCCTTGCGGGGATGGTGGCtgtg                        >  2:906813/1‑149 (MQ=255)                                                                                                     gagCGACCATTTTGTGGCGG‑TGGCGGGGGGGAAAATTTACC‑ACCGGGAGCCCCCCCCTTAGGAAGAGGGGTGCGGCCACCCAGCCGCGGAGCGGTCAACCCCTCCTGGTAAAACAGCACTGGCACCCCTTCGCGGGGATGGTGGGtgtgt                       >  1:500662/1‑149 (MQ=255)                                                                                                     gagCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAc                      <  2:998418/151‑1 (MQ=255)                                                                                                         gACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAAt                  <  1:75412/151‑1 (MQ=255)                                                                                                          aCCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                 <  2:63231/151‑1 (MQ=255)                                                                                                           ccATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATgg                <  1:522481/151‑1 (MQ=255)                                                                                                            cATTTTGTGGCGG‑TGGCGCGCGGGAAAAATTACA‑CACGGGAGCCCCCCCATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAAACAGCACTGGCACCACCTCCCGGGGATGGGGGGTGTGACCAATgg                >  2:685819/1‑150 (MQ=255)                                                                                                             aTTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc              <  2:592500/151‑1 (MQ=255)                                                                                                               tttGTGGCGG‑TGGCGCGCGGGAAAAATTACA‑CACCGGAGCCCCCCCATCAGAAAACGGGGGGGGCACACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAAACAACACTGGCACCACATCGCGGGGATGGGGGGGGTGACCAAAGGGCgg            >  1:523787/1‑151 (MQ=255)                                                                                                                     gcgg‑tggcgCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt      >  2:71621/1‑151 (MQ=255)                                                                                                                           ggcgCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg  <  1:415840/150‑1 (MQ=255)                                                                                                                                    |                                                                                                                                              AGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGCGGAAAAATTCACA‑CACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCG  >  NC_000913/3508527‑3508797

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: ‑