Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	41.7%	89.0 / 16.9	60	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (30/5);  new base C (0/25);  total (30/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

GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACC  >  NC_000913/3556059‑3556357                                                                                                                                                          |                                                                                                                                                       gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTcc                                                                                                                                                      <  1:591421/151‑1 (MQ=255)      aaCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGcc                                                                                                                                                 >  1:879972/1‑151 (MQ=255)       accaccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGcc                                                                                                                                                 >  1:789243/1‑150 (MQ=255)         caccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAg                                                                                                                                               >  1:586249/1‑150 (MQ=255)           ccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAgtgt                                                                                                                                            >  2:102442/1‑151 (MQ=255)                         cgcTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGc                                                                                                                               >  2:585382/1‑150 (MQ=255)                             ttcccccGACAAAAAAAAAGG‑GTCGGGGGATTTTTTCCCTTTCGCCTCCAAGCGAGACGGTTTCCCCCGCCCCCTGGCCGCGCGGCAGG‑TTATA‑AATTTTCTG‑CTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                                            <  1:11091/144‑1 (MQ=255)                              ttcaccgaaacaaaaaaaGCGGTTCGTTGATTTTTCCCCTTGCAGCTCAAGGCAAGAAGGTATTCCCCGCCCCCTGTTCGGGGGGC‑GG‑TGCAA‑GATTTTCTGTTTCATGGAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                                            <  1:475357/137‑1 (MQ=255)                              ttCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                                          >  1:308766/1‑150 (MQ=255)                                 aCCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCGCCGCCAGTGTAGCGATTTCACGCTCAGCAACa                                                                                                                      <  2:451066/151‑1 (MQ=255)                                 aCCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                                       >  2:383861/1‑150 (MQ=255)                                     acgacAAAAAGGC‑GTTCGGTGATTTTTTCCCTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                                    <  2:778493/149‑1 (MQ=255)                                        acAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                               >  2:915538/1‑151 (MQ=255)                                        aaaaaaaaGC‑GTTCGGTGTTATTTTCCCTTCCGCCTTCAAGCGAGACGTTATCCCCCGGCCCCTGGTCCCGCGGCGGT‑TTATG‑AATATTCTG‑TTCGTGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                                <  1:211883/148‑1 (MQ=255)                                         cAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                              >  1:25745/1‑151 (MQ=255)                                         cAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                              >  2:666240/1‑151 (MQ=255)                                         cAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                              >  2:634721/1‑151 (MQ=255)                                          aaaaaaGC‑GTTCGGTTATTTTTTCCCTTCCCCCTCAAAGCGAGAGGGTTTCCCCCGCCCCCTGTCCCCGCGGTGGT‑TTATA‑AATATTCTG‑TCCATGAAGGGAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                                <  1:132650/148‑1 (MQ=255)                                           aaaaaGC‑TTTCGGTGATTTTTTCACTTTCGCCTCAAAGCGAGACGGTATTCCCCGCCCCCTGGCCGCGCGGCGGG‑TTATA‑AATATTCTT‑TTCATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                            <  2:557137/151‑1 (MQ=255)                                           aaaaaGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                             >  1:796385/1‑150 (MQ=255)                                               aGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                                        >  1:721356/1‑151 (MQ=255)                                                 c‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                                        >  1:149166/1‑149 (MQ=255)                                                     tCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                                    >  1:851666/1‑150 (MQ=255)                                                       ggagATTTTTTCTCTTTCGACTCCAAGCGAGCGTGTATCCCCCGACCCCGGGCCGCGCGGCAGG‑TTATA‑AATTTCCTG‑TTAGTGAAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                                  <  1:215581/147‑1 (MQ=255)                                                          gATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                                 >  2:891902/1‑148 (MQ=255)                                                            ttttttCACTTTCGACTTAAGGGGAGATTTTATTCCCCGGCCCCTGGTCGCGCGGCAGT‑TTATA‑AATATTCTT‑TTAAGGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                              <  1:305404/147‑1 (MQ=255)                                                               tttCACTTTCGACTTCAAGCGAGACGTTATTCCCCGCCCCCTGGTCGCGCGGCAGG‑TTATA‑AATTTTCTG‑TTCATGGAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                          <  2:180623/150‑1 (MQ=255)                                                                ttaaCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTCATAGAACTtc                                                                                         >  2:123620/4‑150 (MQ=255)                                                                  cACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                                      >  2:249231/1‑151 (MQ=255)                                                                    cTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                                    >  1:926132/1‑151 (MQ=255)                                                                      ttCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                                  >  1:660518/1‑151 (MQ=255)                                                                          aCTTAAAGCGAGACGGTATTCCCGGCCCCCGGGCCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCAGGCGGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                                  <  1:535549/147‑1 (MQ=255)                                                                           cTTAAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                                <  2:953700/148‑1 (MQ=255)                                                                             tCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                            >  2:191780/1‑150 (MQ=255)                                                                               aggCGAGACGTTATTCCCCGCCCCCGGGCCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGAAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                             <  2:660540/145‑1 (MQ=255)                                                                                    ggtggttatcccccGCCCCCTGTTCGCGGGGC‑GG‑TTTTGAAATTTCCTT‑TGCTTGGAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                                    <  2:531534/140‑1 (MQ=255)                                                                                     gaCGGTATTCCCCGCCCCCTGGTCGCGCGGCAGT‑TTATG‑AATTTTCTG‑TTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                                    <  2:638191/150‑1 (MQ=255)                                                                                        gttatcccccGGCCCCGGGCCGCGCGGCGGT‑TTAAA‑AATATTCTG‑TTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                                 <  1:475869/143‑1 (MQ=255)                                                                                        ggTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                                >  2:297999/1‑151 (MQ=255)                                                                                        ggTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                                >  2:649173/1‑151 (MQ=255)                                                                                             tACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                           >  1:500362/1‑151 (MQ=255)                                                                                               cccGGCCCCGGGTCGCGCGGGAGG‑TTATA‑AATTTTCTT‑TTCATGGAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                           <  1:758030/149‑1 (MQ=255)                                                                                               cccGCCCCCTGGTCGCGCGGCAGT‑TTATG‑AATTTTCTG‑TTCAGGGAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCgcgc                                                             <  1:659034/147‑1 (MQ=255)                                                                                                      ccgggTCGCGC‑GCAGGTTTTTG‑ATTTTTTTT‑TTCATGAAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                                   <  2:721382/147‑1 (MQ=255)                                                                                                        tGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGc                                                 >  1:223698/1‑150 (MQ=255)                                                                                                         ggccgcgggGCAGT‑TTATG‑AATTTTTTG‑TTCATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                               <  2:593330/148‑1 (MQ=255)                                                                                                         ggccgcgggGCAGT‑TTATG‑AATTTTCTT‑TTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                                   <  2:79401/144‑1 (MQ=255)                                                                                                                 ggCAGT‑TTATA‑AATTTTTTG‑TTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtgt                                          <  1:915509/148‑1 (MQ=255)                                                                                                                   cAGT‑TTATG‑AATATTCTG‑TTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaa                                      <  2:691932/150‑1 (MQ=255)                                                                                                                      g‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                                  >  2:336223/1‑151 (MQ=255)                                                                                                                        ttATG‑AATATTCTG‑TTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaa                                      <  2:789270/146‑1 (MQ=255)                                                                                                                               attttCTG‑TCCCGGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                                   <  2:20798/140‑1 (MQ=255)                                                                                                                               atatTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg                           >  1:354529/1‑151 (MQ=255)                                                                                                                                 atTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg                           >  1:613926/1‑149 (MQ=255)                                                                                                                                      ttttCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag                      <  1:410903/148‑1 (MQ=255)                                                                                                                                        ttCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGTTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgaga                     <  2:654286/149‑1 (MQ=255)                                                                                                                                        ttCATGAAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag                      <  1:913429/148‑1 (MQ=255)                                                                                                                                            aggcgggAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc                 <  1:155327/146‑1 (MQ=255)                                                                                                                                               aGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTg            >  2:598497/1‑151 (MQ=255)                                                                                                                                                      aa‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCcaccac     >  2:961153/1‑151 (MQ=255)                                                                                                                                                          tCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCTGTTGCCACCACAcc  >  1:204462/1‑151 (MQ=255)                                                                                                                                                          tCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACAcc  >  2:913073/1‑151 (MQ=255)                                                                                                                                                          tCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACAcc  >  1:792641/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                       GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGC‑GTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGG‑TTATG‑AATATTCTG‑TTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACC  >  NC_000913/3556059‑3556357

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