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

CTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAT  >  NC_000913/1376357‑1376650                                                                                                                                                         |                                                                                                                                                  cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACAcc                                                                                                                                                 >  2:650486/1‑151 (MQ=255)  ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACAccc                                                                                                                                                >  2:337639/1‑151 (MQ=255)          gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGa                                                                                                                                         >  2:914717/1‑150 (MQ=255)                  cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                  <  2:305584/149‑1 (MQ=255)                    ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                              >  1:1034646/1‑151 (MQ=255)                      gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                            >  1:649648/1‑151 (MQ=255)                              tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                     >  2:802112/1‑150 (MQ=255)                              tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                     >  2:213963/1‑150 (MQ=255)                               ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                    >  1:706517/1‑150 (MQ=255)                                   ggCGCCCGAAAATTGGTGAGTGTTTCCCTGGGGGTACGGGAACGGGCCGCCGCGCGT‑ACCGGACCG‑GTGGTGTATCTCGCCCCGTAAAC‑CCGC‑TATCAGGATAATATAAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                                  <  2:379001/148‑1 (MQ=255)                                   ggCGCCCGAAAAAGGGGGATTGTTTCCCTGCGCGTACGGAAACGGGTCGCCGCGCGT‑CCCGGACCC‑GTGGTTTATCTCGACCCTTTAAC‑CCGC‑TACCAGGATAATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                <  1:464435/150‑1 (MQ=255)                                    gCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                               >  1:1085513/1‑150 (MQ=255)                                      gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                            >  1:249488/1‑151 (MQ=255)                                        ccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAGACGATCtt                                                                                                          >  1:1109269/1‑151 (MQ=255)                                           aaaaaTGGTGTTTGTTCCCCTGCGCGTCTAGGACCGGTTCGCCGCGCCC‑CCCGGATCG‑GTGGTGTTCCTCGACCCTTAAAC‑CCCC‑TTACAGAACAATTTGACAAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                             <  2:140463/144‑1 (MQ=255)                                            aaaaGGGTGAGTGTTCCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACG‑GTGGTGTATCTCGAACCGTAAAC‑CCGC‑TATCAGAATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                        <  1:617740/149‑1 (MQ=255)                                               agggggagtgaTTCCCTGCGCGTCCTGGAACAGTTC‑CCGGCCGTCCCCGGAACG‑GGGTTGTATCTCGACCTGTTAAC‑CCCT‑TATCAGGATAATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                    <  1:402804/145‑1 (MQ=255)                                               aTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                   >  2:132907/1‑151 (MQ=255)                                                 ggtgggtgTTCCCCTGCGCGTACGGAAACGGGTCGCCGCGCGT‑ACCGGACCG‑GGGGTGTATCTCGACCCGTAAAC‑CCTC‑TATCAGGATAATATGATCACCCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                   <  2:425070/149‑1 (MQ=255)                                                 ggTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                  >  1:1013357/1‑150 (MQ=255)                                                    gagtgaTCCCCGGCGCGTACTGAAACAGGTCGCCGCGCGT‑CCGGGACCC‑GGGGTTTATCCCGACCCGTTAAC‑CCGC‑TATCAGAATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                              <  1:1030753/151‑1 (MQ=255)                                                    gagcgaTTCCTGGCGCGTCCTGAACCAGTGCGCCGCGGGC‑CCGGCATCG‑GGGGTGTATCCCGACCCTTTAAC‑CCGT‑TACCAGGATCATATAACCACCCCCCCCGCCGACGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                              <  2:789939/151‑1 (MQ=255)                                                       tgaTCCCCTGCGCTTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTTTATCTCAACCCGTTAAC‑CCGT‑TATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                            <  1:873015/150‑1 (MQ=255)                                                           tCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                       >  2:1036038/1‑151 (MQ=255)                                                            cccTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                      >  1:403739/1‑151 (MQ=255)                                                            cccTGCGCGTACTGAAACAGGCCGCCGCGCGT‑ACCGAAACC‑GTGGTGTACCTCAAACCGTTAAC‑CCGC‑TATAAGAATAATATGACCACCCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                       <  2:964245/150‑1 (MQ=255)                                                            cccTGCGCGTACGGGAACAGGCCCCCGCGCGT‑ACCGGACCC‑GTGTTGTATCCCGACCCGTAAAA‑CCCC‑TATCAGGATCATATGATCAACCCCCGCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                      <  2:600155/151‑1 (MQ=255)                                                             ccggcgcGTCCTGGAACGGGTCGCCGCGCGT‑CCGGGACCG‑GTGGTGTATCTCGACCCTTAAAC‑CCTC‑TATCAGGATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                     <  2:1212364/148‑1 (MQ=255)                                                              cTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                    >  2:1158610/1‑151 (MQ=255)                                                                 cgcgTCCTGGAAGGGTGCGCCGGGCG‑‑ACGCGAACC‑GGGGTGTTCCCCCAACCCGTAAA‑CCCCTTAAAAGAATAATAGAATCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                     <  2:536115/147‑1 (MQ=255)                                                                  gcttACTGGAACAGGTCGCCGCGCGT‑CCCGGACCC‑GTGGTGTATCTCGAACCGTAAAC‑CCGC‑TATAAGGATAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                 <  2:937084/147‑1 (MQ=255)                                                                     taggggAACAGGCCGCCGCGCGC‑ACCGGACCC‑GGGGTTTATC‑CCAACCCTTAAACCCGT‑TATCAGAATTATATAATAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                  <  2:739093/142‑1 (MQ=255)                                                                                    ccgcGCGT‑CCCGGAACG‑GTGGTGTATCCCGAACCGTAAAC‑CCGC‑TATAAGGATAATATAAACAACACCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                               <  2:718455/150‑1 (MQ=255)                                                                                      gcgcgT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                              >  2:365828/1‑149 (MQ=255)                                                                                        gcgT‑CCGGGACCC‑GTGGTTTATCTCGACCGTTTAAC‑CCGC‑TTTAAGGATCTTTTAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                          <  1:300864/151‑1 (MQ=255)                                                                                        gcgT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                           >  1:438111/1‑150 (MQ=255)                                                                                         cgT‑ACCGGACCC‑GTGGTGTATCCCGACCCTTTAAC‑CCGT‑TATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                           <  1:966956/149‑1 (MQ=255)                                                                                          gT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                           >  1:17427/1‑148 (MQ=255)                                                                                              ccGGACCT‑GTGGTGTATCTCGACCCTTAAAC‑CCGC‑TATCAGGATAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                       <  1:184833/149‑1 (MQ=255)                                                                                              ccGGACCC‑GTGGTGTTTCCCGACCCTTAAAC‑CCCC‑TATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                          <  1:1001266/146‑1 (MQ=255)                                                                                               cGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGAACATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                    >  2:1054980/1‑151 (MQ=255)                                                                                                ggAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                   >  2:940822/1‑151 (MQ=255)                                                                                                 gAACCGGTGTTGTACCCCGAACCGTTAAC‑CCGT‑TATAAGTATAATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                    <  2:699421/150‑1 (MQ=255)                                                                                                  accc‑GTGGTGTATCTCGACCCGTAAAC‑CCGC‑TATAAGGAAAATATGATCAACCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                    <  1:575884/146‑1 (MQ=255)                                                                                                  accc‑GTGGTGTATCTCAACCCTTTAAA‑CCGC‑TATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                 <  1:73406/149‑1 (MQ=255)                                                                                                   aacgggggtgTATCCCGAACCGTAAAC‑CCTC‑TAACAGGCTCATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                     <  1:903523/141‑1 (MQ=255)                                                                                                       ggggtgTATCCCAAACTTTAAAC‑CCGT‑TTTCGGGATAATTTGAACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                             <  1:543724/149‑1 (MQ=255)                                                                                                         ggtgTATCTCGAACCTTTAAC‑CCGC‑TATCAGGATAATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAg                                           <  1:879386/151‑1 (MQ=255)                                                                                                            gtATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                         <  1:592538/150‑1 (MQ=255)                                                                                                            gtATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                        >  2:35384/1‑151 (MQ=255)                                                                                                                 cccgaCCCTTAAACCCGT‑TATCAGAATAATATAAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                         <  2:330438/142‑1 (MQ=255)                                                                                                                   gAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                 >  2:1080160/1‑151 (MQ=255)                                                                                                                       cGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                              >  2:1205061/1‑150 (MQ=255)                                                                                                                          tAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                          >  1:746502/1‑151 (MQ=255)                                                                                                                           aaa‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCACCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                         <  1:692691/151‑1 (MQ=255)                                                                                                                             a‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                        >  1:658967/1‑150 (MQ=255)                                                                                                                                 cattttaactaTCATATAACCAACCCCATCGCCCCTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                         <  2:924148/137‑1 (MQ=255)                                                                                                                                  c‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                   >  1:1133408/1‑151 (MQ=255)                                                                                                                                         ggATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt             <  1:589050/151‑1 (MQ=255)                                                                                                                                           aTCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc            >  1:339693/1‑150 (MQ=255)                                                                                                                                               tatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg       >  2:256664/1‑151 (MQ=255)                                                                                                                                                atGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc      >  1:375549/1‑151 (MQ=255)                                                                                                                                                 tGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg         <  2:659202/147‑1 (MQ=255)                                                                                                                                                    tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat  >  1:539644/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                  CTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGT‑ACCGGAACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAT  >  NC_000913/1376357‑1376650

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