Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,376,505	0	A	C	42.9%	104.6 / 12.2	63	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (29/7);  new base C (0/27);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.12e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376355‑1376655                                                                                                                                                          |                                                                                                                                                       cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAaca                                                                                                                                                        <  2:923218/151‑1 (MQ=255)    ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                     >  1:417951/1‑151 (MQ=255)      aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCt                                                                                                                                                    <  2:175417/150‑1 (MQ=255)      aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                                   <  2:399618/151‑1 (MQ=255)        cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                                 <  2:780752/151‑1 (MQ=255)          ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                               >  2:468212/1‑151 (MQ=255)               gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGcc                                                                                                                                          >  2:53393/1‑151 (MQ=255)               gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGc                                                                                                                                           <  2:519411/150‑1 (MQ=255)                  cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                       >  2:39621/1‑151 (MQ=255)                  cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                       >  1:922458/1‑151 (MQ=255)                    cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                                      <  1:288950/150‑1 (MQ=255)                            aGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                                             >  1:977555/1‑151 (MQ=255)                             gCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                                             >  2:387976/1‑150 (MQ=255)                                 ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                         >  2:330088/1‑150 (MQ=255)                                  gACGGCGACCTCAAAATGGTGAGGGTTCCCCTGCGGGTACTGGAACC‑GGCGGCCCCCCGTCCCGTAACGGGGG‑TGTCTCCCCAACCT‑TTAACCCCCAACCGGGACAAAAGAATCGACCCCCCCGCCGATCCCCGCCGTTACATCGTCGaaa                                                                                                                       <  1:360034/151‑1 (MQ=255)                                     ggCGCCCGAAAAATGGTGAGTGATTCCCTGCGCGTACGGGAACG‑GGTCGCCGCGCGTACCGAACCCGTGG‑TGTATCTCGAACCG‑TAAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                                       <  2:730982/148‑1 (MQ=255)                                        gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                                 >  2:885364/1‑151 (MQ=255)                                        gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                                 >  2:871379/1‑151 (MQ=255)                                          ccaaaaaaGGGTGGGTGTTTCCCTGCGCGTCCTGGACCA‑GGCCGCCCCGCGCCCCGGACCCGGGG‑TTTATCGCAACCCG‑TTAACCCGTTATCAGGTAAAAATGAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                   <  1:287011/143‑1 (MQ=255)                                          ccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                               >  2:992136/1‑151 (MQ=255)                                          ccGCAAAAGGGTGAGTGATCCCCTGCGCGTACGGAAACG‑GGTCGCCGCGCGCACGGGACCGTTGG‑TGTTCCTCGACCCT‑TTAACCCGTTATCGGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                                <  1:780980/150‑1 (MQ=255)                                          ccGAAAAAGGGTGAGTGATTCCCTGCGCGTACTGAAACG‑GGTCGCCGCGCTTCCCGAAACGGTGG‑TTTATCTCAACCCT‑TTAACCCGCTATGAGGATCATATAATCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                                <  2:508793/150‑1 (MQ=255)                                              aaaaTGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                           >  1:454599/1‑151 (MQ=255)                                              aaaaTGGTAAATATTTCCCGGGGGGTGATGGAGAA‑GGTCGCAGCGAGAACCGGAACCGTGG‑GATGTCCGGACCTG‑TAAATCCCCAATGGGGAAGATATGTGCAAGACCCTCCCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                           <  2:201344/151‑1 (MQ=255)                                              aaaaTGGGGAGTGTTTCCCGGCGCGTACGGAAACG‑GGTGCCCGCGCGTCCCGGACCCGGGG‑TGTATCTCGAACCT‑TTAACCCGCTATCGGTATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                                <  2:835959/146‑1 (MQ=255)                                                   ggtgggtgTTTCCCTGCGCTTCCTGAACCA‑GGTCGCCGCGCGTCCCGGACCGGTGG‑TGTATCTCAACCCT‑TTAACCCGTTATCAGGATAATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                       <  1:613459/150‑1 (MQ=255)                                                    gggagtgtTCCCCTGCGCGTACGGAAACA‑GGTCGCCGCGCGTCCCGGACCCGTGG‑TGTATCTCGACCCG‑TAAACCCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                     <  2:648101/149‑1 (MQ=255)                                                         tgaTTCCCTGCGCGTACGGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGACCCG‑TTAACCCCTTATCAGGATCATTTGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                                  <  2:482810/149‑1 (MQ=255)                                                            ttCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                             >  1:552717/1‑151 (MQ=255)                                                             tCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                            >  2:242269/1‑151 (MQ=255)                                                              cccTGCGTGTACGGGAACA‑GGTCGCGGGGCGT‑CCAGAACCGTGGTTGTATCCCGACCCG‑TAAACCCGTTATAAGTTTAATATGATAAACCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                           <  1:679096/151‑1 (MQ=255)                                                                 tGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                         >  2:546222/1‑150 (MQ=255)                                                                  ggggggtgcggAACAGGGTCGCCGCGCGCCCGGGAACCGTGG‑TTGCCCGCGACCTG‑TAAACCCGCTATCTGGATAATATGAACAACCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                          <  1:807959/140‑1 (MQ=255)                                                                    gcgTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                     >  1:417561/1‑151 (MQ=255)                                                                    gcgTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGAATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGCTGCCCGCCGTTACATCGGCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                     >  2:697906/1‑151 (MQ=255)                                                                       tacggaaaCA‑GGCCGCCGCGCGTACCGAACCGGTGT‑TTTATCTCGAACCT‑TAAACCCGCTATCAGAATAATATAATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                    <  1:599661/143‑1 (MQ=255)                                                                       tACGGGAACA‑GGTCGCCGCGCGTCCCGGACCCGTGG‑TGTATCTCGACCCG‑TTAAACCGTTATCAGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                    <  2:806509/149‑1 (MQ=255)                                                                          tggagaa‑GGTCCCCGCGCGTCCCGAAACGGTGG‑TGTATCTCGACCGG‑TAAACCCTCTATCAGATTATTATATTCACCCCCCCCCCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                                  <  1:463708/142‑1 (MQ=255)                                                                          tGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                                >  1:166723/1‑150 (MQ=255)                                                                                  ggTCGCCGCGCGTCCCGGACCGGTGG‑TGTATCTGTAACCG‑TTAAACCGTTATAAGGATAATATGATCAACACCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                         <  1:865408/150‑1 (MQ=255)                                                                                  ggTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                        >  1:346195/1‑151 (MQ=255)                                                                                  ggTCCCCGCGCGTCCCGGAACCGTGG‑TGTATCTCGACCTT‑TTAAACCCCTAGCAGGATAATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                          <  2:218959/149‑1 (MQ=255)                                                                                    tCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                                      >  1:270517/1‑151 (MQ=255)                                                                                             gTACCGGACCCGTGG‑TGTATCTCAAACCT‑TTAACCCGCTATCAGGATAATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                              <  2:708233/150‑1 (MQ=255)                                                                                               aCCGAACCCGTGG‑TGTATCCCAACCCT‑TTAACCCCCTATAAGGATAATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                             <  1:885331/149‑1 (MQ=255)                                                                                                cgggACCCGGGG‑TGTCCCGCGACCCG‑TAAACCCGCTATCAGAATAATATGAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  1:27324/147‑1 (MQ=255)                                                                                                cgggAACCGTGG‑TGTATCTC‑AACCGTTTAACCCGTTATCAGAATAATATGATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                            <  2:98361/147‑1 (MQ=255)                                                                                                ccGGAACCGTGG‑TTTTTCTCGAACCT‑TTAACCCTCTATCAGAATCATATGAACCACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                               <  1:918808/146‑1 (MQ=255)                                                                                                 aggAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                         >  1:427031/2‑151 (MQ=255)                                                                                                  gaacccGTGG‑TTTATCTCGAACCT‑TAAAACCGCTATCAGGATAATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                           <  2:969094/144‑1 (MQ=255)                                                                                                        gtgg‑tgTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                  >  2:560033/1‑151 (MQ=255)                                                                                                             tgtATCTCGACCCG‑TAAAACCTTTATCAGGATAATATAACCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                                <  2:364239/149‑1 (MQ=255)                                                                                                                       acctgtaaacccGTTATCAGGATAATAAAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                              <  2:720575/132‑1 (MQ=255)                                                                                                                       aacG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATATTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                    >  1:756137/3‑151 (MQ=255)                                                                                                                         cG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                                   >  2:60593/1‑150 (MQ=255)                                                                                                                              aacccGCTATCAGGATAATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                              <  2:311170/148‑1 (MQ=255)                                                                                                                              aacccGCTATCAGAATAATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                              <  2:398444/148‑1 (MQ=255)                                                                                                                                   gCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                         >  1:862586/1‑151 (MQ=255)                                                                                                                                          ggATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt                  >  1:798336/1‑151 (MQ=255)                                                                                                                                          ggATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                    <  2:71010/149‑1 (MQ=255)                                                                                                                                          ggATCAAATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                         <  2:951774/144‑1 (MQ=255)                                                                                                                                               atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt             >  2:945991/1‑151 (MQ=255)                                                                                                                                                  tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc           >  2:945455/1‑150 (MQ=255)                                                                                                                                                       aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg     >  2:518217/1‑151 (MQ=255)                                                                                                                                                          acCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCgg  >  2:387175/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                       CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACA‑GGTCGCCGCGCGTACCGGAACCGTGG‑TGTATCTCGAACCG‑TTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCGG  >  NC_000913/1376355‑1376655

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