Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,640,981	0	C	T	43.4%	51.5 / 36.3	54	F40F (TTC→TTT)	dtpB	dipeptide and tripeptide permease B
Reads supporting (aligned to +/- strand):  ref base C (3/27);  new base T (23/0);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 4.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGT  >  NC_000913/3640832‑3641130                                                                                                                                                        |                                                                                                                                                           tttCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTttc                                                                                                                                                            <  1:979360/150‑1 (MQ=255)  ttCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTttct                                                                                                                                                           <  2:963930/150‑1 (MQ=255)        aTGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTaa                                                                                                                                                  <  2:455792/151‑1 (MQ=255)             aTTTGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTT‑G‑TTaa                                                                                                                                                  >  2:123155/1‑146 (MQ=255)             aTTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTCTTT‑C‑TTTAACaa                                                                                                                                              >  2:938614/1‑149 (MQ=255)                ggAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTGTTCTG‑TCAAACattt                                                                                                                                            >  1:593043/1‑147 (MQ=255)                ggAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑Ctt                                                                                                                                          <  2:575775/150‑1 (MQ=255)                  aaGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTC‑T‑TTAAACAA‑CTTgg                                                                                                                                        >  2:352012/1‑150 (MQ=255)                     gTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTC‑C‑ttta                                                                                                                                                  >  2:356191/1‑136 (MQ=255)                      ttAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTT‑G‑TTTAAAAG‑CTTGGGtttt                                                                                                                                   >  2:534260/1‑149 (MQ=255)                       taataaATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTT‑G‑TTtaaaaacctggggttt                                                                                                                                    >  2:456517/1‑137 (MQ=255)                        aataaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGG‑TTTTTTTC‑T‑TTaaaagcgtgggttttc                                                                                                                                    >  2:49252/1‑137 (MQ=255)                         ataaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGc                                                                                                                                <  1:640983/151‑1 (MQ=255)                               atGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTC‑G‑TTTAAAAG‑CTTGGGTTCTCGCAaga‑g                                                                                                                          >  2:249684/1‑150 (MQ=255)                                     acaacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTT‑G‑TTTAAAAG‑CTTTGGGTCTCCCCAAA‑GCAGGCtt                                                                                                                   >  2:1040258/1‑151 (MQ=255)                                      caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCttt                                                                                                                  <  1:630400/151‑1 (MQ=255)                                      caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCttt                                                                                                                  <  1:994127/151‑1 (MQ=255)                                        acaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTT‑T‑TTTAAAAG‑CCTGGGTTTTCGCTACA‑GCAgggctt                                                                                                                  >  1:589850/1‑145 (MQ=255)                                         caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGt                                                                                                               <  2:1204777/151‑1 (MQ=255)                                          aacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTc                                                                                                              <  1:291300/151‑1 (MQ=255)                                              cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACtt                                                                                                          <  2:452567/151‑1 (MQ=255)                                               ccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTT‑G‑TTTAAAAG‑CTTTGGTTTTTGCCAAA‑GCAGGGCTTTTTTACttt                                                                                                         >  2:132552/1‑151 (MQ=255)                                                   ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTgg                                                                                                      <  1:702314/150‑1 (MQ=255)                                                      gATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCt                                                                                                  <  1:570485/151‑1 (MQ=255)                                                       aTGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCtt                                                                                                 <  2:770500/151‑1 (MQ=255)                                                       aTGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑CTTTTTTC‑GTTTAAAAAC‑CTTGGGTT‑TCGCCAGA‑GCAGGCCTTTTTTTCtttttggg                                                                                                    >  2:613103/1‑144 (MQ=255)                                                         gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTAACTTTTGGTGCttt                                                                                                <  2:757645/150‑1 (MQ=255)                                                              agcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTT‑G‑TTTAAAAA‑CGTGGAGTCTCTCAAGA‑GCAGGGTGTTTTTACTTTTTGGGCtttt                                                                                               >  2:351004/1‑146 (MQ=255)                                                              agcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTtgctgc                                                                                          <  2:729612/151‑1 (MQ=255)                                                                  aCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTg                                                                                      <  1:697971/151‑1 (MQ=255)                                                                  aCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTg                                                                                      <  1:452813/151‑1 (MQ=255)                                                                      cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCt                                                                                  <  1:127085/151‑1 (MQ=255)                                                                        cccATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTC‑T‑TTAAAAAC‑CTGGGGTT‑TCGCCAAA‑GCGGGCGTTTTTTTATTTTTGTGCCTTTGTTTCGCTgggcgtt                                                                                >  1:512240/1‑146 (MQ=255)                                                                         ccATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑Tg                                                                              <  2:937115/151‑1 (MQ=255)                                                                               ttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTT‑G‑TTTAAAAA‑CTTTGGTTTTCCCCAGAGGAAGGCTTTT‑TTACTTTTTGTGGTTTTTCTTCGCCGGGCttttggc                                                                            >  2:905591/1‑142 (MQ=255)                                                                                tCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTT‑G‑TTTAAAAA‑CTTTGATTTTTGCAAAA‑GCAGGCCTTTTTTACTTTTTGTGCTTTTTCTTTGCTGGGCTA‑TGGcc                                                                           >  1:497879/1‑147 (MQ=255)                                                                                tCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCa                                                                        <  1:148778/150‑1 (MQ=255)                                                                                 cATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑CTTTTTTT‑G‑TTTAACAG‑GTTGGGTTTTTGCTAGA‑GCAGGGTTTTGTCCCTTTTTGGGCGTTTTCTGCGCCGGGCTA‑TGGCCTCCtt                                                                      >  2:243554/1‑151 (MQ=255)                                                                                   tGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTT‑G‑TTTAAAAGATTTGGTTTTTGGAAAAG‑GAAGGCTTTTGTCCTTTTTGGGGTTTTTTTTGCGCGGGCTTA‑TGG‑CTCAttt                                                                     >  2:696839/1‑150 (MQ=255)                                                                                        ttttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTT‑G‑TTTAAAAG‑CTTTGTTTTTCGCGAGA‑GCAGGGGTTTTTCCTTTTTTTTGGTTTTTGTGCGCCGGGCTT‑TTGGCTCCttt                                                                     >  1:425719/1‑145 (MQ=255)                                                                                          ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTCTTT‑G‑TTTAAAAA‑CTTTGGTTCTCCCAAAA‑GCAGGGCTTTTTCACTTTTTGTGCTTTTTTTTCGCTGGGCTATTGGCCCCTTTTcctt                                                                 >  1:1110252/1‑146 (MQ=255)                                                                                            ttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTCTTT‑G‑TTTAAAAA‑CTTGGGTTTTCGCCAGA‑GCAGGCCTTTTTCACTTTTTGTGCTTTTTCTGTGCTCGGCTT‑TGGCCTCCTTTCCTTTtgcgg                                                           >  2:61214/1‑151 (MQ=255)                                                                                             tGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTT‑G‑TTAAAAAG‑CTTTGGTTCTCGCAAGA‑GCAGGGTTTTTTCACTTTTTGTGCTTTTTCTTCGCTGGTCTA‑TTGCCTCATTTTCATTggcgg                                                           >  2:854017/1‑150 (MQ=255)                                                                                               tCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCt                                                         <  1:561531/150‑1 (MQ=255)                                                                                                cGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTT‑G‑TTTAACAA‑CTTGGGTTTTCGCAAAA‑GCAAGCTTTTTTTACTTTTTGTGGTTTTTTTTCGCTGGGCTA‑TGGGCCCAttttcctttgggg                                                           >  2:1017369/1‑137 (MQ=255)                                                                                                cGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTT‑G‑TTTAAAAG‑CCTGGGTTTTCGCCAGA‑ACAGGCTTTTTTCCCTTTTTGTGCTTTTTCTGCGCTGGGCTT‑TGGCCTCCttttcctttg                                                              >  2:429455/1‑137 (MQ=255)                                                                                                cGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTAt                                                       <  2:725023/151‑1 (MQ=255)                                                                                                 gAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTC‑T‑TTTAACAG‑CTTGGGTTTTCGCCAGA‑GCAGGCTTTTTTTACCTTTTGTGGTTTTGCTGGGGTGGTGTATTGGCCTCtttttttttgcgg                                                            >  2:492145/1‑137 (MQ=255)                                                                                                     tATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTT‑C‑TTAAAAAA‑CTTGGGTTTTCGCCAGA‑GCAGGGGTTTTTCCCTTTTTGTGCTTTTTCTGCGCTGGTGTT‑TTGCCTCATTTTCAtttgggggctt                                                       >  2:214309/1‑137 (MQ=255)                                                                                                      aTGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCgg                                                  <  2:1226667/150‑1 (MQ=255)                                                                                                          gAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGacc                                             <  2:172876/151‑1 (MQ=255)                                                                                                                    gCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTggg                                   <  2:262639/151‑1 (MQ=255)                                                                                                                          acGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGAcc                               <  1:143856/149‑1 (MQ=255)                                                                                                                           cGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCaa                             <  1:104003/150‑1 (MQ=255)                                                                                                                                    gggCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCAtt                   <  2:845583/151‑1 (MQ=255)                                                                                                                                             ggC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTgg           <  2:379325/150‑1 (MQ=255)                                                                                                                                                      ttcttc‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACtt    >  1:807882/1‑150 (MQ=255)                                                                                                                                                       tcttc‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTtg   <  1:475868/150‑1 (MQ=255)                                                                                                                                                        cttc‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTtgt  >  1:996073/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                           TTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTC‑G‑TTAAACAG‑CTTGGATTCTCGCAAGA‑GCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTA‑TGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGT  >  NC_000913/3640832‑3641130

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