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

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAT  >  NC_000913/1376355‑1376650                                                                                                                                                              |                                                                                                                                                  cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAaca                                                                                                                                                   <  1:646013/151‑1 (MQ=255)   cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACAcc                                                                                                                                                 >  2:551339/1‑151 (MQ=255)     tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACAcc                                                                                                                                                 <  2:922665/149‑1 (MQ=255)        cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGc                                                                                                                                            >  2:605416/1‑151 (MQ=255)        cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGc                                                                                                                                            >  1:739103/1‑151 (MQ=255)         gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGcc                                                                                                                                           <  2:151199/151‑1 (MQ=255)             aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATg                                                                                                                                       >  2:767387/1‑151 (MQ=255)                         cgTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                            >  1:498925/1‑150 (MQ=255)                            agccgggcgggcacccGCAAAAT‑GGGGGGTGATCCCCGGCGCTTACTGGAACGGG‑TCGCCGCGCGTCCGGGAA‑CCGTGG‑TGTATCCCAAACCTTAAAA‑CGCC‑TATCAGGATC‑ATATGCCCACCCCCCCCGCCGATGCCCGCCGTTACAtc                                                                                                                         <  1:408300/136‑1 (MQ=255)                             gCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                                        >  1:103084/1‑150 (MQ=255)                                 ggcgggcgcccaaaaaaT‑GGTGATTGTTCCCCTGCGCTTCCTGGAACGGG‑CCGCCGCGCGTACCGGAA‑CCGGGT‑TTTTCCTCGACCCTTAAAC‑CCGC‑TATCAGGATC‑ATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                    <  1:1019054/137‑1 (MQ=255)                                 ggcgggCGACCGAAAAGT‑GGTGAGTGTTTCCCTGCGCGTACTGGAACAGT‑TCCCCGCGCGT‑CCGGGA‑CCCTGGTTTTATC‑CCAACCCTTAAACCCGT‑TTTAAGAATA‑ATATGAAAACCCCCCCCGCCGTTCCCCGCCGTTACATCGTCGaa                                                                                                                   <  1:997170/147‑1 (MQ=255)                                 ggACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                    >  2:467904/1‑150 (MQ=255)                                    cGGCGCCGGCAAATT‑GGTGGTTGTTCCCCGGCGGGTACTGGAACGGG‑TCGCCGCGCGTACCGGAA‑CCGGGG‑TTTTCCTCGACCCGTAAAC‑CCGC‑TATCAGGATA‑ATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                                  <  1:383481/149‑1 (MQ=255)                                     ggggcccaaaaaaG‑GGTGGGTGATCCCCTGCGCGTACGGGAACGGG‑TCGCCGCGCGTCCGGGAC‑CCGGGG‑TTTATCCCGAACCGTAAAC‑CCGC‑TTTCAGAATC‑ATATAACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                <  1:294559/141‑1 (MQ=255)                                     ggCGACCGAAAAGT‑GGGGAGTAATCCCCGGCGCGTACGGGACCGGG‑TCCCCGCGCGTCCCGGAA‑CGGTGG‑TTTATCCCGAACC‑TTAAA‑CCCCTTATCAGGATA‑ATATGATCAACCCCCCCGCCGAGGCCCGCCGTTACATCGTCGaaaa                                                                                                                 <  1:123275/149‑1 (MQ=255)                                     ggCGACCGAAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CGGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                               <  2:804594/151‑1 (MQ=255)                                       cgcccGAAAAAG‑GGTGATTGATTCCCTGCGCTTACGGGAACGGT‑TCGCCGCGCGTCCCGGAC‑CGGTGG‑TGTATCTCGAACCGTTAAC‑CCGC‑TATCAGGATA‑ATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                              <  2:388898/147‑1 (MQ=255)                                        gACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                            >  1:743817/1‑151 (MQ=255)                                         aCCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                           >  1:34327/1‑151 (MQ=255)                                           ccgaaaaTGGGTGAGTGATCCCCGGGGGGTAGGGAACCGGG‑CCGCCGCGCGT‑CCGGAACCCGGGG‑TGTCTCTCGAACCGTAAAC‑CCGC‑TATCAGGATC‑ATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                             <  1:582760/145‑1 (MQ=255)                                              aaaaT‑GGTGCGGGATTCCCGGCGCGTACAGAATCAGG‑TCCCCGCGTGTACCGGAC‑CCGTGG‑TGTACCGCGAA‑CGTTAAACCCCC‑TAACAGGCTA‑ATATGACCAACCCCCCCGCCTATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                          <  1:403043/147‑1 (MQ=255)                                              aaaaG‑GGTCATTGTTTCCTTGGGCGTCCAGGACGAGG‑CCGCCGCGCGTCCAGGAC‑CGGTGG‑TGTATCTCG‑ACCGTTAAA‑CCCCTTATCGGGATC‑AAAAGACAAACCCCCGCCCGGTTGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                             <  1:383966/144‑1 (MQ=255)                                                    gtggattctttCCCTGCGCGTACGGGAACAGC‑GCCCCGCGGGTCCCGGAA‑CGGTGG‑TGTTTCCCGAACCTTTAAA‑CCGC‑TATCAGATTC‑ATAGGAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                  <  2:557478/141‑1 (MQ=255)                                                    ggTGAGTCCTTCCTGGCGTTTATTGAAACAGG‑CCGCCGCGTTTCCTAAAC‑CCGGGT‑TGTATCCAGATCTTTTAAC‑CCGT‑TATAGGGTTC‑ATATAAACAACCCCCTCGCAGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                    <  2:697871/148‑1 (MQ=255)                                                     gtgagtgtTTCCCGGCGCGCA‑GGGAACCGGCCCGCCGGGCGTCCGGGAC‑CCGGGG‑TTTCTCGCAACCCTTAAAA‑CCCC‑TAT‑AGGAACAATATACTCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                                 <  1:794145/150‑1 (MQ=255)                                                              tCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                       >  1:133893/1‑151 (MQ=255)                                                               cccTGCGCGTACTGGAACGGG‑TCGCCGCGCGTCCCGGAC‑CCGTGG‑TGTATCTCGACCCTTTAAC‑CCGC‑TATCAGGATC‑ATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                       <  2:987930/150‑1 (MQ=255)                                                                ccTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                     >  2:981804/1‑151 (MQ=255)                                                                       gttcTGGAACAGG‑TCGCCGCGCGTCCGGAAA‑CCGGGT‑TTTATCTCGACCCTTAAAC‑CCTC‑TATAAGGATA‑ATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                                   <  2:258980/143‑1 (MQ=255)                                                                       gTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                              >  2:221791/1‑151 (MQ=255)                                                                            ggAACGGG‑TCGCCGCGCGTACCGGAA‑CGGTTG‑TGTATCTCGAACCTTTAAA‑CCGC‑TATCAGGATA‑ATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                          <  2:373331/150‑1 (MQ=255)                                                                            ggAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                         >  1:705366/1‑151 (MQ=255)                                                                              aaCAGG‑CCGCCGCGCGTCCCGGAC‑CCGTGG‑TGTATCTC‑AACCCTTAAACCCGC‑TATCAGGATC‑ATATAACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                             <  2:655030/145‑1 (MQ=255)                                                                                 aGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                     >  2:691554/1‑150 (MQ=255)                                                                                 aGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                     >  2:943510/1‑150 (MQ=255)                                                                                 aGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                     >  1:807189/1‑150 (MQ=255)                                                                                  agttCGCCGCGCGTACGGGAA‑CCGTGG‑TGTATCCCAAACCGTAAAC‑CCGC‑TATCAGGTTC‑ATAGGCTCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                      <  2:207477/146‑1 (MQ=255)                                                                                   g‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                    >  1:250433/1‑149 (MQ=255)                                                                                            gggTACGGGAC‑CGGTGG‑TGTTTCCCGACCCTTTAAC‑CCCT‑TTAAAGGAAC‑ATATTATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                           <  1:1016199/148‑1 (MQ=255)                                                                                             ccggcctaaaccggtgg‑ttTACCACGAACCGTT‑AA‑CCGCTTTTCAGGATC‑ATAAACCCACCTCCCCCGATGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                           <  1:623293/137‑1 (MQ=255)                                                                                              gtcccGGAC‑CCGGGG‑TTTATCTCGACCCGTTAAC‑CCGC‑TTTCAGGATC‑ATATGATCAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                           <  2:941165/145‑1 (MQ=255)                                                                                                acgggAC‑CCGTGG‑TTTATCTCGCAGTGTAAAC‑CCGC‑TATAAGGATA‑AGATTATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                      <  2:454297/148‑1 (MQ=255)                                                                                                 ccGGAC‑CCGTGG‑TGTATCCCGACCCTTTAAC‑CCGT‑TACCGGGATA‑TTATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                        <  2:300269/132‑1 (MQ=255)                                                                                                      a‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                >  1:309768/1‑151 (MQ=255)                                                                                                           tgg‑tgTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                             >  1:777825/1‑150 (MQ=255)                                                                                                                        aaCCGTTAAC‑CCCC‑TATCAGAATA‑ATATAACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                 <  2:267129/150‑1 (MQ=255)                                                                                                                         aCCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                               >  1:1005569/1‑151 (MQ=255)                                                                                                                              tAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                          >  1:466399/1‑151 (MQ=255)                                                                                                                               aaa‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCTGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                          <  2:840254/150‑1 (MQ=255)                                                                                                                                   ccGC‑TACCAGTATA‑ATATGATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGa                       <  2:259394/150‑1 (MQ=255)                                                                                                                                     gC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                    >  2:175756/1‑151 (MQ=255)                                                                                                                                      c‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                   >  1:153988/1‑151 (MQ=255)                                                                                                                                        tATCAGGATC‑ATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                    <  1:132897/149‑1 (MQ=255)                                                                                                                                            agtaTC‑ATATAACCACCCCCCCCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                   <  1:792728/143‑1 (MQ=255)                                                                                                                                            aGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg              >  2:219447/1‑151 (MQ=255)                                                                                                                                                tC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt          >  2:463825/1‑151 (MQ=255)                                                                                                                                                tC‑ATATGACAACCCCCCGCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg              <  1:707921/147‑1 (MQ=255)                                                                                                                                                     atGATCAACCCCCTCGCCGCTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg       <  1:707415/150‑1 (MQ=255)                                                                                                                                                     atGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc      >  2:783160/1‑151 (MQ=255)                                                                                                                                                        aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat    >  1:739196/1‑150 (MQ=255)                                                                                                                                                         tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat  >  1:287971/1‑151 (MQ=255)                                                                                                                                                              |                                                                                                                                                  CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAAT‑GGTGAGTGATTCCCTGCGCGTACTGGAACAGG‑TCGCCGCGCGTACCGGAA‑CCGTGG‑TGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATC‑ATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAT  >  NC_000913/1376355‑1376650

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