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% 331.9 / 52.0 139F47V (TTT→GTT) yddHflavin reductase like‑protein
Reads supporting (aligned to +/- strand):  ref base A (55/62);  new base C (8/13);  total (64/75)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

Only 100 of 144 total aligned reads displayed.

CAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAC  >  NC_000913/1533614‑1533811
                                                                                                   |                                                                                                  
cAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCaaa                                                                                                    >  1:1263654/1‑100 (MQ=255)
cAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCaaa                                                                                                    <  2:644773/100‑1 (MQ=255)
 aaCGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAAc                                                                                                   <  1:3138613/100‑1 (MQ=255)
   cGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTc                                                                                                 >  1:87469/1‑100 (MQ=255)
       cGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCAcc                                                                                             <  1:1240783/100‑1 (MQ=255)
         ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcg                                                                                           <  2:2860764/100‑4 (MQ=255)
          gATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgt                                                                                          <  1:1761454/100‑5 (MQ=255)
          gATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGc                                                                                          <  2:3951886/100‑1 (MQ=255)
           aTAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtg                                                                                         <  1:225555/100‑6 (MQ=255)
           aTAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtg                                                                                         <  2:1505966/100‑6 (MQ=255)
            tAACGCTGCCAAATTTACCGTTATGTTCAATCAACTCCCTTGTCCATGTCGATTTATCTACCACAATCGCCTCGCGTTGAGGTTCAAACTc                                                                                                 >  1:225295/1‑91 (MQ=255)
            tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtgg                                                                                        <  2:3221922/100‑7 (MQ=255)
             aaCGATGCCAAATTTACCGTTATGTTTAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCAtt                                                                                       <  1:3419306/100‑1 (MQ=255)
              accctGCCTAATTTACCGTTTCGTTCAATCCACTCCCTTCTCCATCTCTATTTACCTACCACAATCGCCACGCGTGAAGGTTCAAACTCCACCgg                                                                                           <  1:3751407/91‑1 (MQ=255)
               cGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGa                                                                                     >  2:3481082/1‑100 (MQ=255)
                 aTGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAcc                                                                                   <  2:4135163/100‑1 (MQ=255)
                   gCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGGTCAAACTCCACCGGCATTGACCAg                                                                                 >  2:1631505/1‑100 (MQ=255)
                    ccgggtttACCTTTATGTTTATTCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCTTGCGTGGAGGTTCAAACTTCACCGGCATTGACCAg                                                                                 >  1:2852231/6‑99 (MQ=255)
                    ccAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAgg                                                                                >  2:1281299/1‑100 (MQ=255)
                    ccAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAgg                                                                                >  1:3027931/1‑100 (MQ=255)
                    ccAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAgg                                                                                >  1:371286/1‑100 (MQ=255)
                     cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGc                                                                               <  1:2903688/100‑1 (MQ=255)
                           tACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGc                                                                         >  1:1884236/1‑100 (MQ=255)
                            aCCGTTATGTTCAATCAACTCCCTGTTCTATGTCGATTTATCTACCATAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGc                                                                         <  1:205498/99‑1 (MQ=255)
                             ccGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCa                                                                       >  2:3247732/1‑100 (MQ=255)
                               gTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTCGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATa                                                                     <  1:2406142/100‑1 (MQ=255)
                                   tGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGt                                                                 <  2:1818812/100‑1 (MQ=255)
                                    gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                >  2:1196971/1‑100 (MQ=255)
                                    gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                >  1:1044078/1‑100 (MQ=255)
                                    gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCACGCTGCAGCCATAATGtt                                                                >  1:65376/1‑100 (MQ=255)
                                     ttCAATCGACTCCCTGGTCCATGTCGATTTATCTACCACAATCTCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                               <  1:1386544/100‑1 (MQ=255)
                                     ttCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                               <  2:2451505/100‑1 (MQ=255)
                                     ttCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTa                                                               >  2:3579852/1‑100 (MQ=255)
                                      tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                              >  2:2624590/1‑100 (MQ=255)
                                      tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                              <  1:272248/100‑1 (MQ=255)
                                      tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                              >  2:4127521/1‑100 (MQ=255)
                                      tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                              <  2:1818517/100‑1 (MQ=255)
                                      tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                              <  2:2076279/100‑1 (MQ=255)
                                       caatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcg                                                             >  2:1465657/1‑100 (MQ=255)
                                       caatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcg                                                             <  2:416106/100‑1 (MQ=255)
                                       caatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcg                                                             >  2:50596/1‑100 (MQ=255)
                                       caatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcg                                                             >  2:1845038/1‑100 (MQ=255)
                                         atcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTGcgcc                                                           <  1:2480894/100‑1 (MQ=255)
                                         atcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgcc                                                           >  1:425733/1‑100 (MQ=255)
                                          tcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgccg                                                          >  2:2480363/1‑100 (MQ=255)
                                           caaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgccgc                                                         >  2:939684/1‑100 (MQ=255)
                                            aaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCt                                                        >  2:3901703/1‑100 (MQ=255)
                                              cTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                                      >  1:4294750/1‑100 (MQ=255)
                                               tCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTggg                                                     <  1:3009820/100‑1 (MQ=255)
                                                cccTGGTCCATGTCGATTTATCTACCATAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGa                                                    >  2:4261943/1‑100 (MQ=255)
                                                  cTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGAtt                                                  <  1:1423702/100‑1 (MQ=255)
                                                   tGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTg                                                 >  1:1152945/1‑100 (MQ=255)
                                                   tGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTg                                                 >  2:2997158/1‑100 (MQ=255)
                                                    gttCCATGTCTATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACGGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                <  1:377958/98‑1 (MQ=255)
                                                    ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                <  2:4377944/100‑1 (MQ=255)
                                                    ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                <  2:863921/100‑1 (MQ=255)
                                                      tCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTc                                              <  1:773862/100‑1 (MQ=255)
                                                      tCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTc                                              <  2:846377/100‑1 (MQ=255)
                                                       ccATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCa                                             >  2:1485501/1‑100 (MQ=255)
                                                       ccATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCa                                             >  2:1003199/1‑100 (MQ=255)
                                                         aTGTCGATTTATCTACCACAATCGCCACTTGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                           <  2:2615745/100‑1 (MQ=255)
                                                         aTGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                           >  2:1983898/1‑100 (MQ=255)
                                                           gTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGa                                         <  1:3500608/100‑1 (MQ=255)
                                                             cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGaaa                                       >  2:2210274/1‑100 (MQ=255)
                                                             cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGaaa                                       >  2:2876968/1‑100 (MQ=255)
                                                                tttATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTg                                    <  2:1177052/100‑1 (MQ=255)
                                                                 ttATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTgg                                   <  2:3000278/100‑1 (MQ=255)
                                                                  tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                  <  1:1196966/100‑1 (MQ=255)
                                                                   aTCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTg                                 <  2:1918920/100‑1 (MQ=255)
                                                                   aTCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTg                                 >  2:790536/1‑100 (MQ=255)
                                                                    tCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGa                                >  1:532454/1‑100 (MQ=255)
                                                                     cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                               <  2:3526848/100‑1 (MQ=255)
                                                                     cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                               <  1:1657006/100‑1 (MQ=255)
                                                                      tACCACAATCTCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATc                              <  1:833124/100‑1 (MQ=255)
                                                                       tccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                             <  1:2043854/99‑1 (MQ=255)
                                                                       aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                             <  1:3172114/100‑1 (MQ=255)
                                                                       aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                             <  2:2504484/100‑1 (MQ=255)
                                                                        ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGCCCAGGCTGCAGCCATAATGTTACTCCGCTGGGATTGCTCATCGAAACTGGTGATCa                             >  1:575378/1‑99 (MQ=255)
                                                                        ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCaa                            >  2:2042595/1‑100 (MQ=255)
                                                                        ccACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCaa                            <  1:4395235/100‑1 (MQ=255)
                                                                           caATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGAc                         <  2:1075603/100‑1 (MQ=255)
                                                                           caATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGAc                         <  2:2316260/100‑1 (MQ=255)
                                                                            aaTCGCCACGCGTTGAGGTTCAAACTCCACCGGCATTGACCAGTCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGAc                         <  1:742072/99‑1 (MQ=255)
                                                                             gtcaatgccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCTAAACTGGTGATCAAGACCg                       >  1:421712/11‑100 (MQ=37)
                                                                              tCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGt                      <  1:1636398/100‑1 (MQ=255)
                                                                                gCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg                    <  2:371287/100‑1 (MQ=255)
                                                                                 ccACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCTGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTgg                   >  1:2388729/1‑100 (MQ=255)
                                                                                 ccACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTgg                   <  1:1064385/100‑1 (MQ=255)
                                                                                   gccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc                 >  2:3545743/5‑100 (MQ=255)
                                                                                   aCGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGc                 <  2:4352874/100‑1 (MQ=255)
                                                                                    ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGGTTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  2:3085825/4‑100 (MQ=255)
                                                                                    ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                >  1:4125802/4‑100 (MQ=255)
                                                                                       gTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATg             >  1:2297614/1‑100 (MQ=255)
                                                                                        tGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGa            >  2:1838843/1‑100 (MQ=255)
                                                                                            ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa        <  1:1983890/100‑1 (MQ=255)
                                                                                            ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa        <  1:1281294/100‑1 (MQ=255)
                                                                                               tCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACa     >  1:2517138/1‑100 (MQ=255)
                                                                                               tCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACa     >  1:65974/1‑100 (MQ=255)
                                                                                                 aaaCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGa   >  2:4329828/1‑100 (MQ=255)
                                                                                                  aaCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAc  <  2:2519979/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGAC  >  NC_000913/1533614‑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: