Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,533,713 A→C 15.8% 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.8% 322.4 / 49.5 131F47V (TTT→GTT) yddHflavin reductase like‑protein
Reads supporting (aligned to +/- strand):  ref base A (58/52);  new base C (10/11);  total (68/63)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

AACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAC  >  NC_000913/1533615‑1533811
                                                                                                  |                                                                                                  
aaCGCCCGGGATAACGCTGCCAAATTTACCGTTATGTTCAATCAACGCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAAc                                                                                                   <  2:2821554/100‑1 (MQ=255)
aaCGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAAc                                                                                                   >  2:686495/1‑100 (MQ=255)
aaCGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAAc                                                                                                   <  1:190387/100‑1 (MQ=255)
   gCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTcc                                                                                                >  2:1656833/1‑100 (MQ=255)
     ccGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCAc                                                                                              <  1:3522727/100‑1 (MQ=255)
      cGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcg                                                                                             <  2:306283/100‑2 (MQ=255)
        ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcg                                                                                           <  1:1221368/100‑4 (MQ=255)
        ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCgg                                                                                           <  1:1539342/100‑1 (MQ=255)
        ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCgg                                                                                           <  1:161040/100‑1 (MQ=255)
         gATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgt                                                                                          <  1:3122271/100‑5 (MQ=255)
          aTAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtg                                                                                         <  1:433431/100‑6 (MQ=255)
          aTAACGATGCCAAATTTACCGTTATGCTCAATCAACTCCCTGGTCCATGTCGATTTACCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtg                                                                                         <  2:1220030/100‑6 (MQ=255)
           tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCAt                                                                                        <  2:1898901/100‑1 (MQ=255)
            aaCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtgg                                                                                        <  2:2829295/99‑7 (MQ=38)
             aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTg                                                                                      <  1:3021508/100‑1 (MQ=255)
             aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTg                                                                                      >  1:2042000/1‑100 (MQ=255)
             aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTg                                                                                      <  2:2631975/100‑1 (MQ=255)
             aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTg                                                                                      >  2:2499503/1‑100 (MQ=255)
                aTGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAcc                                                                                   <  1:1146215/100‑1 (MQ=255)
                    cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAACCTCCACCGGCCTTGACCAGGc                                                                               >  1:1152017/1‑100 (MQ=255)
                    cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGc                                                                               >  1:1272316/1‑100 (MQ=255)
                    cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGc                                                                               >  1:2239936/1‑100 (MQ=255)
                    cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGc                                                                               <  1:2592617/100‑1 (MQ=255)
                      aaTTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACAGGCATTGACCAGGCTg                                                                             >  2:2305841/1‑100 (MQ=255)
                         ttACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAg                                                                          <  1:979916/100‑1 (MQ=255)
                         ttACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAg                                                                          >  2:626738/1‑100 (MQ=255)
                          tACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGc                                                                         <  2:3134102/100‑1 (MQ=255)
                           aCCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGcc                                                                        <  1:567314/100‑1 (MQ=255)
                            ccGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCa                                                                       <  2:2661713/100‑1 (MQ=255)
                                 aTGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATg                                                                  <  2:1196293/100‑1 (MQ=255)
                                 aTGTTCAAACAACTCCCTGGTCCATGTCTATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAAt                                                                   <  2:2308301/99‑1 (MQ=255)
                                   gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                >  1:1140733/1‑100 (MQ=255)
                                   gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                >  1:1998361/1‑100 (MQ=255)
                                   gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGACATTGACCAGGCTGCAGCCATAATGtt                                                                >  2:2634779/1‑100 (MQ=255)
                                    ttCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                               <  1:724477/100‑1 (MQ=255)
                                       aatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgc                                                            >  2:199709/1‑100 (MQ=255)
                                       aatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgc                                                            >  1:3739238/1‑100 (MQ=255)
                                           aaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCt                                                        >  2:3249350/1‑100 (MQ=255)
                                            cttCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTg                                                       <  1:492767/98‑1 (MQ=255)
                                             cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCCCGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATATTACGCCGCTgg                                                      >  2:3297542/1‑100 (MQ=255)
                                             cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                                      >  2:754217/1‑100 (MQ=255)
                                             cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                                      >  1:571191/1‑100 (MQ=255)
                                             cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                                      <  2:3811592/100‑1 (MQ=255)
                                              tCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTggg                                                     <  2:2915018/100‑1 (MQ=255)
                                              tCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgag                                                     >  2:492765/1‑98 (MQ=255)
                                               cccTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGa                                                    >  2:1940005/1‑100 (MQ=255)
                                               cccTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGa                                                    <  2:3178092/100‑1 (MQ=255)
                                                   ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                >  2:278451/1‑100 (MQ=255)
                                                   ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                >  2:278076/1‑100 (MQ=255)
                                                    gcccATGTCTATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                <  1:3290291/97‑1 (MQ=255)
                                                       cATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCAt                                            <  1:2385388/100‑1 (MQ=255)
                                                         tGTCGATTTATCTACTACAATCGCCACGCGTTGAGGTTTAAACTCCACCGGCATTAACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCg                                          <  2:1489905/100‑1 (MQ=255)
                                                           tggATTTATCTACCATAATCGCCACGTGGGGAGGTTCAAACTCCGCCGGAATTGACCAGGCTGCGGCCATAATGTTACGCCGCTGGGAt                                                   <  2:3909251/87‑1 (MQ=255)
                                                             gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                      >  1:1933111/1‑100 (MQ=255)
                                                             gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                      >  2:1330765/1‑100 (MQ=255)
                                                             gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                      >  2:3850358/1‑100 (MQ=255)
                                                                ttATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTgg                                   <  1:2578832/100‑1 (MQ=255)
                                                                 tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                  <  1:2834724/100‑1 (MQ=255)
                                                                 tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                  <  1:310327/100‑1 (MQ=255)
                                                                  aTCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTg                                 >  1:2046072/1‑100 (MQ=255)
                                                                   tCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTg                                 <  2:1264269/99‑1 (MQ=255)
                                                                   tCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGa                                >  1:367204/1‑100 (MQ=255)
                                                                    cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                               >  2:939561/1‑100 (MQ=255)
                                                                    cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                               <  2:1933102/100‑1 (MQ=255)
                                                                    cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                               <  1:2004696/100‑1 (MQ=255)
                                                                     tACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATc                              >  2:908423/1‑100 (MQ=255)
                                                                      aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                             <  2:303462/100‑1 (MQ=255)
                                                                      aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                             <  2:571188/100‑1 (MQ=255)
                                                                      aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                             <  1:1865507/100‑1 (MQ=255)
                                                                       ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCaa                            <  1:2869711/100‑1 (MQ=255)
                                                                       ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCaa                            >  1:3440897/1‑100 (MQ=255)
                                                                       ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCaa                            >  1:137245/1‑100 (MQ=255)
                                                                       ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGACATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCaa                            <  2:722586/100‑1 (MQ=255)
                                                                        cacaATCGCCACGCGTTGAGGTTCAAACGCCACCGGCATTGACCAGCCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAg                           <  2:3905727/100‑1 (MQ=255)
                                                                        cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAg                           >  1:3202918/1‑100 (MQ=255)
                                                                         acaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGCAACTGGTGATCTAg                           >  2:3475892/1‑99 (MQ=255)
                                                                         acaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGa                          <  1:1656842/100‑1 (MQ=255)
                                                                           aaTCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGAcc                        <  1:3443370/100‑1 (MQ=255)
                                                                            aTCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCg                       <  2:2840596/100‑1 (MQ=255)
                                                                              caatgccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCCGGGATTGCTCATCGAAACTGGTGATCAAGACCGtt                     >  2:1146660/9‑100 (MQ=38)
                                                                               ggcACGCCTGTAGGTTCCAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTGCGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCATTg                    <  2:3412381/98‑1 (MQ=255)
                                                                               aatgccggTGGAGTTTGAACCTCCAACGGCATTGACCAGGCTGCAGCCATAATATAACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCg                       >  2:2201226/8‑97 (MQ=37)
                                                                                 tgccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTggg                  >  2:3080669/6‑100 (MQ=255)
                                                                                  gccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc                 >  2:620847/5‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  2:840279/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  2:1324382/4‑100 (MQ=255)
                                                                                   ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  1:1473449/4‑100 (MQ=255)
                                                                                        ggAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAt           >  2:1589610/1‑100 (MQ=255)
                                                                                         gAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAtt          >  1:159426/1‑100 (MQ=255)
                                                                                          aGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAttt         <  1:63059/100‑1 (MQ=255)
                                                                                          aGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAttt         <  2:2528141/100‑1 (MQ=255)
                                                                                           ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa        >  1:698156/1‑100 (MQ=255)
                                                                                           ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa        >  1:1094763/1‑100 (MQ=255)
                                                                                           ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGACTTa        >  1:1912999/1‑100 (MQ=255)
                                                                                             ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAAc      <  2:1380718/100‑1 (MQ=255)
                                                                                               cAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg    >  2:205141/1‑100 (MQ=255)
                                                                                               cAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg    >  1:2636872/1‑100 (MQ=255)
                                                                                               cAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg    >  1:3620800/1‑100 (MQ=255)
                                                                                               cAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg    >  2:3771548/1‑100 (MQ=255)
                                                                                               cAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg    >  1:2889664/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: