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.5%	86.9 / 29.0	68	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (32/3);  new base C (0/33);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 2.60e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376357‑1376648                                                                                                                                                         |                                                                                                                                                cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACAcc                                                                                                                                               >  1:815948/1‑151 (MQ=255)   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACAccc                                                                                                                                              <  2:1197410/150‑1 (MQ=255)          gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGAt                                                                                                                                      >  1:271361/1‑151 (MQ=255)              aCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGcca                                                                                                                                  <  2:871653/151‑2 (MQ=255)                    ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                            >  2:1040351/1‑151 (MQ=255)                             cTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                   >  1:572414/1‑151 (MQ=255)                               ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                  >  2:790436/1‑150 (MQ=255)                               ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                  >  2:421845/1‑150 (MQ=255)                                   ggcgcccaaaaaaTGGGGATTGTTTCCCTGCGCGTATGGAAACGGTTCGCCGCGCGTCCCGAAACG‑GGGGTTTATC‑CCAACCCTTAAACCCGT‑AAAAAGGTT‑AATATAACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                              <  2:1281652/141‑1 (MQ=255)                                   ggCGCCCGAAAAAGGTGGAGTGATCCCCTGCGCGTGCGGAAACGGGTCGCCGCGCGTACGGGACCG‑GGGGTGTTTCCCGACCCTTTAAC‑CCGC‑TATCAGGAT‑AATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGaaaa                                                                                                               <  2:401568/149‑1 (MQ=255)                                     cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                           >  2:1055741/1‑151 (MQ=255)                                       aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                         >  2:847646/1‑151 (MQ=255)                                        ccGAAAAATGGTGAGTGTTCCCCTGCGCTTCCTGGAACGGGTCCCCGCGCGCCCCGGACCC‑GTGTTTTATCTCGACCCTTAAAC‑CCGC‑TATCAGTAT‑AATATAAACACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                          <  1:1298453/149‑1 (MQ=255)                                           aaaaaGGGGGAGTGTTCCCCGGCGCGTCCTGGAACGGGCCGCCGCGCGTCCGGGACCC‑GGGGTGTATCCCGACCCGTAAAC‑CCGC‑‑AACCGGATCCATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                       <  1:551105/148‑1 (MQ=255)                                            aaatgggggagtgaTCCCCGGCGCTTGCGGAAAGAGGTCGCCGGGCGTACCGGAACG‑GTGTTTTATCCCGACCCGTAAAC‑CCGC‑TATGAGATT‑CAAAGGAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                      <  1:1099262/141‑1 (MQ=255)                                             aaagggtagttgtttCCCTGCGTGTCCGGGACCGGGCCGCCCCGCGCCCCGGACCG‑GTGGTTTTTCCCGACCCTTAAAC‑CCGT‑TATCAGATT‑AATATGATAAACCCCCGCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                   <  1:842017/138‑1 (MQ=255)                                             aaaTGGTGGGTGTTTCCCGGCGCGTTCTGGAGCGGGCCGCCGCGCTTCCCGGAACG‑GTGGTTTATCTCGACCCTTAAAA‑CCGC‑TATCAGGAT‑CATAGGATAAACACCCTCCCGGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                   <  2:41702/151‑1 (MQ=255)                                              aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                  >  2:839451/1‑151 (MQ=255)                                                 ggggagtgtTTCCCTGCGCGTACGGGACCGGGCCGCCGCGCGTCCCGAACCG‑GTGGTGTATCTCGACCCTTAAAC‑CCGT‑TATAAGGAT‑AATATAACCAACCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                 <  1:523170/146‑1 (MQ=255)                                                      gggaTCCCCGGCGCGTACTGAAAAAGGTCGCCCCGCGTACCGAACCG‑GGGGTGTACCTCAAGCCGTAAAC‑CCGC‑TATCAGGAT‑AAAATGATCAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                          <  1:430416/149‑1 (MQ=255)                                                           tCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                     >  1:116657/1‑151 (MQ=255)                                                                  gcgTACTGGAACGGTCCGCCGGGCGTCCCGGACCC‑GTGGTTTATCTCGACCCGTAAAC‑CCGC‑TATAAGGAT‑CATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                <  2:882686/149‑1 (MQ=255)                                                                  gcgTACTGGAACGGGTCGCCGCGCGTCCCGGAACC‑GTGGTGTATCTCAACCCGTAAAC‑CCTC‑TACCAGGAT‑CATATGAAAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                               <  2:52420/150‑1 (MQ=255)                                                                   cgTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                             >  2:986866/1‑151 (MQ=255)                                                                   cgTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                             >  1:299442/1‑151 (MQ=255)                                                                     tACGGGAACAGGTCGCCGCGCGTACCGGAACG‑GTGGTGTATCTCGACCCGTTAAC‑CCGC‑TATCAGGAT‑CATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                             <  2:274837/149‑1 (MQ=255)                                                                      aCTGGACCAGGTCGCCGCGCGTCCCGGACCC‑GTGGTTTATCTCGACCCTTAAAA‑CCGC‑TACCAGGAT‑CATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                           <  2:918094/150‑1 (MQ=255)                                                                      aCTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                          >  1:412889/1‑151 (MQ=255)                                                                        ggaaccgggccgccgcGCGTCCCGGACCG‑GTGTTTTTCCTCG‑ACCCTTAAA‑CCCCTTATAAGAAT‑AATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                        <  1:268128/141‑1 (MQ=255)                                                                         ggAACGGGCCGCCGCGCGTACCGGAACC‑GTGGTGTTCCTCGACCCTTTAAA‑CCGC‑TATCAGGAT‑AATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                          <  2:1280134/148‑1 (MQ=255)                                                                          aaaCAGGCCGCCGCGAGTA‑CGGAACCGGTGGTGTA‑CCCAACCCTATAAACCCGT‑TATCAGAAT‑CTTATGATCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                          <  1:669693/146‑1 (MQ=255)                                                                              aGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                   >  2:1026322/1‑150 (MQ=255)                                                                               ggTCGCCGCGCGTCCCGGACCG‑GTGGTGTATCTCGACCCGTTAAC‑CCGC‑TATCAGGAT‑CATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                  <  2:1199501/150‑1 (MQ=255)                                                                               ggTCGCCGCGCGTCCCGAACCC‑GTGGTGTATC‑CCAACCCTTAAACCCGC‑TATCAGGAT‑AATATAACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                    <  2:61145/148‑1 (MQ=255)                                                                                  cgcggcGCGTCCGGGACCG‑GTGGTTTATCTCAACCCTTTAAC‑CCGT‑TATCAGGAT‑AATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                              <  2:868371/151‑1 (MQ=255)                                                                                     cgcgcgTCCCGGAACC‑GTGGTGTATCCAGACCCGTTAAC‑CCGC‑TATCAGGAT‑AATATGAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                           <  2:683462/151‑1 (MQ=255)                                                                                      gcgcgTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                            >  2:946536/1‑149 (MQ=255)                                                                                       cgcgTCCCGGACCG‑GTGGTTTATCTAGACCCTTTAAC‑CCTC‑TATCAGTAT‑CATATGACAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         <  2:1255646/151‑1 (MQ=255)                                                                                        gcgTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         >  1:334631/1‑150 (MQ=255)                                                                                        gcgTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         >  1:222996/1‑150 (MQ=255)                                                                                        gcgTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         >  1:193183/1‑150 (MQ=255)                                                                                             cctggacctttggtgTCCCTCGCCCCTTAAAA‑CCAT‑TATAAGAAA‑CAAATGACAAACCCCCTCCCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                       <  2:532078/138‑1 (MQ=255)                                                                                                 aaCC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                               >  1:378922/1‑151 (MQ=255)                                                                                                 aaCC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                               >  1:570984/1‑151 (MQ=255)                                                                                                       tggtgTTTCTCAACCAATAAAA‑CCAA‑TAACAGAAT‑CACATAATCAACCCCCACGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                           <  2:51250/150‑1 (MQ=255)                                                                                                            tATCTCGAACCTTAAAC‑CCGC‑TATCAGGTC‑AATATAACCAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                      <  2:1004010/150‑1 (MQ=255)                                                                                                               ctcGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                  >  2:22305/1‑151 (MQ=255)                                                                                                               cacaaccctttAAC‑CCAC‑TACCACGAC‑CATAAGATAACCCCCCTCGCTAATACCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                    <  2:247672/140‑1 (MQ=255)                                                                                                                  gAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                               >  1:1110478/1‑151 (MQ=255)                                                                                                                        ttAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                         >  1:769383/1‑151 (MQ=255)                                                                                                                          aaa‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                       >  1:566156/1‑151 (MQ=255)                                                                                                                           aacccGC‑TATCAGGAT‑TATATAATCAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                         <  2:1024962/146‑1 (MQ=255)                                                                                                                           aacccGC‑TATAAGGAT‑AAAAGAATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                        <  1:1241784/147‑1 (MQ=255)                                                                                                                           aa‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                      >  2:897804/1‑151 (MQ=255)                                                                                                                            aaccccttATCAGTAT‑AATAGGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACaa                           <  1:193703/140‑1 (MQ=255)                                                                                                                                 c‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 >  1:262391/1‑151 (MQ=255)                                                                                                                                   tATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 >  1:103998/1‑150 (MQ=255)                                                                                                                                      gcagattaatatGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc              <  1:960646/143‑1 (MQ=255)                                                                                                                                        ggAT‑AATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                <  1:1125772/146‑1 (MQ=255)                                                                                                                                        gaaT‑CATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg             <  2:943139/147‑1 (MQ=255)                                                                                                                                         gAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc          >  1:176877/1‑151 (MQ=255)                                                                                                                                              atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt      >  2:299530/1‑151 (MQ=255)                                                                                                                                                atGATCAACCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCa         <  1:946473/146‑1 (MQ=255)                                                                                                                                                atGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc    >  2:228723/1‑151 (MQ=255)                                                                                                                                                atGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc          <  2:101817/145‑1 (MQ=255)                                                                                                                                                 tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc    >  1:250550/1‑150 (MQ=255)                                                                                                                                                  gATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCa         <  2:36653/144‑1 (MQ=255)                                                                                                                                                   aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  1:617282/1‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                CTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAT‑CATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376357‑1376648

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