Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,533,713 A→C 15.6% 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.6% 374.2 / 64.4 155F47V (TTT→GTT) yddHflavin reductase like‑protein
Reads supporting (aligned to +/- strand):  ref base A (70/60);  new base C (14/10);  total (84/71)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

Only 100 of 157 total aligned reads displayed.

AACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAC  >  NC_000913/1533615‑1533811
                                                                                                  |                                                                                                  
aaCGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAAc                                                                                                   <  1:367684/100‑1 (MQ=255)
     ccGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCAc                                                                                              >  1:2210425/1‑100 (MQ=255)
      cGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCATAATCGCCACGCGTGGAGGTTGAACCTCCAcg                                                                                             <  2:2965536/100‑2 (MQ=255)
       gggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgc                                                                                            <  2:927476/100‑3 (MQ=255)
        ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCgg                                                                                           <  1:1669894/100‑1 (MQ=255)
        ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCgg                                                                                           <  2:4010370/100‑1 (MQ=255)
         gATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgt                                                                                          <  2:1110937/100‑5 (MQ=255)
         gATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGc                                                                                          >  2:3260836/1‑100 (MQ=255)
          aTAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtg                                                                                         <  2:2014376/100‑6 (MQ=255)
          aTAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcacgt                                                                                          <  1:3790125/99‑5 (MQ=255)
           tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtgg                                                                                        <  2:4037267/100‑7 (MQ=255)
           tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCAt                                                                                        <  2:482564/100‑1 (MQ=255)
            aaCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtggc                                                                                       <  2:763712/100‑8 (MQ=38)
             aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTg                                                                                      >  2:615232/1‑100 (MQ=255)
             aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTg                                                                                      >  2:16378/1‑100 (MQ=255)
             aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTg                                                                                      >  2:3371748/1‑100 (MQ=255)
              cGATGCCGAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGa                                                                                     >  2:2310161/1‑100 (MQ=255)
              cGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGa                                                                                     >  1:3058377/1‑100 (MQ=255)
                aTGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCGCCGGCATTGAcc                                                                                   <  2:2465601/100‑1 (MQ=255)
                aTGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAcc                                                                                   <  2:2134316/100‑1 (MQ=255)
                aTGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAcc                                                                                   <  2:119782/100‑1 (MQ=255)
                  gCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAg                                                                                 <  1:3179652/100‑1 (MQ=255)
                  gCCAAATCTATCGCTATGTTCAATCAACTCCCTGGTCCAGGTCGTTGTATCTCCGACAATCGCCACTCGTCGATCTCCAGCCTCCACCCGcctgtaac                                                                                   >  1:2090981/1‑91 (MQ=255)
                    cAAATTTACCGTTATGTTCAATCAACTCCCTGTTCCATGTCGGTTTTTCTACCAAAATCGCCATGCGGGGAGGTTCAAACTCCACCGTCATTGACCAGGc                                                                               >  2:1099666/1‑100 (MQ=255)
                    cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCCAACTCCACCGGCATTGACCAGGc                                                                               >  1:4263087/1‑100 (MQ=255)
                    cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGc                                                                               <  1:3223502/100‑1 (MQ=255)
                    cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGCGGAGGTTCAAACTCCACCGGCATTgcc                                                                                    >  1:2805704/1‑93 (MQ=255)
                     aaaTTTACCGTTATGTTGAATCATCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCt                                                                              <  2:3854531/100‑1 (MQ=255)
                      aaTGTACCGTTATGGTCAATCAACTCCCTGGTCCATGTCAATTTATCTACCACAATCGCCACTCGTTGAGGTTCAAACTCCACCGGCATTAACCAGGCTg                                                                             <  1:2348861/100‑1 (MQ=255)
                        tttACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGATCAAACTCCACCGGCATTGACCAGGCTGCa                                                                           >  1:1732665/1‑100 (MQ=255)
                          tACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGc                                                                         <  1:4015501/100‑1 (MQ=255)
                            ccGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCa                                                                       <  2:4382123/100‑1 (MQ=255)
                                  tGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGt                                                                 <  1:2690862/100‑1 (MQ=255)
                                   gctcaatcaaCTCGCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                >  2:2986062/3‑100 (MQ=255)
                                   gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                >  2:4237330/1‑100 (MQ=255)
                                   gTTCAATCAACTCCCTGGTCCATGTCGATTTATCCACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                <  1:3217007/100‑1 (MQ=255)
                                    ttCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                               <  2:1317877/100‑1 (MQ=255)
                                    ttCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                               >  1:1624613/1‑100 (MQ=255)
                                      caatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcg                                                             >  2:3720017/1‑100 (MQ=255)
                                       aatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgc                                                            >  1:3110453/1‑100 (MQ=255)
                                        atcaaCTCCCTGGTCCATGTAGATTTATCTACGACAATCGCCACGCGTGGAGGTACAAACTCCACCGGCATTCACCAGGCGGCAGCCATAACGTTAcgc                                                            >  1:274481/1‑99 (MQ=255)
                                             cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                                      <  2:1442869/100‑1 (MQ=255)
                                             cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                                      <  1:4153403/100‑1 (MQ=255)
                                              tCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTggg                                                     >  1:547227/1‑100 (MQ=255)
                                              tCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTg                                                       <  1:3898625/98‑1 (MQ=255)
                                                  tGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACGCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTg                                                 <  1:1359849/100‑1 (MQ=255)
                                                   ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                <  2:4343204/100‑1 (MQ=255)
                                                   ccgatctGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                <  2:2890678/94‑1 (MQ=255)
                                                     tCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCt                                               >  2:1327470/1‑99 (MQ=255)
                                                      ccATGTCGTTTCATCTCCCACTATCGCGACGCGTGGATGTTCAAATGCCAACGGGATTGACCAGGCTGCAGTCTTCTTGTTACGCCGCTGGGAttcg                                                <  1:1020172/97‑3 (MQ=255)
                                                        aTGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                           >  2:426143/1‑100 (MQ=255)
                                                         tGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCg                                          <  1:971744/100‑1 (MQ=255)
                                                         tGTCGATTTATCTACCACAATCGCCACGCGAGGGGGTTCAAACTCCACCAGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCg                                          >  1:726902/1‑100 (MQ=255)
                                                            cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGaaa                                       >  1:1014527/1‑100 (MQ=255)
                                                             gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                      <  1:3465013/100‑1 (MQ=255)
                                                             gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                      >  2:2149858/1‑100 (MQ=255)
                                                              aTTTATCTACCACACTCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                      <  1:70294/99‑1 (MQ=255)
                                                              aTTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACt                                     <  1:3904794/100‑1 (MQ=255)
                                                                 tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGCAACTGGt                                  >  1:2385196/1‑100 (MQ=255)
                                                                 tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                  >  2:2465933/1‑100 (MQ=255)
                                                                 tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                  >  2:184309/1‑100 (MQ=255)
                                                                  aTCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTg                                 <  1:1804234/100‑1 (MQ=255)
                                                                  aTCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTg                                 >  2:3173208/1‑100 (MQ=255)
                                                                  aTCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTg                                 <  2:2055008/100‑1 (MQ=255)
                                                                     tACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATc                              <  1:16378/100‑1 (MQ=255)
                                                                     tACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATc                              <  1:795086/100‑1 (MQ=255)
                                                                       ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCaa                            >  1:644727/1‑100 (MQ=255)
                                                                        cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAg                           <  2:3428017/100‑1 (MQ=255)
                                                                        cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAg                           <  2:695059/100‑1 (MQ=255)
                                                                        cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAg                           <  2:1503138/100‑1 (MQ=255)
                                                                          caATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGAc                         <  1:3303648/100‑1 (MQ=255)
                                                                            gtcaatgccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCg                       >  2:4403303/11‑100 (MQ=37)
                                                                            aTCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCg                       >  1:2918114/1‑100 (MQ=255)
                                                                               gCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg                    >  2:2190744/1‑100 (MQ=255)
                                                                               gCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg                    >  1:1939147/1‑100 (MQ=255)
                                                                               gCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg                    >  2:1758043/1‑100 (MQ=255)
                                                                                ccACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTgg                   <  2:1559954/100‑1 (MQ=255)
                                                                                 cACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAACCTGGTGATCAAGACCGTTggg                  >  1:2481797/1‑100 (MQ=255)
                                                                                 cACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTggg                  >  1:4087715/1‑100 (MQ=255)
                                                                                  gccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc                 >  1:1829717/5‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCGTTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAATCTAGTGATCAAGACCGTTGGGc                 >  2:3530247/4‑99 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGCAACTGGTGATCAAGACCGTTGGGcc                >  1:1424668/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  1:2255167/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  2:2856001/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  2:3206869/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  2:1757965/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  1:990615/4‑100 (MQ=255)
                                                                                    gcgTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCa               <  1:826730/100‑1 (MQ=255)
                                                                                     cgTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCAt              >  1:1614117/1‑100 (MQ=255)
                                                                                        ggAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAt           >  1:2437893/1‑100 (MQ=255)
                                                                                           ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGGTCAAGACCGTTGGGCCATGATTTa        >  1:2876217/1‑100 (MQ=255)
                                                                                           ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa        <  2:2422287/100‑1 (MQ=255)
                                                                                           ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa        >  1:3945225/1‑100 (MQ=255)
                                                                                            gTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTaa       >  1:3070794/1‑100 (MQ=255)
                                                                                              tCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTaata     >  1:2680050/1‑98 (MQ=255)
                                                                                              tCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACa     <  1:2190735/100‑1 (MQ=255)
                                                                                               cAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg    <  2:3644448/100‑1 (MQ=255)
                                                                                               cAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg    >  2:2826696/1‑100 (MQ=255)
                                                                                                 aaCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAc  >  1:3728554/1‑100 (MQ=255)
                                                                                                 aaCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAc  >  2:1474700/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
AACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAC  >  NC_000913/1533615‑1533811

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: