Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,001,164	0	A	C	42.4%	92.1 / 20.2	59	T274P (ACC→CCC)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base A (29/5);  new base C (0/25);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 5.49e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTC  >  NC_000913/3001015‑3001313                                                                                                                                                       |                                                                                                                                                       cATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAAc                                                                                                                                                       <  2:274312/151‑1 (MQ=255) cATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAAc                                                                                                                                                       >  2:751227/1‑151 (MQ=255)     cGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGt                                                                                                                                                   <  2:986256/151‑1 (MQ=255)            gTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGcc                                                                                                                                           >  1:46946/1‑151 (MQ=255)                tcaAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGc                                                                                                                                       >  2:817142/1‑151 (MQ=255)                 caAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCt                                                                                                                                      >  2:496357/1‑151 (MQ=255)                      cATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTAc                                                                                                                                 >  1:406269/1‑151 (MQ=255)                        tttgggggTGGGGGTTTTGGTAATAAACAGGTTGTACGGAAGGGGCCAAGGGCGGCATCCCTGACCAGAAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCGGTAAAGGGTGTTTCCCCGCACCCCGTC‑CCCGCCACGCTTTTAc                                                                                                                                 <  1:326759/143‑1 (MQ=255)                          tGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCAtt                                                                                                                             >  1:325091/1‑151 (MQ=255)                              gggggCGTTTTTGTAAAAAAACGGGATTTACGGGAAGAGCCAAGGGCGCTATCCCTACCCAGAAGGCTGGGCGGCATCCCGGTGAAAGTTTCCC‑TTAGCTGTAAAGTGTTTTCCCCCGCACCCCGTA‑CCCGCCACGCTTTTACCATTGACg                                                                                                                         <  1:808567/148‑1 (MQ=255)                                 ggCGGTTTTGATAAAAAACAGGATTTACGGGAGGAGCAAATGGCGGCATCCCGGACCACCAAGTTGGGCGCCTTCCCGGTAAAATTTTCCC‑TTAGCCGTAAAGAGTTTTCCCCCGCACCCCGTA‑CCCGCCACGCTTTTACCATTGACgg                                                                                                                        <  2:174415/149‑1 (MQ=255)                                     gTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCaa                                                                                                                    >  1:679791/1‑149 (MQ=255)                                     gTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAAt                                                                                                                  >  1:304323/1‑151 (MQ=255)                                     gTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAAt                                                                                                                  >  1:879226/1‑151 (MQ=255)                                     gTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAAt                                                                                                                  >  1:882815/1‑151 (MQ=255)                                      ttttGGAAATAAACAGGATGTAGTGGAAGAGCAAATGGGGGCTTCCCTGACCAGCAAGCTTGGGGGCATTCCGGAAAAAGTTCCCT‑TTAGCCGTAAAGATTTTTCCCTCGCACCCCTTA‑CCCGCCACGCTTTTACCATTGACGGGCaaa                                                                                                                   <  2:761982/149‑1 (MQ=255)                                               aacacggtaTGTACGGGAAGGGCAAGTGGCGGCTTCCCGCCCCAGCAAGCTTGGCGCCATTCGGGTGAATGTTCCCC‑TGCGCGGTGAAGAGTTTTTCCCCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAc                                                                                                        <  1:588785/143‑1 (MQ=255)                                                aaaCAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAcc                                                                                                       >  1:889131/1‑151 (MQ=255)                                                 aaCAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcg                                                                                                      >  2:670104/1‑151 (MQ=255)                                                  aCAGGATGTACTGGAAGAGCCAATGGGGCCATTCCTGCCAAGCAAGCTTGGCGGCATTCCGGTGAAATTTTCCC‑TTAGCCGTAAAGGGTGTTTCCCCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                                                                                     <  2:167446/151‑1 (MQ=255)                                                  aCAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                                                                                     >  1:932908/1‑151 (MQ=255)                                                      gATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACg                                                                                                 >  2:201733/1‑151 (MQ=255)                                                       aTGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACgg                                                                                                >  2:655134/1‑151 (MQ=255)                                                           aCTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAAc                                                                                             >  1:666891/1‑150 (MQ=255)                                                             tGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                          >  1:293354/1‑151 (MQ=255)                                                              ggAAGGGCAAAGGGCGCCATCCCTGCCCAGCAAGTTGGGGGGCATCCGGGTAAAATTTTCCC‑TTAGCGGTAAAGAGTTTTCCCCCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACa                                                                                            <  2:156054/148‑1 (MQ=255)                                                              ggAAGACCAAAGGGCGGCTTCCCTCACCAGCAAGTTGGGCGGCTTTCCGGAAAAATTTCCCT‑TTCCCCGTGAAGGGTTTTCCCCCGCACCCCTTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTg                                                                                         <  2:620718/151‑1 (MQ=255)                                                              ggAAGACCAAAGGGCGGCTTCCCGGACCAGCATTCTGGGCGGCTTTCGGGAGAAATTTCCCC‑TTAGCGGTAAAGGGTGTTCCCCCGCAACCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAAc                                                                                             <  1:825542/147‑1 (MQ=255)                                                                aaGAGCAAATGGGGGCTTCCCTCCCCAGAAGGCTTGGCGGCTTCCCGGTAAAAGTTCCCC‑TTACCGGTAAAGGTTGTTCCCCCGCACCCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                          <  2:12400/148‑1 (MQ=255)                                                                    gaaaagggcggcTTTCCGGACCAGCAAGCTGGGGGGC‑TTCCGGGAGAATTTTCCCTTTAGCGGGGAAGGGTGTTCCCTCGCACCCCGCA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa                                                                                      <  1:417800/142‑1 (MQ=255)                                                                       aagggcggcTTTCCTGCCCAGCAAGCTTGGCGCCTTTCCGGTAAAATTTCCCC‑TTAGCCGTGAAGAGTTTTCCCCCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAgg                                                                                    <  2:285846/144‑1 (MQ=255)                                                                       aaTGGCGGCCTTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                                >  1:950842/1‑151 (MQ=255)                                                                          ggcggcATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGt                                                                             >  1:416549/1‑151 (MQ=255)                                                                                    cggACCAGCATTGTTGGGGCCATCCGGGTAAAAGTTTCCC‑GTCTCGGGGAAGAGTTTTCCCGCCCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGa                                                                        <  1:1004459/144‑1 (MQ=255)                                                                                       aCCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGt                                                                >  2:50410/1‑151 (MQ=255)                                                                                        ccAGCAAGCTGGGCGGCATCCCGGTGAATGTTTCCC‑TTAGCGGGAAAGAGTTTTCCCTCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTg                                                                 <  2:864531/149‑1 (MQ=255)                                                                                          aGCAATCTGGGCGGCTTTCGGGTGAAAGTGTCCC‑TTGGCCGTAAAGGGTTTTCCGCCGAACCCCGCA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTa                                                             <  1:754247/151‑1 (MQ=255)                                                                                           gCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTa                                                             >  1:148469/1‑150 (MQ=255)                                                                                             aagttgggcggcTTTCGGGTGAAAGTTCCCC‑TTAGCGGTGAAGAGTGTTCCCTCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAaca                                                          <  2:1007476/145‑1 (MQ=255)                                                                                             aaGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAac                                                           >  2:560573/1‑150 (MQ=255)                                                                                                cTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCg                                                       >  2:326296/1‑151 (MQ=255)                                                                                                   gggggcATTCGGGTGAAAGTTCCCC‑TGGGCGGGGAGGGGTTTTTCCTCGCACCCGGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAaca                                                          <  1:485188/142‑1 (MQ=255)                                                                                                   ggcggcATCCGGGAAAAAGTTTCCC‑TTACCCGTAAAGAGTGTTCCCTCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCg                                                       <  2:890673/148‑1 (MQ=255)                                                                                                           ccgggaaaaaTTTTCCC‑TTAGCCGTAAAGTGTTTTCCCCCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCa                                            <  1:857496/144‑1 (MQ=255)                                                                                                            cgggTGAAAGTTCCCC‑TTAGCCGAGAATTTTTTTCCCTCCCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGc                                             <  1:925925/147‑1 (MQ=255)                                                                                                            ccGGTAAAATTTCCCC‑TTAGCCGTGAAGGGTTTTCCCCCCCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAt                                           <  1:824759/151‑1 (MQ=255)                                                                                                               gTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                        >  1:701746/1‑151 (MQ=255)                                                                                                                  aaaGTTCCCC‑TTCGCCGTGAAGGGTTTTTCCCCCCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACg                                     <  1:163747/151‑1 (MQ=255)                                                                                                                  aaaGTTCCCC‑TTAGCCGTGAAGGGTGTTCCCCCCCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCa                                       <  2:1018025/149‑1 (MQ=255)                                                                                                                   atgTTCCCC‑TTAGCTGTAAAAAGTGTTCCCGCGCACCCCTCC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCa                                       <  2:45276/146‑1 (MQ=255)                                                                                                                   aattttCCCTTTGGCCGGGAAGAGTTTTTCCCCGC‑ACCCGTACCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACg                                     <  1:836858/148‑1 (MQ=255)                                                                                                                   aaGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACgg                                    >  2:58143/1‑151 (MQ=255)                                                                                                                   aaGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACgg                                    >  2:979597/1‑151 (MQ=255)                                                                                                                                   gTGAAGAGTTTTCCCCCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATc                          <  2:744906/146‑1 (MQ=255)                                                                                                                                     gAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTg                   >  2:51853/1‑151 (MQ=255)                                                                                                                                            gtTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggggg            >  2:559998/1‑151 (MQ=255)                                                                                                                                                ccccGCACCCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggggg           <  2:416592/145‑1 (MQ=255)                                                                                                                                                 cccGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATa       <  2:177656/149‑1 (MQ=255)                                                                                                                                                     cAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAg    >  2:490957/1‑150 (MQ=255)                                                                                                                                                      aaCCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTc  >  2:491657/1‑151 (MQ=255)                                                                                                                                                      aaCCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTc  >  2:905646/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTC  >  NC_000913/3001015‑3001313

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