Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,533,713 A→C 15.1% F47V (TTT→GTT)  yddH ← flavin reductase like‑protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,533,7130AC15.1% 384.8 / 58.1 150F47V (TTT→GTT) yddHflavin reductase like‑protein
Reads supporting (aligned to +/- strand):  ref base A (72/55);  new base C (13/10);  total (85/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 152 total aligned reads displayed.

AACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAC  >  NC_000913/1533615‑1533811
                                                                                                  |                                                                                                  
  cGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAAc                                                                                                   <  1:784033/98‑1 (MQ=255)
  cGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTc                                                                                                 >  2:596269/1‑100 (MQ=255)
  cGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTc                                                                                                 >  2:1497420/1‑100 (MQ=255)
  cGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTc                                                                                                 >  1:1545687/1‑100 (MQ=255)
    cccGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCa                                                                                               >  2:3030662/1‑100 (MQ=255)
    cccGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCa                                                                                               >  2:1050244/1‑100 (MQ=255)
      cGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcg                                                                                             <  2:2990676/100‑2 (MQ=255)
      cGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcg                                                                                             <  1:944198/100‑2 (MQ=255)
        ggATTACGAGGCCAAATTTACCGTTATGTTCAATCATCTCCCTGGTCTATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcg                                                                                           <  1:4021529/100‑4 (MQ=255)
        ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcg                                                                                           <  1:2645729/100‑4 (MQ=255)
        ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCgg                                                                                           <  2:314597/100‑1 (MQ=255)
         gATAACGATGCCAAATTTACCGTTATGTTCAATCAGCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgt                                                                                          <  1:1189526/100‑5 (MQ=255)
         gATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgt                                                                                          <  2:574087/100‑5 (MQ=255)
         gATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgt                                                                                          <  1:3611827/100‑5 (MQ=255)
         gATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCGGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTGGAACCTCCAcgcg                                                                                           <  1:1911668/99‑4 (MQ=255)
              cGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGa                                                                                     >  1:2579311/1‑100 (MQ=255)
               gATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAc                                                                                    >  1:2842102/1‑100 (MQ=255)
                 tGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCa                                                                                  <  2:1473688/100‑1 (MQ=255)
                 tGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCa                                                                                  >  1:1511435/1‑100 (MQ=255)
                  gCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAg                                                                                 >  2:4187726/1‑100 (MQ=255)
                   ccAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAgg                                                                                >  1:4153719/1‑100 (MQ=255)
                   ccAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAgg                                                                                <  2:1179030/100‑1 (MQ=255)
                     aaaTTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCAGAATCGCCACGCGTGGAGGTTCAAACTCCAGCGGCATTGACCAgg                                                                                <  1:3598671/98‑1 (MQ=255)
                          tACCGTTATGTTCAATCAACTCCCTGGTCCATGTCTATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGc                                                                         <  1:330812/100‑1 (MQ=255)
                          tACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGc                                                                         <  2:1253042/100‑1 (MQ=255)
                          tACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGc                                                                         <  2:2538971/100‑1 (MQ=255)
                            ccGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCa                                                                       >  1:2395084/1‑100 (MQ=255)
                            ccGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCa                                                                       <  2:3736032/100‑1 (MQ=255)
                             cGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCAt                                                                      >  2:201587/1‑100 (MQ=255)
                               ttATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATaa                                                                    >  2:1877286/1‑100 (MQ=255)
                                tATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCTACCGGCATTGACCAGGCTGCTGCCATAAt                                                                   >  2:2377648/1‑100 (MQ=255)
                                  tGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGt                                                                 >  2:888826/1‑100 (MQ=255)
                                   gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                >  1:709384/1‑100 (MQ=255)
                                    ttCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                               <  2:3971590/100‑1 (MQ=255)
                                    ttCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                               <  1:4187704/100‑1 (MQ=255)
                                     tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                              <  1:1349909/100‑1 (MQ=255)
                                      caatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcg                                                             >  1:2339182/1‑100 (MQ=255)
                                       aatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgc                                                            <  1:4055726/100‑1 (MQ=255)
                                        atcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgcc                                                           >  2:1317627/1‑100 (MQ=255)
                                        atcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgc                                                            >  2:2548357/1‑99 (MQ=255)
                                          caaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgccgc                                                         >  2:3903531/1‑100 (MQ=255)
                                          caaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGATCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgccgc                                                         >  2:2700723/1‑100 (MQ=255)
                                             cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                                      >  1:4069825/1‑100 (MQ=255)
                                                  tGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTg                                                 <  1:1913573/100‑1 (MQ=255)
                                                   ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                >  1:3377994/1‑100 (MQ=255)
                                                   ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                <  1:3126571/100‑1 (MQ=255)
                                                   ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                <  1:1332583/100‑1 (MQ=255)
                                                     tCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACGAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTc                                              <  1:2282291/100‑1 (MQ=255)
                                                     tCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTc                                              <  2:2630371/100‑1 (MQ=255)
                                                      ccATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCa                                             >  1:2197384/1‑100 (MQ=255)
                                                        aTGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                           <  1:704050/100‑1 (MQ=255)
                                                        aTGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                           <  2:322404/100‑1 (MQ=255)
                                                          gTCGATTTATCTACCTCAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGCCCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGa                                         >  1:4165389/1‑100 (MQ=255)
                                                             gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                      >  1:1608349/1‑100 (MQ=255)
                                                              aTTTATCTACCACTATCGCCACGCGTGTCGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATCGCTCATCGAAACt                                     >  1:2640505/1‑100 (MQ=255)
                                                               tttATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACt                                     >  2:4184713/1‑99 (MQ=255)
                                                                 tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                  <  2:3496475/100‑1 (MQ=255)
                                                                 tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                  >  1:1487989/1‑100 (MQ=255)
                                                                 tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                  <  1:3550620/100‑1 (MQ=255)
                                                                 tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                  >  1:345169/1‑100 (MQ=255)
                                                                  aTCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTg                                 >  2:1864685/1‑100 (MQ=255)
                                                                   tCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGa                                >  2:467854/1‑100 (MQ=255)
                                                                    cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                               >  2:3980227/1‑100 (MQ=255)
                                                                    cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                               >  2:845182/1‑100 (MQ=255)
                                                                    cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                               >  1:1610005/1‑100 (MQ=255)
                                                                    cTACCACAATCGCCACGCATGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                               >  2:1753901/1‑100 (MQ=255)
                                                                    cTACCACAATCGCCACGAGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATCATGTTACGCCGCTCGGCTTGCTCATCGAAACTGGTGa                                >  1:1262772/1‑99 (MQ=255)
                                                                      aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                             <  2:1593769/100‑1 (MQ=255)
                                                                      aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                             >  2:2290617/1‑100 (MQ=255)
                                                                       ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCaa                            >  1:1499769/1‑100 (MQ=255)
                                                                       ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCaa                            >  1:457633/1‑100 (MQ=255)
                                                                         acaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGa                          >  2:29924/1‑100 (MQ=255)
                                                                          caATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGAc                         >  2:1106682/1‑100 (MQ=255)
                                                                            aTCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCg                       <  2:2481677/100‑1 (MQ=255)
                                                                               gCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg                    >  1:3860591/1‑100 (MQ=255)
                                                                                ccACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTgg                   <  1:4020099/100‑1 (MQ=255)
                                                                                ccACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTgg                   <  2:2855010/100‑1 (MQ=255)
                                                                                  gccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc                 >  2:2647272/5‑100 (MQ=255)
                                                                                  aCGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc                 <  2:677130/100‑1 (MQ=255)
                                                                                   ncggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  1:3801825/4‑100 (MQ=255)
                                                                                   cgcgTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                <  2:2337061/100‑1 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCACGACCGATGGGcc                >  1:1042798/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  1:322593/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  2:790243/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  1:419506/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  2:489054/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  2:1004008/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  2:3600007/4‑100 (MQ=255)
                                                                                    cggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCa               >  2:3092595/3‑100 (MQ=255)
                                                                                     cgTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCAGAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCAt              <  1:1747456/100‑1 (MQ=255)
                                                                                        ggAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAt           <  1:636951/100‑1 (MQ=255)
                                                                                          aGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAttt         >  1:1076341/1‑100 (MQ=255)
                                                                                           ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa        <  1:1679174/100‑1 (MQ=255)
                                                                                             ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAAc      <  2:1289265/100‑1 (MQ=255)
                                                                                             ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAAc      >  2:2014945/1‑100 (MQ=255)
                                                                                             ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAAc      <  1:596267/100‑1 (MQ=255)
                                                                                             ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAAc      >  1:4120228/1‑100 (MQ=255)
                                                                                             ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAAc      >  2:2622115/1‑100 (MQ=255)
                                                                                                aaaCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGa   >  2:699621/1‑100 (MQ=255)
                                                                                                 aaCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGCTGGGCCATGATTTAACAGAc  <  1:2617737/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
AACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAC  >  NC_000913/1533615‑1533811

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