Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	345,328	0	C	A	46.0%	91.7 / 43.1	63	intergenic (+337/‑76)	yahL/yahM	uncharacterized protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (33/1);  new base A (0/29);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 3.95e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCC  >  NC_000913/345179‑345477                                                                                                                                                        |                                                                                                                                                      aaaaaaTGATCTCATGCAGATGTTTTGTGAATGAGTTGGTTGACATTTATATGAAAAAATTCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGGTCATATc                                                                                                                                                       >  2:1250588/1‑150 (MQ=255)    aaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCaaa                                                                                                                                                    >  1:729602/1‑150 (MQ=255)    aaaTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCaaa                                                                                                                                                    >  2:813495/1‑150 (MQ=255)      aTGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAAc                                                                                                                                                 >  2:1028220/1‑151 (MQ=255)                  aGAAGTTTTGTGAATGTGTTGGTAGACATTCATATGAAAACAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑GGTATCTTACGCATCGCAGTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCACAAACCGCCATTAGCtc                                                                                                                                     >  2:786732/1‑151 (MQ=255)                     tGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtcatc                                                                                                                                  >  1:362997/1‑151 (MQ=255)                      gTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCtcatc                                                                                                                                  >  1:1153419/1‑150 (MQ=255)                          tgtgAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGa                                                                                                                              >  1:24851/1‑150 (MQ=255)                               aTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGc                                                                                                                         >  2:558196/1‑150 (MQ=255)                                 gtgtTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGc                                                                                                                         >  2:1033908/1‑148 (MQ=255)                                  ttttggttgACATTAATAAAAAAAAAAACAAAATCCCATCATGTTTGTGGGGGTTCTTTTC‑TGTATCTTACCCCCCCCACTAAACCCTA‑AAAAAAAATATGTGGAAGGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCa                                                                                                                        <  2:1005607/146‑1 (MQ=255)                                    tggtttgACATTTATAAAAAAAAAATCATATTCCCATCTTTTTTGTGGGGGTTTTTTTT‑TTTATCTTCCGCACCGCCCTAAACCTTA‑AAAAAAAATATGTTGAAAGCCTCAAATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAga                                                                                                                      <  2:2239/144‑1 (MQ=255)                                         tgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                              >  2:425863/1‑151 (MQ=255)                                         tgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc                                                                                                              >  1:62349/1‑151 (MQ=255)                                            cATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                           >  1:294272/1‑151 (MQ=255)                                             aTTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAgg                                                                                                           >  1:985731/1‑150 (MQ=255)                                              ttaatatAAAAAAAATCATAATCCCATCATTTTTGTGTGGGGTCTTTTT‑TTTATCTTACCCATCGCACCCAACCCTG‑ACAAAAAATATGCTGTAAGCCTAATACAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAg                                                                                                            <  2:934858/145‑1 (MQ=255)                                                 atatGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt                                                                                                      >  2:350770/1‑151 (MQ=255)                                                   ataaaaaaaaTCATAATCCCTTCATTTTTG‑GGGGGGTTTTTTCTTTTATCTAACGCACCGCACCCAAGCCTA‑ACAAAAAATATGCTGAAAGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAt                                                                                                         <  2:59883/143‑1 (MQ=255)                                                      aaaaaaaTTATAATCCCACCTTTTTTGTGTGGTTTCTTTTC‑TGTACTTTACACACCGCA‑AAAAACCTGCAAAAAAAATTTCCTTGAAGGCAAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                                 <  1:725228/151‑1 (MQ=255)                                                      aaaaaaaTAATAATCCCATCTTGTTTGTGGGGGTTTTTTTT‑TGTATCTTACCCACCCCACTCAAGCCTA‑AAAAAAAATATGTTGAAAGGCTCATAACAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt                                                                                                 <  1:350566/151‑1 (MQ=255)                                                        aaaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                               >  2:566877/1‑151 (MQ=255)                                                        aaaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATgg                                                                                               >  1:841940/1‑151 (MQ=255)                                                         aaaaTCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGt                                                                                              >  2:548737/1‑151 (MQ=255)                                                               aTAATCCCATCTTTTTTGTGTGGGTTTTTTTC‑TGTTTCTACCGCACCGCCCTCAAGCTTA‑AAAAAAAATATGCGGTAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACggg                                                                                         <  2:790375/150‑1 (MQ=255)                                                                 aaTTCCATCATTTTTGTGGGGGTTTTTTTT‑TTTATCTTACGCACCCCCCTCAACCCTG‑AAAAAAAATATGCGGAAGGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                                                      <  2:957735/151‑1 (MQ=255)                                                                 aaTTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGtt                                                                                      >  1:4274/1‑151 (MQ=255)                                                                 aaTCCCATCATTTTTGTGTGGGGTTTTTTC‑TGTATCTTCCGCACCGCACTCAACCCTG‑AAAAAAAATATGCTGAAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg                                                                                           <  1:580369/146‑1 (MQ=255)                                                                 aaTCCCATCATGTTTGTGGGGGGTTTTTTT‑TGTTTCTTCCGCATCCCACTCAAGCCTG‑AAAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACgggg                                                                                        <  1:1007855/149‑1 (MQ=255)                                                                    cccATCTTTTTTGTGGGGGGTTTTTTC‑TTTATCTTACCCACCCCCCACAAGCCAA‑AAAAAAAATTTGTTTAAAGCCTCATAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGt                                                                                       <  2:796248/146‑1 (MQ=255)                                                                     ccatattGTTGGTGGGGGGTTTTTTT‑TTTATCTTACGCACCGCACTCAACCCTG‑AAAAAAAATTTGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTc                                                                                     <  2:976566/142‑1 (MQ=255)                                                                       atcatGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGc                                                                                >  1:442139/1‑151 (MQ=255)                                                                        ttttttttGTGTGGGTTTTTTTC‑TTTATCTTCCGCATCCCACCCAAGCCTG‑AAAAAAAATATGTTGAAGGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTc                                                                                     <  2:1004318/140‑1 (MQ=255)                                                                          atGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTc                                                                             >  1:472412/1‑151 (MQ=255)                                                                           tGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCg                                                                            >  1:43332/1‑151 (MQ=255)                                                                             ttggggtgGGTTTTTTTT‑TGTTTCTTCCGCTCCACACTCAAGCAAG‑AAAGAAAATAGGCTGTAGGGTTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg                                                                                 <  1:143858/139‑1 (MQ=255)                                                                                  gtgGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAACTATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCgg                                                                     >  2:1081621/1‑151 (MQ=255)                                                                                   tgGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGt                                                                    >  1:567563/1‑151 (MQ=255)                                                                                      gtttttttt‑tGTATCTCACGCATCGCACTCAACCCTG‑AAAGAAAATATGCTGTAGGCCTAAAATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTcc                                                                  <  1:73547/146‑1 (MQ=255)                                                                                        tCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAAc                                                               >  1:998624/1‑151 (MQ=255)                                                                                            ttC‑TGTATCTTCCGCTCCGCACCCACCCCTA‑AAAAAAAATATGCTGAATGGTACATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACttt                                                            <  2:985486/150‑1 (MQ=255)                                                                                             tC‑TTTATCTTACGAACCGCACACAGCCCGG‑ACAAAAAATATGCTGTAAGCCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtttt                                                           <  2:247923/150‑1 (MQ=255)                                                                                              tttGTATCTTACCCACCCCACTCAAGCCTG‑AAAAAAAATATGCTGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtttt                                                           <  1:946962/148‑1 (MQ=255)                                                                                                  tATCTTACGCATCGCACCCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaa                                                       >  1:829975/1‑150 (MQ=255)                                                                                                   aTCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaaa                                                      >  2:989390/1‑150 (MQ=255)                                                                                                       tACGCATCCCACTCAACCCTG‑AAAAAAAATATGCTGTAAGGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAc                                                     <  1:975224/147‑1 (MQ=255)                                                                                                           cATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTg                                              >  1:629206/1‑150 (MQ=255)                                                                                                              cGCACTCAACCCTA‑AAAAAAAATATGCTGTAAGGCTCAAAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa                                             <  2:665302/148‑1 (MQ=255)                                                                                                                  cTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                          >  1:451021/1‑147 (MQ=255)                                                                                                                   tCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                          >  1:532543/1‑146 (MQ=255)                                                                                                                     aaGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATccc                                   >  1:620473/1‑151 (MQ=255)                                                                                                                     aaGCCTA‑AAAAAAAATATGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATc                                          <  1:837876/144‑1 (MQ=255)                                                                                                                      aGCCTG‑ACAGAAAGTATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcccc                                  >  1:1074873/1‑151 (MQ=255)                                                                                                                         cctaaaaaaaaaTATGCTGTAAGCCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt                                      <  2:509860/137‑1 (MQ=255)                                                                                                                             aaaaaaaaTATGCTGTAAGCCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCttatt                             <  2:294269/146‑1 (MQ=255)                                                                                                                               aaaaaaTATGCGGAAAGGCACATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCtt                                <  2:1223535/143‑1 (MQ=255)                                                                                                                                   aaTATGCTGTAAGGCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATcc                      <  1:16710/151‑1 (MQ=255)                                                                                                                                   aaTATGCTGTAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATcc                      <  2:1215771/151‑1 (MQ=255)                                                                                                                                   aaTATGCTGAAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATcc                      <  1:38393/151‑1 (MQ=255)                                                                                                                                      atGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTg                   >  1:171127/1‑151 (MQ=255)                                                                                                                                       tGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTg                   <  1:1048222/150‑1 (MQ=255)                                                                                                                                           gaaaacctgaaagcAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGAGTTCGCGTCCTCGATTGCGGCCCAACTGTGCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCt                    <  2:961226/132‑1 (MQ=255)                                                                                                                                             aaGTCTCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTg            <  2:895304/151‑1 (MQ=255)                                                                                                                                               ggCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGa          <  1:530514/151‑1 (MQ=255)                                                                                                                                                gCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGa          >  2:1237578/1‑150 (MQ=255)                                                                                                                                                        cAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGcc  >  1:141394/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                      AAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGTGTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTG‑ACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCAAAGTGGAACATTGCC  >  NC_000913/345179‑345477

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