Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,386,365 T→A 28.6% F169Y (TTC→TAC)  yfbL → putative M28A family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,3650TA28.6% 226.2 / 115.5 151F169Y (TTC→TAC) yfbLputative M28A family peptidase
Reads supporting (aligned to +/- strand):  ref base T (62/45);  new base A (21/22);  total (83/68)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 154 total aligned reads displayed.

GCGGTGTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGAC  >  NC_000913/2386266‑2386462
                                                                                                   |                                                                                                 
gCGGTGTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCctt                                                                                                   >  1:1055770/1‑100 (MQ=255)
     gtgGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCcttccttc                                                                                              <  2:3685525/100‑1 (MQ=255)
     gtgGCAGGATGACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTc                                                                                              <  2:1961540/100‑1 (MQ=255)
      tgGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCg                                                                                             >  1:3052186/1‑100 (MQ=255)
      tgGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCg                                                                                             <  1:995492/100‑1 (MQ=255)
      tgGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCg                                                                                             >  1:3800848/1‑100 (MQ=255)
       ggCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGt                                                                                            >  1:597214/1‑100 (MQ=255)
        gCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTa                                                                                           >  1:3821146/1‑100 (MQ=255)
         cAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAg                                                                                          >  1:582902/1‑100 (MQ=255)
           ggATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAg                                                                                          >  1:804891/1‑98 (MQ=255)
             aTTACTCGAACTGGCACGTTTGTTACAACAGCAAGTACCGAAAACAGGCGTGCAGCTGTCCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGAt                                                                                      <  2:3190851/100‑1 (MQ=255)
               tACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGa                                                                                    <  1:3092222/100‑1 (MQ=255)
                aCTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATg                                                                                     >  2:3472298/1‑98 (MQ=255)
                 cTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGaaa                                                                                  >  1:519990/1‑100 (MQ=255)
                 cTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGaaa                                                                                  <  2:1502192/100‑1 (MQ=255)
                 cTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGaa                                                                                   >  2:1604278/1‑99 (MQ=255)
                 cTCGAACTGGCACGTTTGTTACATAAGCAAGTACCGAAAACGGGCGTGCAGCTGGTCGCATATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGa                                                                                    >  2:500759/1‑98 (MQ=255)
                  tCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTTGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAAt                                                                                 <  2:2945332/100‑1 (MQ=255)
                  tCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAAt                                                                                 <  1:3250067/100‑1 (MQ=255)
                   cGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATg                                                                                >  1:1609060/1‑100 (MQ=255)
                    gAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATgg                                                                               >  2:856469/1‑100 (MQ=255)
                      aCTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATgggg                                                                             <  1:3543514/100‑1 (MQ=255)
                         ggCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAgc                                                                          >  1:2024110/1‑100 (MQ=255)
                           cACGTTTGTTACATCAGCAAGTACCGAAAACCGGCGTGCAGCTGGTCGCCGATGCGACGGAAGAACACCCCTTCTTCCGTAGCGATGAAATGGGGa                                                                            >  2:1430592/1‑96 (MQ=255)
                             cGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCgg                                                                      <  2:1909818/100‑1 (MQ=255)
                             cGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCgg                                                                      >  1:4031460/1‑100 (MQ=255)
                                ttGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgc                                                                   <  2:3156210/100‑1 (MQ=255)
                                ttGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgc                                                                   <  2:1020259/100‑1 (MQ=255)
                                 tGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgca                                                                  <  2:1805462/100‑1 (MQ=255)
                                 tGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgca                                                                  <  1:530014/100‑1 (MQ=255)
                                 tGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgca                                                                  <  2:421709/100‑1 (MQ=255)
                                  gTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgca                                                                  >  2:670795/1‑99 (MQ=255)
                                  gTTACATCAGCAAGTACCGAAAACAGGCGTCCAACTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcat                                                                 >  2:2718720/1‑100 (MQ=255)
                                    tACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgc                                                               >  2:22616/1‑100 (MQ=255)
                                    tACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgc                                                               >  1:1375577/1‑100 (MQ=255)
                                    tACATCAGCAAGTACCGAAAACAGGCATGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgc                                                               <  1:3009209/100‑1 (MQ=255)
                                     aCATCAGCAAGTTCCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCACTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgca                                                              <  2:2016403/100‑1 (MQ=255)
                                     aCATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgca                                                              >  2:2605935/1‑100 (MQ=255)
                                     aCATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgca                                                              >  1:1210802/1‑100 (MQ=255)
                                     aCATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGAGGAAATGGGGAGCGCGGtgcatgca                                                              >  1:1957930/1‑100 (MQ=255)
                                     aCATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgca                                                              <  1:3310268/100‑1 (MQ=255)
                                     aCATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgca                                                              <  1:4124646/100‑1 (MQ=255)
                                        tCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCt                                                           >  1:1469377/1‑100 (MQ=255)
                                         cAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCtt                                                          <  2:2569749/100‑1 (MQ=255)
                                          aGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTc                                                         >  2:1382167/1‑100 (MQ=255)
                                           gCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCg                                                        >  2:2321864/1‑100 (MQ=255)
                                           gCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCg                                                        <  2:2187887/100‑1 (MQ=255)
                                            cAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCCGCTTCGc                                                       >  2:3654631/1‑100 (MQ=255)
                                             aaGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGTAGCGCGGTGCATGCAGCTTCGCt                                                      <  1:3502112/100‑1 (MQ=255)
                                              aGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCtt                                                     <  2:1375568/100‑1 (MQ=255)
                                               gTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTg                                                    >  2:1143343/1‑100 (MQ=255)
                                                 aCCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAg                                                  >  1:213437/1‑100 (MQ=255)
                                                 aCCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAg                                                  <  1:2854859/100‑1 (MQ=255)
                                                  ccGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGc                                                 >  2:3542486/1‑100 (MQ=255)
                                                   cGAAAACAGGCGTGCAGCTGGTCGCGTATGCGTCGGAAGAACCGCCCTGCTTTCGGAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAg                                                  <  2:170589/98‑1 (MQ=255)
                                                     aaaaCAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTc                                              >  2:741770/1‑100 (MQ=255)
                                                     aaaaCAGACGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATAGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTc                                              >  2:2166032/1‑100 (MQ=255)
                                                      aaaCAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTcc                                             <  1:3771550/100‑1 (MQ=255)
                                                      aaaCAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTcc                                             <  1:3355506/100‑1 (MQ=255)
                                                       aaCAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCa                                            <  1:242707/100‑1 (MQ=255)
                                                        aCAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAg                                           <  1:424521/100‑1 (MQ=255)
                                                         cAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGt                                          >  1:3240275/1‑100 (MQ=255)
                                                         cAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGt                                          >  2:2938669/1‑100 (MQ=255)
                                                         cAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGt                                          >  2:1038349/1‑100 (MQ=255)
                                                            gCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGaa                                       >  1:4122147/1‑100 (MQ=255)
                                                            gCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGaa                                       <  2:2086182/100‑1 (MQ=255)
                                                               tGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGCGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAAtt                                    >  1:3641150/1‑100 (MQ=255)
                                                                  aGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAAt                                 <  1:2279205/100‑1 (MQ=255)
                                                                  aGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAAt                                 <  1:441433/100‑1 (MQ=255)
                                                                   gCTGGTCGCCTATGGGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATg                                <  2:1533836/100‑1 (MQ=255)
                                                                   gCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATg                                >  1:640010/1‑100 (MQ=255)
                                                                   gCTGGTCGCCGATGCGTCGGAAGAACCGCCCTTCTGTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATg                                <  2:2727062/100‑1 (MQ=255)
                                                                     tGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTGTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGa                               <  2:574782/99‑1 (MQ=255)
                                                                     tGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGAt                              <  1:249829/100‑1 (MQ=255)
                                                                      ggTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATa                             >  1:1830973/1‑100 (MQ=255)
                                                                      ggTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATa                             >  1:2657883/1‑100 (MQ=255)
                                                                      ggTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATa                             >  1:111072/1‑100 (MQ=255)
                                                                         cGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCa                          >  2:1843596/1‑100 (MQ=255)
                                                                         cGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCa                          >  1:2295216/1‑100 (MQ=255)
                                                                          gCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCAc                         >  1:1861065/1‑100 (MQ=255)
                                                                           ccTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATGAATGATAGCACt                        >  1:1847025/1‑100 (MQ=255)
                                                                            cTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTg                       >  1:2061016/1‑100 (MQ=255)
                                                                            cTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGGGAAATTAATGATAGCACTg                       >  2:3961437/1‑100 (MQ=255)
                                                                              aTGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGga                     <  1:1682409/100‑1 (MQ=255)
                                                                              aTGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGga                     >  1:2883669/1‑100 (MQ=255)
                                                                              aTGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGAGTCAAGTGAAATTAATGATAGCACt                        >  2:379058/1‑97 (MQ=255)
                                                                              aTGCGTCGGAAGAACCGCCCTTCTCTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGga                     <  2:2594634/100‑1 (MQ=255)
                                                                                 cGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgat                  >  1:1120177/1‑100 (MQ=255)
                                                                                 cGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgat                  >  1:1223088/1‑100 (MQ=255)
                                                                                  gTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatg                 <  2:299510/100‑1 (MQ=255)
                                                                                  gTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatg                 >  2:1499194/1‑100 (MQ=255)
                                                                                   tCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatga                >  2:573236/1‑100 (MQ=255)
                                                                                   tCGGAAGAACCGCCAGTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatga                <  2:2179303/100‑1 (MQ=255)
                                                                                   tCGGAAGAACCGCACTTCTTTCGGAGCGATGAAATGGGAAGCGCGGGGCATGCAGCTTCGCTTGAGCGGACAGTGAAATTAATGATAGCACTAga                     >  2:2548883/1‑95 (MQ=255)
                                                                                    cGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatgat               >  2:1609417/1‑100 (MQ=255)
                                                                                          aaCCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTa         <  2:1916129/100‑1 (MQ=255)
                                                                                             cGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTa      >  1:2445246/1‑100 (MQ=255)
                                                                                             cGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTa      >  1:105101/1‑100 (MQ=255)
                                                                                              gCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTGCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTAc     <  2:2260135/100‑1 (MQ=255)
                                                                                              gCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTAc     <  1:741776/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GCGGTGTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTACGAC  >  NC_000913/2386266‑2386462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: