Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,978 A→C 25.0% intergenic (‑188/+120) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9780AC25.0% 515.5 / 167.6 243intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (69/101);  new base C (47/2);  total (140/103)
Fisher's exact test for biased strand distribution p-value = 2.88e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-03

Only 100 of 247 total aligned reads displayed.

AAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705880‑3706077
                                                                                                  |                                                                                                           
  gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                          <  1:2897167/100‑1 (MQ=255)
  gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                          <  1:280471/100‑1 (MQ=255)
   tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                         <  2:1559957/100‑1 (MQ=255)
    aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                        >  2:1773937/1‑100 (MQ=255)
    aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                        <  1:2523573/100‑1 (MQ=255)
     ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                       <  2:2869453/100‑1 (MQ=255)
      tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                      >  1:1022670/1‑100 (MQ=255)
       gctgcCGGATGACGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                     <  1:1093019/100‑1 (MQ=255)
        ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                    >  2:3712696/1‑100 (MQ=255)
            cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                >  1:803558/1‑100 (MQ=255)
             ggATGGCGCTACACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTCTACGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                               <  1:2461596/100‑1 (MQ=255)
              gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                              <  1:727678/100‑2 (MQ=37)
               aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                             <  2:382467/100‑1 (MQ=255)
                   cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                        <  2:1804473/100‑1 (MQ=255)
                      tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                     >  1:2976622/1‑100 (MQ=255)
                       gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGTGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                    >  2:1201301/1‑100 (MQ=255)
                        cACGCAGAGGCACAGCGCCACCGGGTATACATATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                   <  2:505950/100‑1 (MQ=255)
                         aCGCAGCGCCACAGCGCCACCGGGTCTACAGATTGTTGACAAAGCGCGCATCCGTTTCACCCGATGCGACGTTAACGCCTTATCCCGCCTCC‑AACATc                                                                                    >  1:1735371/1‑98 (MQ=16)
                          cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                 <  2:1893263/100‑1 (MQ=255)
                           gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                <  2:2958851/100‑1 (MQ=255)
                            cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCTAC‑AAAATCGTGCa                                                                               >  1:3459258/1‑100 (MQ=255)
                               agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                            <  1:2733075/100‑1 (MQ=255)
                                   acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                                        <  2:3577927/100‑1 (MQ=255)
                                    caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                       <  2:2370798/100‑1 (MQ=255)
                                    caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                       <  2:401821/100‑1 (MQ=255)
                                        gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                   >  1:1187763/1‑100 (MQ=255)
                                         ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                  >  1:112978/1‑100 (MQ=255)
                                           cgatcttataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                <  2:1431508/94‑1 (MQ=255)
                                             cGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                                           <  1:998186/100‑1 (MQ=255)
                                             cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑g                                                           >  2:1520647/1‑100 (MQ=255)
                                             cGGGTATACAGATTGTTGACAAAGTGCGCTTGGTTTATGCCGGATGAGGAGTTAACGCCTTAGCCATCCTAC‑AAAATCGTGCGAATTGAATACATTGc                                                                >  1:3643028/1‑98 (MQ=17)
                                                 tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑t                                                     >  2:3153657/1‑100 (MQ=255)
                                                     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑g                                               <  2:646015/100‑1 (MQ=255)
                                                     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑g                                               <  2:1187356/100‑1 (MQ=255)
                                                     cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑g                                               <  1:2832051/100‑1 (MQ=255)
                                                      aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑Gt                                              >  2:3621596/1‑100 (MQ=255)
                                                       gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTa                                             <  1:1220134/100‑1 (MQ=255)
                                                       gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTa                                             <  2:2976643/100‑1 (MQ=255)
                                                        aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAg                                            <  1:2373010/100‑1 (MQ=255)
                                                          tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGc                                          >  2:2109402/1‑100 (MQ=255)
                                                             tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTg                                       >  2:282792/1‑100 (MQ=255)
                                                              gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGa                                      >  1:2021679/1‑100 (MQ=255)
                                                                 aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATaa                                   <  1:848009/100‑1 (MQ=255)
                                                                   aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGc                                 <  1:1661067/100‑1 (MQ=255)
                                                                     tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                               >  1:3081429/1‑100 (MQ=255)
                                                                      gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGCCCTGATAAGCGGa                              <  2:2538387/100‑1 (MQ=255)
                                                                       cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                             <  1:3000541/100‑1 (MQ=255)
                                                                        gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                            <  1:1773930/100‑1 (MQ=255)
                                                                         cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                           >  1:372536/1‑100 (MQ=255)
                                                                           ttGTGTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAGAGTGCTAATTAAATACATTGC‑‑T‑GT‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                               <  1:440247/94‑1 (MQ=17)
                                                                             gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                       <  1:1275269/100‑1 (MQ=255)
                                                                              tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      <  2:3181627/100‑1 (MQ=255)
                                                                              tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      <  1:1201295/100‑1 (MQ=255)
                                                                              tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                      >  2:500551/1‑100 (MQ=255)
                                                                              tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  1:3605245/6‑100 (MQ=18)
                                                                              tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                      >  1:961340/6‑100 (MQ=17)
                                                                              tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGC‑‑A‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  1:1468470/6‑100 (MQ=21)
                                                                               ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                     <  1:1171033/100‑1 (MQ=255)
                                                                               gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     <  1:2937904/96‑1 (MQ=17)
                                                                               gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  1:494504/5‑100 (MQ=18)
                                                                                 aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑TA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:3550985/3‑100 (MQ=17)
                                                                                 aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑‑A‑TG‑GG‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:2568609/3‑100 (MQ=17)
                                                                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  1:3307256/3‑100 (MQ=17)
                                                                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:3500583/3‑100 (MQ=18)
                                                                                 aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑AT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                  >  1:1938436/3‑99 (MQ=21)
                                                                                  ggtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑‑AGGA‑AC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  1:1978380/11‑100 (MQ=21)
                                                                                  ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1296614/2‑100 (MQ=17)
                                                                                  ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:2468600/2‑100 (MQ=17)
                                                                                  ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:2696508/2‑100 (MQ=17)
                                                                                  gaccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑‑A‑GA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:3456266/11‑100 (MQ=18)
                                                                                   gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 <  2:3267846/100‑1 (MQ=255)
                                                                                   gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 >  2:17402/1‑100 (MQ=255)
                                                                                   gCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTGC‑‑A‑GA‑GA‑TTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  2:2427896/1‑99 (MQ=14)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑‑A‑‑A‑TA‑TGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                  >  2:2232253/1‑97 (MQ=21)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  2:1675406/1‑99 (MQ=21)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGC‑‑A‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  2:1096657/1‑99 (MQ=21)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑A‑GG‑GA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag                 >  2:2581391/1‑99 (MQ=18)
                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑‑A‑GA‑GA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  1:2602598/1‑99 (MQ=17)
                                                                                    ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  2:746243/9‑98 (MQ=17)
                                                                                    ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                >  2:2291139/9‑98 (MQ=18)
                                                                                    ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑‑A‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                >  1:2325546/9‑98 (MQ=25)
                                                                                    ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑‑A‑AG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                >  1:318330/9‑98 (MQ=17)
                                                                                    ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                >  2:1962587/1‑100 (MQ=255)
                                                                                    ccGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑‑A‑GA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  1:1059579/1‑98 (MQ=17)
                                                                                    ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                >  2:87263/1‑98 (MQ=18)
                                                                                     cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt               >  2:3379908/8‑97 (MQ=18)
                                                                                      ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑‑ATGG‑GC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt              >  2:1349059/7‑100 (MQ=18)
                                                                                      ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtat              >  1:2271363/7‑96 (MQ=17)
                                                                                      ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt              <  1:1389415/100‑1 (MQ=255)
                                                                                        aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg            <  1:1847113/100‑1 (MQ=255)
                                                                                        aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg            <  2:774892/100‑1 (MQ=255)
                                                                                         tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc           >  1:3397427/1‑100 (MQ=255)
                                                                                           cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         >  2:3567587/1‑100 (MQ=255)
                                                                                             gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt       <  2:1649900/100‑1 (MQ=255)
                                                                                              cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  2:2026018/100‑1 (MQ=255)
                                                                                              cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  2:2004902/100‑1 (MQ=255)
                                                                                              cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  1:2109391/100‑1 (MQ=255)
                                                                                              cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAGTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg      <  1:2169642/100‑1 (MQ=255)
                                                                                                 tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa   <  2:1385288/100‑1 (MQ=255)
                                                                                                  aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAGATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt  >  2:2664900/1‑100 (MQ=255)
                                                                                                  |                                                                                                           
AAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705880‑3706077

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