Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,2710GA55.3% 1.0 / 53.5 47H108H (CAC→CAT) ybcVDUF1398 domain‑containing protein YbcV
Reads supporting (aligned to +/- strand):  ref base G (13/8);  new base A (8/18);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 4.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATC  >  NC_000913/579026‑579471
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                        
gTATTCCCGGATGAACTATTACCACAGCCCTGACGGGGAACTTCTCTGAGGGAGGGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTAGTATGCCAACACCCCGTTGCTGACACGGAAGAAACCGGACGTTATGATTTATCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGGAATGAATTATCAACGGTATAGTAATATCTTTTATGTTCGTgta                                                                                                                                                                                                      <  1:19596/250‑3 (MQ=255)
                          ncccTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATCATTAAAAACGATGCACCAAGGGTTTAGCCCGTACATGTATTGTATTATGCCAACACCCCGGCGCTGACACGGAAGAAACCGGACGTTCTGATTTAGCGTGGAACGTTTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTg                                                                                                                                                                                                        >  1:38657/2‑222 (MQ=38)
                                    ggAACTTATCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGATGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCa                                                                                                                                                                   >  1:362897/1‑249 (MQ=255)
                                               gCGGGAGTGTCCGGGAATAATTAAAA‑CGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTGGCAGTAAATAGTAATGAATTATCGAATGTATAGTAATATCTTTGATGTGCATTGATATTTGTAGCCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTg                                                                                                                                                       <  1:328578/249‑1 (MQ=255)
                                                                       aaaCGATGCACACAGGGTTTAACGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCATGGAAAGATTTGTGTAGTGTTCGGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTACGTTCATGGAGATTGGCAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAAtgtg                                                                                                                                        <  1:434888/241‑1 (MQ=25)
                                                                                acacaGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATTAATTATCAAAGGGATAGTAATATCTTTTATGTTCATTGATATTTGTAACCCATCGTAAAACTCCGGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAAcc                                                                                                                     <  1:66993/251‑1 (MQ=255)
                                                                                   naGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGaaa                                                                                                                                                                                  >  1:109513/2‑187 (MQ=255)
                                                                                         ttAGCGCGTACACGTATTGCATGATGCCAACTCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGATCACAGTAAATAGTAATTAATGATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTAGTAACCCATCGGAAAACTCCTGCGTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCGATCATAg                                                                                                             <  1:315674/250‑1 (MQ=25)
                                                                                           aGCGCGTACTCGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGAc                                                                                                          <  1:402202/251‑1 (MQ=255)
                                                                                           aGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGAc                                                                                                          <  1:141167/251‑1 (MQ=255)
                                                                                                           gcattatGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATc                                                                                                                 <  1:16642/226‑1 (MQ=255)
                                                                                                               tatGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGAGTTCTCTTGATTTCAACNt                                                                                                                    <  1:348450/221‑1 (MQ=255)
                                                                                                                 ngCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGt                                                                                                                                                      >  1:275728/2‑185 (MQ=255)
                                                                                                                   ccAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCtt                                                                                  >  1:94623/1‑251 (MQ=255)
                                                                                                                           ccGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCCACGGTATAGTAATATCTTTTCTGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACATATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAAttt                                                                          >  1:283504/1‑251 (MQ=255)
                                                                                                                                gCTGACACGGAAGNAACCGGACGTTATGATTTAGCGTGGAAAGATNTNTGTAGTGTTCTTAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTNAAGTAATATCNTGTNNATTANTNGNTATTTGTAATCAATCGGAAAACTCCTGCTNTAGCAAGATTTTCCCTGTATTGCNNAANNNTGNNTtctc                                                                                                                                 <  2:334168/191‑1 (MQ=21)
                                                                                                                                gCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGACTGCTCTCAGTAAATAGTGATGAATTATCAAAGGCATAGTAATATCTTTTATGTTCGTGGATATTAGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCACTGTATTGCTGACATGTGATTtctc                                                                                                                                 >  1:334168/1‑191 (MQ=255)
                                                                                                                                 cTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGCATGCTATAAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACAtt                                                                    >  1:298857/1‑251 (MQ=255)
                                                                                                                                   gACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTaacc                                                                                                                                                                                           <  1:421279/130‑3 (MQ=255)
                                                                                                                                     cacGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTg                                                                                                                              >  1:300059/1‑189 (MQ=255)
                                                                                                                                          aagaaACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTNg                                                                                 <  1:347445/229‑1 (MQ=255)
                                                                                                                                                nnnnnnnnnntGNNNNNGCGNGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGt                                                                                                                                                      <  2:275728/144‑1 (MQ=18)
                                                                                                                                                 cGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGAGGCGtg                                                                                        <  1:102563/215‑2 (MQ=255)
                                                                                                                                                 cGGACGCTATGATATAGCGTGGCAAGNTTTGTGNNGANTTCTGAATGCTCTCAGTAANTANNNANNNANTATCCNAGNNATAGTAATANCTTTTATGTTCATGGATATTTGTCACCCATCGGAAAACTCCTGCTTTAGAAAGNTTTTCCCTGNNTTGCTGAAATGTGATTTCTCTNGATNNCNNCNTATNAtn                                                                                                              >  2:102563/1‑192 (MQ=18)
                                                                                                                                                  ggACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATGTACAACCTTTTTAAGTc                                                   <  1:339799/251‑1 (MQ=255)
                                                                                                                                                  ggACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATGTCTTGAGAATTTAACATTTACAACCTTTTTAAGTc                                                   <  1:339804/251‑1 (MQ=255)
                                                                                                                                                   gACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGAGTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAgg                                                                                                            <  1:68288/193‑1 (MQ=255)
                                                                                                                                                      gTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCa                                                                                                                                                                   >  1:202440/1‑135 (MQ=255)
                                                                                                                                                           gATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGANTTTNNNNNTannnnnn                                                                                                                                               >  2:268918/1‑144 (MQ=255)
                                                                                                                                                           gATTTAGCGAGGAAAGATTTGTGTAGTGTTCTGCATGCTCTCAGTAAATAGAAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGAAACACATCGGAAAACTACTGCTTTCGCACGATTTTCCCTGTATTACTGAAATGTGATTTCTCTTGATTTCAACCTATCATAAGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTAAAACATTTTTAAGTCCTTttatt                                           >  1:338421/1‑250 (MQ=25)
                                                                                                                                                               tAGCGTGGAAAGATTTTTGTATTGTTCTGAAGGCTCTCAGTAAATAGTAATGAATTGTCAAAGGTATAGTAATATCTGTTATGTTCATGGATATTTGTAACCCATCGGAAAACGCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAAtg                                                                                                                                          <  1:386201/151‑1 (MQ=255)
                                                                                                                                                                 ncGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGTATTTTCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACgg                                    >  1:82995/2‑251 (MQ=255)
                                                                                                                                                                 gCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACgg                                    >  1:329986/1‑251 (MQ=255)
                                                                                                                                                                     ggAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAAt                                                                                                                                           >  1:390519/1‑144 (MQ=255)
                                                                                                                                                                               gtgtAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCATCCTATCATAGTACGTTTCTATAAGACGCGTGTTTCTTGAGAATTTAACATTTACAGCCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTc                       <  1:198866/250‑1 (MQ=37)
                                                                                                                                                                                      gtTCTNAATGCNCNNAGTAAATAGTAATGAATTANCAAAGGTATAGTAATNTCTTNTATGTNCGNg                                                                                                                                                                                                        >  2:177881/1‑66 (MQ=17)
                                                                                                                                                                                      gtTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATACATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATGTCAACCTAGCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTAc                                                                 <  1:177881/201‑1 (MQ=255)
                                                                                                                                                                                           nnaTGCTCNNAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCa                                                                                                                                                                   <  2:9627/96‑1 (MQ=255)
                                                                                                                                                                                           gAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGAtgt            <  1:377904/249‑1 (MQ=255)
                                                                                                                                                                                             aTGCTCTCAGTAAATAGTCATGAATTATCAAAGGTATCGTAATATCTTTTCTGTTCATGGATATTTCCACACCATCGGAAATCTCCTGCTTTCGCAAGATTTTCGCTGTCTTGCTGTAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATCAGATGc                                                                                           >  1:411547/1‑168 (MQ=255)
                                                                                                                                                                                               gCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACGCCTGCTTTAGCAAGATTTGTCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAATCTATCATATGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACAGTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAata       <  1:227100/250‑1 (MQ=255)
                                                                                                                                                                                                 tctcAGTAAATAGTAATGAATTATCANAGGTATNNTANTATCTTTTATGTTCATGGANATNNNTNNNCNCTCGGNAANNTCCTGCTTTNGCAAGATTTTCCCTGTATTGCTGAAATGTGATTtctc                                                                                                                                 >  2:95622/1‑126 (MQ=18)
                                                                                                                                                                                                 tctcAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTTAAATGTGATTtctc                                                                                                                                 <  1:95622/126‑1 (MQ=255)
                                                                                                                                                                                                    cAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCGAACACAATGTGAATATTATc  <  1:344274/250‑1 (MQ=255)
                                                                                                                                                                                                        aaaTAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTga                                                                                >  1:862/1‑168 (MQ=255)
                                                                                                                                                                                                                  gNATNNTCAANGGTATAGTAATATCTGTTATGTTCNTGGATATTTGTAACCCATCGGAAAACTCCTGNTNTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCANCCTATCATAGNACGNNTCTANANGNTGCGTGTTTCTTGAGAATTTAACATTTACNACCTTTTTAAGTCCTTTTATTAANNCGNNNTTNNCGttt                         <  2:338421/213‑1 (MQ=37)
                                                                                                                                                                                                                                       tatCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTgag                                                                               <  1:261814/138‑1 (MQ=255)
                                                                                                                                                                                                                                            tttATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTAc                                                                 <  1:156978/147‑1 (MQ=255)
                                                                                                                                                                                                                                                    cGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTAc                                                                 >  1:58977/1‑139 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                        
GTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGGAAAACTCCTGCTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCATAGGACGTTTCTATAAGATGCGTATTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATC  >  NC_000913/579026‑579471

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