Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 26.8% intergenic (‑194/+114) 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,9840TG26.8% 352.3 / 147.8 168intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (60/63);  new base G (43/2);  total (103/65)
Fisher's exact test for biased strand distribution p-value = 2.91e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

Only 100 of 170 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTC  >  NC_000913/3705885‑3706082
                                                                                                   |                                                                                                      
ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                        <  1:2509244/100‑1 (MQ=255)
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                      <  1:1498460/100‑1 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                    <  2:566474/100‑1 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                    >  1:696283/1‑100 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                    >  2:115415/1‑100 (MQ=255)
      ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                  <  1:2852104/100‑1 (MQ=255)
       cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                 <  1:766363/100‑1 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                               <  1:999751/100‑2 (MQ=37)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                              <  2:2245566/100‑1 (MQ=255)
             gcgtTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                          <  2:1280378/100‑1 (MQ=255)
             gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                          <  1:496771/100‑1 (MQ=255)
             gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                          >  1:926053/1‑100 (MQ=255)
               gcTTCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGAATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                        <  2:724490/100‑1 (MQ=255)
                cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                                       >  1:2132347/1‑100 (MQ=255)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                      <  2:1324200/100‑1 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                     >  1:1004754/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                    <  2:1202658/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                  >  2:1334618/1‑100 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                  >  2:1457113/1‑100 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGTCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                 <  2:1370391/100‑1 (MQ=255)
                         gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCGTTGTTTATGCCGGATACGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                              >  2:2998643/1‑100 (MQ=255)
                          agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                             <  1:2757140/100‑1 (MQ=255)
                           ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                            >  1:1995156/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                          >  2:269383/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                          >  2:2913057/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                          >  1:556115/1‑100 (MQ=255)
                              acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                          <  1:451369/99‑1 (MQ=255)
                                aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                       <  2:2414643/100‑1 (MQ=255)
                                aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                       <  2:93137/100‑1 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                   <  1:147234/100‑1 (MQ=255)
                                     cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                  <  1:674036/100‑1 (MQ=255)
                                        cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                              >  1:2685563/1‑100 (MQ=255)
                                         gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga                                                             <  2:3058533/100‑1 (MQ=255)
                                            tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑t                                                        >  1:1843147/1‑100 (MQ=255)
                                               aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTc                                                     <  2:1673182/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGt                                                   >  2:22436/1‑100 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGt                                                   >  2:829813/1‑100 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTa                                                  <  1:2061772/100‑1 (MQ=255)
                                                    ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GC‑AT‑TTTTGTAgg                                                >  1:1179582/1‑100 (MQ=255)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑tt                                                       >  2:2555142/1‑91 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCGTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCt                                             >  2:767657/1‑99 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTg                                            >  1:2236491/1‑100 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCGTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GT‑AT‑TTTCGTAGGCCt                                             >  2:1258312/1‑99 (MQ=255)
                                                         gccAAAGTGCGCTTCGTTTATGCCGGATGCGGAGTTAACGCCTTATCGGGCCCGC‑AAAATCGTGCCAGTTCCTTCCATTGCA‑GC‑CT‑TTTCATCGGcc                                              >  2:2754636/3‑97 (MQ=9)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGc                                      <  1:863583/100‑1 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCgg                                    <  1:759182/100‑1 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCgg                                    <  1:2861921/100‑1 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAg                                  <  1:2396036/100‑1 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAg                                  >  2:2482128/1‑100 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgc                                 >  2:2254556/1‑100 (MQ=255)
                                                                   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgc                                 <  2:791689/100‑1 (MQ=255)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                                <  1:538876/100‑1 (MQ=255)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                                >  2:2401238/1‑100 (MQ=255)
                                                                     tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgcgc                               >  2:566416/1‑100 (MQ=255)
                                                                      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCa                              <  1:2076635/100‑1 (MQ=255)
                                                                       tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCAt                             >  1:114102/8‑100 (MQ=17)
                                                                       cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑TGACGTAGGCCTGATAAGCGTAGCGCATc                            >  1:471439/8‑100 (MQ=18)
                                                                        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATc                            >  1:596718/7‑100 (MQ=18)
                                                                          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                          <  2:2419844/100‑1 (MQ=255)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAg                          >  2:2800360/5‑100 (MQ=21)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAg                          >  1:1050866/5‑100 (MQ=17)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑AT‑ACGTGTAGGCCTGATAAGCGTAGCGCATCAg                          >  1:1975596/5‑100 (MQ=18)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑GA‑TTATGTAGGCCTGATAAGCGTAGCGCATCAg                          >  2:1903096/5‑100 (MQ=17)
                                                                          gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                          >  2:705336/5‑100 (MQ=18)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                         >  2:322928/1‑100 (MQ=255)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCA‑TG‑GG‑CCATGTAGGCCTGATAAGCGTAGCGCATCAGGc                        >  2:1996690/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTAACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCAGAACGCATCa                           >  1:1765018/3‑97 (MQ=17)
                                                                            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                        >  1:702360/1‑100 (MQ=255)
                                                                            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                        <  1:2919397/100‑1 (MQ=255)
                                                                             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                       <  1:2658495/100‑1 (MQ=255)
                                                                             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGACTGATAAGCGGAGCGCATCAAGCa                       >  2:1884736/1‑100 (MQ=255)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:50535/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                       >  1:2474904/2‑100 (MQ=18)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                      >  2:62791/1‑100 (MQ=255)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                      >  2:536358/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GA‑GA‑TTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                      >  2:2520172/1‑99 (MQ=14)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                      >  1:1727568/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                      >  1:1248746/1‑99 (MQ=25)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:2154051/9‑98 (MQ=18)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                     >  1:2688128/9‑98 (MQ=25)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                     >  2:2864190/9‑98 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                  >  2:1850411/7‑98 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑TCTGGTAGGCCTGATAAGCATCGCGCATCAGGCaagct                   >  2:2666410/7‑96 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑TTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacgc                    >  2:205698/7‑96 (MQ=17)
                                                                                 ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGATCGCATCAGGCAtttt                   >  1:3058388/1‑100 (MQ=255)
                                                                                  gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                  >  1:261136/6‑100 (MQ=17)
                                                                                  gataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                  >  1:2372691/6‑100 (MQ=18)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                 <  1:1431348/100‑1 (MQ=255)
                                                                                    taaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                >  1:1566176/4‑100 (MQ=18)
                                                                                     gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc               >  1:818485/1‑100 (MQ=255)
                                                                                      aggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑GA‑CTATGTAGGCCTGATAAGCATAGCGCATCAGGCagctttacgt              >  1:1349939/2‑91 (MQ=21)
                                                                                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt             >  2:40406/1‑100 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttn           <  2:1755328/100‑2 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg           <  1:570961/100‑1 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc          >  1:1927671/1‑100 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc          <  1:322930/100‑1 (MQ=255)
                                                                                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTCTTTGCCATc      >  2:2061358/1‑99 (MQ=255)
                                                                                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg    >  2:89427/1‑100 (MQ=255)
                                                                                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAtagat   >  1:3072196/1‑96 (MQ=255)
                                                                                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg    <  1:2710850/99‑1 (MQ=255)
                                                                                                   |                                                                                                      
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTC  >  NC_000913/3705885‑3706082

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