Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,970 G→A 29.4% 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,9700GA29.4% 474.9 / 221.5 276intergenic (‑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 (98/97);  new base A (76/5);  total (174/102)
Fisher's exact test for biased strand distribution p-value = 2.00e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-02

Only 100 of 291 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069
                                                                                                   |                                                                                                       
cAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                                                         >  1:3979227/1‑100 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                       >  2:3560538/1‑100 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                       <  1:3449022/100‑1 (MQ=255)
   ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                      >  2:1842173/1‑100 (MQ=255)
     aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                    <  2:895158/100‑1 (MQ=255)
       tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                  <  2:148511/100‑1 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                 <  1:1814597/100‑1 (MQ=255)
        gAAGGATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                 <  1:1359184/100‑1 (MQ=255)
         aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                <  2:1540494/100‑1 (MQ=255)
          aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                               >  2:3772074/1‑100 (MQ=255)
          aGTATTGATGCCGGCTGGCGCTGCACGCAGAGGCACAGCGCCCCCGGGTATACAGATTGTTGACAAGGATAGCTTTGTTTATACGGACTGCGGCGTTaa                                                                                                >  1:1876403/1‑99 (MQ=255)
           gTATTGCTGCCGGATGGCGCTGCACGTAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                              <  1:1825075/100‑1 (MQ=255)
           gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                              <  2:3619328/100‑1 (MQ=255)
            tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                             >  2:2251846/1‑100 (MQ=255)
              ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                           <  1:2338699/100‑1 (MQ=255)
               tgctgcCGGATGGCGCTGCGCGCAGGGGCACAGCGCGACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                           <  1:959391/99‑1 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                         >  1:1662104/1‑100 (MQ=255)
                 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                        >  2:1117413/1‑100 (MQ=255)
                 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                        >  2:2558601/1‑100 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                      >  2:1480318/1‑100 (MQ=255)
                      ggATGGCGCTGTACGCAGAGGCACAGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                   <  1:1011624/100‑1 (MQ=25)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                   <  1:2156974/100‑1 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                  <  1:3486378/100‑2 (MQ=37)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                  <  2:4163479/99‑1 (MQ=255)
                          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                               <  1:2259561/100‑1 (MQ=255)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                              >  2:552468/1‑100 (MQ=255)
                            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                             >  2:2459558/1‑100 (MQ=255)
                             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                            >  1:1625780/1‑100 (MQ=255)
                             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                            >  1:3112902/1‑100 (MQ=255)
                              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTACGCCGGATGCGGCGTTAACGCCTTACCCGGCCTACaaaa                                                                           <  1:2716160/100‑1 (MQ=255)
                               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                          >  2:1259397/1‑100 (MQ=255)
                               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                          <  1:2884159/100‑1 (MQ=255)
                               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTAGGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaa                                                                          <  2:4276910/100‑10 (MQ=21)
                                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                         >  1:437807/1‑100 (MQ=255)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                        <  1:3231746/100‑1 (MQ=255)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                        <  1:410782/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGCCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTCATCCGGCCTACAAAATCGTg                                                                      >  2:2690664/1‑100 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                     <  2:2045090/100‑1 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                     <  2:300097/100‑1 (MQ=255)
                                      agagGCACAGCGCCACCGGGTATACAGATGGTTGACAAAGTGCGCCTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaa                                                                   >  1:540435/1‑100 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                               >  2:4059046/1‑100 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                               >  2:3054584/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                              >  2:1326473/1‑100 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                           <  1:1556991/100‑1 (MQ=255)
                                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                          <  1:1349932/100‑1 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                       >  2:3730608/1‑100 (MQ=255)
                                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                      <  2:2727898/100‑1 (MQ=255)
                                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑Ga‑a                                               >  1:1493618/1‑100 (MQ=255)
                                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAAATTGCA‑Ga‑a                                               >  2:3682449/1‑100 (MQ=255)
                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACCAAATCGTGCAAATTCAATACATTGCA‑GA‑Attt                                            >  2:2488582/1‑99 (MQ=255)
                                                              cAGATTTTTGACAAAGTGCGCTTTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑Atttt                                           <  1:3551303/98‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑g                                       <  2:2125560/100‑1 (MQ=255)
                                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑Gt                                      >  1:2812120/1‑100 (MQ=255)
                                                                  gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAgg                                   >  1:3687242/7‑100 (MQ=18)
                                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGc                                  >  2:2514801/1‑100 (MQ=255)
                                                                   cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGc                                  >  2:2929313/6‑100 (MQ=18)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCTGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGcc                                 <  1:3843477/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGcc                                 >  2:2644443/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGcc                                 <  2:1356045/100‑1 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTg                               >  1:2002146/1‑100 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATa                            <  1:3562274/100‑1 (MQ=255)
                                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATaa                           >  1:254817/1‑100 (MQ=255)
                                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCgg                       >  1:4209601/1‑100 (MQ=255)
                                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCgg                       <  2:2403523/100‑1 (MQ=255)
                                                                                gattttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAg                     <  2:1141612/97‑1 (MQ=255)
                                                                                ctcTTTGTTTATGCCAGATGCGGCTTCAACGTCTTATCCTGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAg                     <  1:407485/98‑1 (MQ=17)
                                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAg                     <  1:7141/100‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAgc                    <  2:1756668/100‑1 (MQ=255)
                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                 <  2:3966997/100‑1 (MQ=255)
                                                                                    ctttgtgggtcGGATAAGGCGTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCA‑GG‑AACGAT‑‑GTAGGCCTGATAAGCGTAGCGCa                 >  1:2015075/11‑100 (MQ=21)
                                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GT‑TTATTG‑‑GTAGGCCTGATAAGCGCAGCGCAt                >  2:449426/8‑100 (MQ=21)
                                                                                     cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATc               >  2:4157011/8‑100 (MQ=18)
                                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa              <  2:3407641/100‑1 (MQ=255)
                                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa              <  2:977857/100‑1 (MQ=255)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCA‑TG‑GGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  1:2607893/5‑100 (MQ=17)
                                                                                        gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg             >  1:1720797/5‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAATCGTGCTAATTCAATATATTGCA‑GA‑AACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:4184301/4‑100 (MQ=17)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GG‑AGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:2083808/4‑100 (MQ=17)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCA‑AT‑CTCCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:1048740/4‑100 (MQ=17)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAACTTTGCAAATTCAATATATTGC‑‑‑A‑ATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:3139000/4‑100 (MQ=18)
                                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg            >  2:93083/1‑100 (MQ=255)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGG‑AAATTC‑‑‑‑AGGCCTGATAAGCGTAGCGCATCAGGct          >  2:1468175/3‑99 (MQ=25)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCA‑GAGATTAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:2828031/3‑100 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTACGCCTGCTAAGCGGAgcccctc               >  1:819110/3‑91 (MQ=17)
                                                                                           tGCCGGATGCGGCGTTAACTACATATCCGGCCTACAAAATCGAGCAAATTCAATACATTGCA‑GA‑ACTTTC‑‑GTACGCCTGATAAGCGAAGCGCAACAGGc           >  1:3194132/1‑99 (MQ=17)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  1:1295317/1‑100 (MQ=255)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  2:3325062/1‑100 (MQ=255)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  2:519332/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CC‑ATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:753226/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GG‑AGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:63589/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GG‑AGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:3386122/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCA‑AT‑CTCCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:614531/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑A‑ATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:2769529/2‑100 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAATTTCAATATATTGCA‑GA‑GATTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  1:3056236/1‑99 (MQ=14)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTG‑GG‑CTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:4198592/1‑99 (MQ=21)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GT‑TTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa         >  2:4355208/1‑99 (MQ=25)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CC‑ATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  1:1085725/1‑99 (MQ=21)
                                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt      >  2:1594833/1‑100 (MQ=255)
                                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt      <  2:1972368/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  2:2364972/100‑1 (MQ=255)
                                                                                                   |                                                                                                       
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069

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