Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,970 G→A 25.1% intergenic (‑180/+128) 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,9700GA25.1% 422.0 / 135.7 207intergenic (‑180/+128)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (66/89);  new base A (47/5);  total (113/94)
Fisher's exact test for biased strand distribution p-value = 5.69e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

Only 100 of 217 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705873‑3706069
                                                                                                 |                                                                                                         
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCAAAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  2:3563220/1‑100 (MQ=255)
 ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                        >  2:655458/1‑100 (MQ=255)
 ggAATGAAGTATAGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                        >  2:2461609/1‑100 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                       >  2:2154484/1‑100 (MQ=255)
    aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                     >  2:3315026/1‑100 (MQ=255)
     tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                    >  2:427861/1‑100 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                <  1:280471/100‑1 (MQ=255)
          tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTAACAAAGTGCGATTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                               >  1:836643/1‑100 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              <  1:3643778/100‑1 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              <  1:2523573/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                             <  2:2869453/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                             <  1:982211/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                             <  2:3154222/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                             <  1:2139904/100‑1 (MQ=255)
             tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                            >  1:1022670/1‑100 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCGCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                           <  1:1093105/100‑1 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                           >  2:23952/1‑100 (MQ=255)
                tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                         >  1:3016339/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  2:2066277/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  2:1998442/1‑100 (MQ=255)
                    ggATGGCGCTACACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTCTACGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                     <  1:2461596/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                    <  1:727678/100‑2 (MQ=37)
                      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                   <  1:775701/100‑1 (MQ=255)
                      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                   <  2:382467/100‑1 (MQ=255)
                      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                   <  1:3000303/100‑3 (MQ=25)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                   <  2:1302352/97‑3 (MQ=25)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑aa                                                                             <  2:1804473/100‑1 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑aa                                                                             >  2:3261348/1‑100 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGTGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATc                                                                         >  2:1201301/1‑100 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACATATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                        <  2:505950/100‑1 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGACCGCTGCGACTTTAACGCCTTCTCCGCCAT‑AC‑AAAAACg                                                                        >  1:121753/1‑100 (MQ=21)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTg                                                                      >  1:2574668/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTg                                                                      <  2:1893263/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCACCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATACCGGATGCGGCGTTAACGCCTTATCCGACCT‑AC‑AAAATCGTg                                                                      >  1:1348269/1‑100 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGc                                                                     <  2:2958851/100‑1 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCT‑AC‑AAAATCGTGCa                                                                    >  1:3459258/1‑100 (MQ=255)
                                     gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCaaa                                                                  >  1:2307030/1‑100 (MQ=255)
                                      agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAAt                                                                 <  1:2733075/100‑1 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAAtt                                                                <  1:442051/100‑1 (MQ=255)
                                          acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCaa                                                             >  2:3575216/1‑100 (MQ=255)
                                             gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATAc                                                          <  1:1746054/100‑1 (MQ=255)
                                               gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACAt                                                        <  2:2264695/100‑1 (MQ=255)
                                               gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACAt                                                        >  1:1187350/1‑100 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACAtt                                                       >  1:112978/1‑100 (MQ=255)
                                                 cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCAGGATGTGGCGTTAACTCTTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACAtt                                                       <  1:883397/99‑1 (MQ=18)
                                                  cgatcttataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGc                                                     <  2:1431508/94‑1 (MQ=255)
                                                    cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAg                                                   <  2:3048563/100‑1 (MQ=255)
                                                     gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGa                                                  >  1:3203892/1‑100 (MQ=255)
                                                         ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑Att‑‑‑t                                          <  1:3570604/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCg                                       <  2:1187356/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCg                                       <  1:2832051/100‑1 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGt                                      >  2:3621596/1‑100 (MQ=255)
                                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTa                                     <  1:3621568/100‑1 (MQ=255)
                                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTa                                     <  2:2976643/100‑1 (MQ=255)
                                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTa                                     >  2:414958/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTg                               >  2:282792/1‑100 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGa                              >  1:2021679/1‑100 (MQ=255)
                                                                         aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAg                          >  2:3732283/1‑100 (MQ=255)
                                                                          aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGc                         <  1:1249325/100‑1 (MQ=255)
                                                                          aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGc                         <  1:1661067/100‑1 (MQ=255)
                                                                          aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGc                         <  1:1894304/100‑1 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCgg                       >  2:2369793/1‑100 (MQ=255)
                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGCCCTGATAAGCGGa                      <  2:2538387/100‑1 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑CC‑ACAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAg                     <  1:308987/100‑1 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAg                     >  2:3616599/1‑100 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAgcg                   >  1:372536/1‑100 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAgcg                   >  1:1975033/1‑100 (MQ=255)
                                                                                 tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAgcgc                  >  1:2856861/1‑100 (MQ=255)
                                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                 >  1:2205716/1‑100 (MQ=255)
                                                                                    ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑ACAAAAATCTTGCCAATTCAATATATTGCAGG‑ACC‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATc               >  2:251765/7‑100 (MQ=18)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATc               >  2:1916605/1‑100 (MQ=255)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATc               >  2:2665227/1‑100 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa              <  2:3181627/100‑1 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCa              <  1:1201295/100‑1 (MQ=255)
                                                                                     tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAg             >  1:1468470/6‑100 (MQ=21)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg             <  1:1171033/100‑1 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAg             <  1:2937904/96‑1 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:2937924/5‑100 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAC‑AAAAACCTGCAAATTCAATATATTGCAGG‑ATA‑‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCa              >  2:1528971/5‑100 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑AC‑AAAAACCTGCAAATTCAATATATTGCAGG‑ATA‑‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAg             >  1:494504/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTGCAGT‑TTA‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAgg            >  1:3339126/4‑100 (MQ=21)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACCTGCAAATTCAATATATTGCAGG‑ATA‑‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:1036689/4‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTGCAGT‑TTA‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGc           >  2:761071/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑AC‑AAAAACATGCAAATTCAATAGATTGCAGA‑GAT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:2357275/3‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAATGTGCAAATTCAATAAATTGC‑‑A‑ATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc           >  2:1584065/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCATGAAAATTCAATACATTGCAAG‑ATT‑‑‑TTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:2696508/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:1510802/2‑100 (MQ=17)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  2:17402/1‑100 (MQ=255)
                                                                                          gCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑AC‑AAAAGCACGCAAATTCAATATATTGCAGA‑AAT‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGACaa         >  2:177405/1‑99 (MQ=11)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTGCAGG‑ACC‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:3478606/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTGCACC‑ATC‑‑‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  2:1096657/1‑99 (MQ=21)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTGCACC‑ATC‑‑‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  2:1675406/1‑99 (MQ=21)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAACCTTGCTAATTCAATATATTGCAGG‑GAC‑‑‑TATGTAGGCCTGATAAGCATAGCGCATCAGGCag         >  2:2581391/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACCTGCAAATTCAATATATTGCAGG‑ATA‑‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  2:3526736/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACATGCAAATTCAATAGATTGCAGA‑GAT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  1:2602598/1‑99 (MQ=17)
                                                                                           ccGGATAAGGCGTTTACGCCGCATCTGGCCT‑AC‑AAATTCTTCCAAATTCAATATATTGCAGA‑AAT‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCaat        >  1:1059579/1‑98 (MQ=17)
                                                                                            cGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt       >  2:237122/1‑100 (MQ=255)
                                                                                             ggATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt      <  1:1389415/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  1:1847113/100‑1 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  1:2951913/100‑1 (MQ=255)
                                                                                                 |                                                                                                         
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705873‑3706069

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: ‑