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% 295.9 / 47.2 121F47V (TTT→GTT) yddHflavin reductase like‑protein
Reads supporting (aligned to +/- strand):  ref base A (53/49);  new base C (10/8);  total (63/58)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

Only 100 of 126 total aligned reads displayed.

CAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACG  >  NC_000913/1533614‑1533812
                                                                                                   |                                                                                                   
cAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCaaa                                                                                                     >  1:2237960/1‑100 (MQ=255)
cAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCaaa                                                                                                     >  1:2365975/1‑100 (MQ=255)
 aaCGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAAc                                                                                                    <  1:279008/100‑1 (MQ=255)
 aaCGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAAc                                                                                                    >  1:2858156/1‑100 (MQ=255)
  aCGCCCGGGACAACGATTCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAGTCGCCACGCGTGGAGGTTCAAAc                                                                                                    <  2:3420554/99‑1 (MQ=255)
   cGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTc                                                                                                  <  2:2155966/100‑1 (MQ=255)
       cGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcg                                                                                              <  2:1951278/100‑2 (MQ=255)
       cGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCAcc                                                                                              >  2:2950045/1‑100 (MQ=255)
       cGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCAcc                                                                                              >  1:2490915/1‑100 (MQ=255)
        gggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgc                                                                                             <  1:621216/100‑3 (MQ=255)
         ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcg                                                                                            <  2:2444828/100‑4 (MQ=255)
         ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcg                                                                                            >  2:468461/1‑97 (MQ=255)
         ggATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCgg                                                                                            <  1:1931926/100‑1 (MQ=255)
          gATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgt                                                                                           >  1:3475951/1‑96 (MQ=255)
            tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtgg                                                                                         <  1:1407664/100‑7 (MQ=255)
            tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGTTTGAACCTCCAcgcgtgg                                                                                         <  2:2924988/100‑7 (MQ=255)
            tAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCAt                                                                                         >  1:1327845/1‑100 (MQ=255)
              aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTg                                                                                       <  2:589751/100‑1 (MQ=255)
              aCGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTg                                                                                       <  1:2412363/100‑1 (MQ=255)
               cGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGa                                                                                      <  1:3280228/100‑1 (MQ=255)
               cGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGa                                                                                      >  2:926136/1‑100 (MQ=255)
                gATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAc                                                                                     <  2:2800009/100‑1 (MQ=255)
                 aTGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAcc                                                                                    <  2:376700/100‑1 (MQ=255)
                 aTGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGAcc                                                                                    <  1:926145/100‑1 (MQ=255)
                   gCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAg                                                                                  >  2:407225/1‑100 (MQ=255)
                   gCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAg                                                                                  >  1:114018/1‑100 (MQ=255)
                     cAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGc                                                                                >  1:124293/1‑100 (MQ=255)
                            aCCGTTATGTTCAATCAACTCCCTGGTCGATGTCGTTTTATCTACCATAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGc                                                                          <  2:67149/99‑1 (MQ=255)
                              cGTTATGTCCAATCAACTCCCTGGTCCATTTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCAt                                                                       <  1:249789/100‑1 (MQ=255)
                               gTTATGTTCAATCACCTACCTGGTCAATGTAGGTTTTTCTAAAACAATCGCCACACGTGGAGGGTCAAAGTTTACCGGCATTGAACAGGcggaag                                                                           >  2:3381513/1‑90 (MQ=255)
                               gTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATa                                                                      <  2:3793849/100‑1 (MQ=255)
                                 tATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAAt                                                                    <  2:457301/100‑1 (MQ=255)
                                  aTGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATg                                                                   <  1:122279/100‑1 (MQ=255)
                                    gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                 >  2:2835192/1‑100 (MQ=255)
                                    gTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGtt                                                                 >  1:72313/1‑100 (MQ=255)
                                      tcaatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAc                                                               >  2:124042/1‑100 (MQ=255)
                                       caatcaaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcg                                                              >  2:1531772/1‑100 (MQ=255)
                                            aaCTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCt                                                         >  2:1462397/1‑100 (MQ=255)
                                            aaCTCCCTGGTCCATGTCGATTTATCTACCACAATAGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTAcgccgc                                                          >  2:3637636/1‑99 (MQ=255)
                                                cccTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGa                                                     >  1:975013/1‑100 (MQ=255)
                                                cccTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGCGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTgg                                                       >  2:1599581/1‑98 (MQ=255)
                                                 ccTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGAt                                                    <  2:721595/100‑1 (MQ=255)
                                                 ccTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGAt                                                    >  2:3236034/1‑100 (MQ=255)
                                                    ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                 <  1:552317/100‑1 (MQ=255)
                                                    ggTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGc                                                 >  2:2007510/1‑100 (MQ=255)
                                                     gTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCt                                                <  1:2950066/100‑1 (MQ=255)
                                                        cATGTCGATTTATCTATCACAATCGCCACGCGTGTAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCCCAt                                             <  1:938257/100‑1 (MQ=255)
                                                         aTGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                            >  2:2241085/1‑100 (MQ=255)
                                                           gTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGa                                          >  1:1637688/1‑100 (MQ=255)
                                                           gTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGa                                          <  2:815443/100‑1 (MQ=255)
                                                           gTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGa                                          <  2:3776278/100‑1 (MQ=255)
                                                             cGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCACTGACCAGGCTGCAGCCATAATGTTACGCCGCTGTGATTGCTCATCGaaa                                        >  2:2247036/1‑100 (MQ=255)
                                                              gATTTATCTACCCCGGTCGCCGCGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATc                                            <  2:3722344/95‑1 (MQ=255)
                                                              gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                       >  1:1100763/1‑100 (MQ=255)
                                                              gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                       <  2:243048/100‑1 (MQ=255)
                                                              gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                       >  1:1722181/1‑100 (MQ=255)
                                                              gATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAAc                                       >  1:2537849/1‑100 (MQ=255)
                                                                 ttATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTgg                                    <  1:797228/100‑1 (MQ=255)
                                                                  tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                   >  1:1828001/1‑100 (MQ=255)
                                                                  tATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGt                                   >  2:64162/1‑100 (MQ=255)
                                                                   aTCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGtn                                  <  2:2436661/100‑2 (MQ=255)
                                                                    tCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGa                                 <  1:2036733/100‑1 (MQ=255)
                                                                     cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCGGCTTGGATTGCTCATCGGGACTgatgat                                >  2:3023933/1‑100 (MQ=255)
                                                                     cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                                <  2:3483743/100‑1 (MQ=255)
                                                                     cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                                >  2:2777887/1‑100 (MQ=255)
                                                                     cTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGAt                                >  2:527368/1‑100 (MQ=255)
                                                                      tACCACAATCGCCACGCGTGTAGGTTCAATCTCGGCCGGTATTGGCCAGGCTGGAGCCATAATGTTACGCCGCTGGGATTGCTCCTCGCACCTGGTGATc                               <  1:2224016/100‑1 (MQ=255)
                                                                       aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                              >  2:1303786/1‑100 (MQ=255)
                                                                       aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                              <  1:1135758/100‑1 (MQ=255)
                                                                       aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                              >  1:2670951/1‑100 (MQ=255)
                                                                       aCCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCa                              >  2:2900185/1‑100 (MQ=255)
                                                                         cacaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAg                            >  1:1406435/1‑100 (MQ=255)
                                                                          acaATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGa                           >  1:2797605/1‑100 (MQ=255)
                                                                           caATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGAc                          >  1:2041255/1‑100 (MQ=255)
                                                                              tCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGt                       <  1:1078426/100‑1 (MQ=255)
                                                                                gCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg                     <  1:2605934/100‑1 (MQ=255)
                                                                                gCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg                     <  2:822919/100‑1 (MQ=255)
                                                                                gCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg                     <  2:1619421/100‑1 (MQ=255)
                                                                                aatgccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTg                     >  1:2402933/8‑100 (MQ=38)
                                                                                 ccACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTgg                    >  1:529616/1‑100 (MQ=255)
                                                                                 ccACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTgg                    <  1:1198111/100‑1 (MQ=255)
                                                                                    cgcgTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                 >  1:2862436/1‑100 (MQ=255)
                                                                                    ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTggg                   >  2:3282082/4‑98 (MQ=255)
                                                                                    ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                 >  1:1638781/4‑100 (MQ=255)
                                                                                    ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                 >  1:2471613/4‑100 (MQ=255)
                                                                                    ccggTGGAGTTTGAACCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGcc                 >  1:2753580/4‑100 (MQ=255)
                                                                                          gAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAtt           >  1:1271872/1‑100 (MQ=255)
                                                                                          gAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAtt           <  2:1793240/100‑1 (MQ=255)
                                                                                           aGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAttt          <  1:2459941/100‑1 (MQ=255)
                                                                                           aGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGAttt          <  1:410235/100‑1 (MQ=255)
                                                                                            ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATg              >  2:1068712/1‑95 (MQ=255)
                                                                                            ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa         >  2:3496732/1‑100 (MQ=255)
                                                                                            ggTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTa         <  2:61633/100‑1 (MQ=255)
                                                                                              ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCTAAACTGGTGATCAAGACCGTTGGGCCATGATTTAAc       <  2:2460750/100‑1 (MQ=255)
                                                                                              ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAAc       >  1:3152185/1‑100 (MQ=255)
                                                                                              ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAAc       >  2:2085352/1‑100 (MQ=255)
                                                                                              ttCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAAc       <  2:1221997/100‑1 (MQ=255)
                                                                                               tCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACa      <  1:1231716/100‑1 (MQ=255)
                                                                                                cAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg     <  1:1780062/100‑1 (MQ=255)
                                                                                                 aaaCTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAg     >  2:178315/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
CAACGCCCGGGATAACGATGCCAAATTTACCGTTATGTTCAATCAACTCCCTGGTCCATGTCGATTTATCTACCACAATCGCCACGCGTGGAGGTTCAAACTCCACCGGCATTGACCAGGCTGCAGCCATAATGTTACGCCGCTGGGATTGCTCATCGAAACTGGTGATCAAGACCGTTGGGCCATGATTTAACAGACG  >  NC_000913/1533614‑1533812

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: