Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 17.9% intergenic (‑258/+50) 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,706,0480GT17.9% 303.7 / 36.3 111intergenic (‑258/+50)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (44/45);  new base T (17/2);  total (64/47)
Fisher's exact test for biased strand distribution p-value = 1.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 115 total aligned reads displayed.

GCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705952‑3706147
                                                                                                  |                                                                                                   
gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgc                                                                                                  >  1:1298443/1‑100 (MQ=255)
   ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCa                                                                                               >  1:1352490/1‑100 (MQ=255)
   atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:802141/9‑100 (MQ=18)
    tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATCCACGTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:2243968/8‑100 (MQ=21)
    tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCAt                                                                                              <  2:1751471/100‑1 (MQ=255)
     gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATc                                                                                             <  2:555853/100‑1 (MQ=255)
     gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATc                                                                                             <  1:1495464/100‑1 (MQ=255)
      tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCa                                                                                            >  1:1741771/1‑100 (MQ=255)
      tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCa                                                                                            <  2:658747/100‑1 (MQ=255)
       ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           <  1:1190594/100‑1 (MQ=255)
       ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           <  1:1481371/100‑1 (MQ=255)
       gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  1:569326/5‑100 (MQ=18)
        tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                          <  1:69866/100‑1 (MQ=255)
         aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:1405477/3‑100 (MQ=17)
         aggCCGGATAAGGCGTTAACGCCGCATCCGGCATACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                         >  1:2129630/3‑100 (MQ=18)
          ggtcggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:1725949/11‑100 (MQ=17)
          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  1:1238634/2‑100 (MQ=25)
          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  1:1651901/2‑100 (MQ=25)
          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:697073/2‑100 (MQ=17)
           gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                       >  1:780720/1‑100 (MQ=255)
           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                       >  1:1545436/1‑99 (MQ=21)
            ccggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                      >  1:1237990/9‑98 (MQ=18)
            ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                      >  2:1319939/9‑98 (MQ=17)
            ccGGATAAGGCGTTCACGTCGCATCGGGTCTACAAAACCATGAAAATTCAATATATTGCAGCAGCTGTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                      >  2:1190519/1‑98 (MQ=14)
              ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                                                   >  2:681972/7‑98 (MQ=21)
              ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                    >  2:816145/7‑100 (MQ=17)
              ggataaggcgTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGcgat                                                                                      >  1:32945/7‑94 (MQ=21)
              ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGATCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:410041/7‑100 (MQ=17)
              ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                    >  2:1806817/7‑100 (MQ=14)
              ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                                    >  2:1623479/7‑96 (MQ=18)
              ggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGCCTGATAAGCGTAGCGAATCAGGCaacta                                                                                    >  1:70169/7‑96 (MQ=21)
              ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                                                    <  2:1907756/100‑1 (MQ=255)
               gataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGCCTGATAAGCGTAGCGCATCAGGCaactac                                                                                   >  2:1801121/6‑95 (MQ=21)
               gataaggcgTTAACGCCGCATACAACCTACAAAAACGTGCAAATTCAATAACTTGCAGCATATACGTAGTCCTGATAAACGGAGCGCCTCAAGCAttttt                                                                                   >  2:1538096/6‑100 (MQ=9)
               gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                                                   <  1:1152144/100‑1 (MQ=255)
                  gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                                >  1:1770233/1‑100 (MQ=255)
                  aaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaattctgtg                                                                                <  2:140165/98‑9 (MQ=25)
                   cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               <  2:1997017/100‑1 (MQ=255)
                   cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  2:1648661/1‑100 (MQ=255)
                    ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                              >  2:33715/1‑100 (MQ=255)
                     gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                                                             >  1:1727334/1‑100 (MQ=255)
                     gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGAGCATCAGGCATTtttgccttt                                                                             >  1:1279050/1‑100 (MQ=255)
                        ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATATGCGTAGCGCATCAGGCATTtttgcctttgcc                                                                          <  1:1554062/100‑1 (MQ=255)
                           aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       <  2:319192/100‑1 (MQ=255)
                           aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       <  2:140808/100‑1 (MQ=255)
                           aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       <  2:604345/100‑1 (MQ=255)
                            ctccTTATCGGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       <  2:1764582/97‑1 (MQ=255)
                             gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                                      >  2:1727823/1‑99 (MQ=255)
                             gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                                                      >  1:468479/1‑99 (MQ=255)
                               cTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                                                   <  2:1987316/100‑1 (MQ=255)
                                 tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                                                 <  1:1220327/100‑1 (MQ=255)
                                  aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  1:1141984/100‑1 (MQ=255)
                                     cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGATCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                             <  2:802222/100‑1 (MQ=255)
                                     cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                             >  1:1435362/1‑100 (MQ=255)
                                     cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                             >  2:1212377/1‑100 (MQ=255)
                                     cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                             <  2:2168804/100‑1 (MQ=255)
                                     cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                             <  1:25037/100‑1 (MQ=255)
                                      ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                            >  1:918291/1‑100 (MQ=255)
                                        ccTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                          <  1:1717279/100‑1 (MQ=255)
                                          tACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                        >  1:371736/1‑100 (MQ=255)
                                           aCAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       >  1:1628180/1‑100 (MQ=255)
                                              aaaTCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                    >  1:2119284/1‑100 (MQ=255)
                                               aaTCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                   <  1:454014/100‑1 (MQ=255)
                                                aTCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                                                  <  1:2243155/100‑1 (MQ=255)
                                                    tGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  2:1123622/100‑1 (MQ=255)
                                                      cAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            >  1:338686/1‑100 (MQ=255)
                                                      cAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            >  1:2005117/1‑100 (MQ=255)
                                                       aaaTTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                           >  1:739278/1‑100 (MQ=255)
                                                        aaTTCAATACATTGCAGAATTTTTGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGTTATCAGTCTTGTATGGCTTTTAATTATTCGATa                                          <  2:1063336/100‑1 (MQ=255)
                                                         aTTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                         <  1:1638917/100‑1 (MQ=255)
                                                           tCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                       <  2:637189/100‑1 (MQ=255)
                                                                 cATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                 >  2:2006118/1‑100 (MQ=255)
                                                                 cATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt                                 >  1:2044313/1‑100 (MQ=255)
                                                                  aTTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt                                >  2:10387/1‑100 (MQ=255)
                                                                  aTTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt                                <  2:1352484/100‑1 (MQ=255)
                                                                    tGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTGGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                              >  2:1246056/1‑100 (MQ=255)
                                                                     gCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  1:1839352/1‑100 (MQ=255)
                                                                         gAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          >  2:1280182/1‑100 (MQ=255)
                                                                           aTTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                        >  1:63250/1‑100 (MQ=255)
                                                                            ttttCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                        <  1:2148615/99‑1 (MQ=255)
                                                                                  gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  >  2:2206982/1‑100 (MQ=255)
                                                                                  gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  <  2:63250/100‑1 (MQ=255)
                                                                                  gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                  <  2:1407615/100‑1 (MQ=255)
                                                                                   tAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGc                 >  2:930808/1‑100 (MQ=255)
                                                                                    aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc                <  1:403495/100‑1 (MQ=255)
                                                                                     ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               <  2:1951618/100‑1 (MQ=255)
                                                                                         tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt            >  2:1255277/1‑99 (MQ=255)
                                                                                         tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  2:2078198/100‑1 (MQ=255)
                                                                                          gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAtag          >  2:1358611/1‑98 (MQ=255)
                                                                                          gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg          >  1:506578/1‑100 (MQ=255)
                                                                                          gATAAGCGGAGCGCATCAGGCACTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg          >  2:1297320/1‑100 (MQ=255)
                                                                                            tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGACACGTTTTCGAAGTGATGTTTGCCTAATGGAt        >  2:637760/1‑100 (MQ=255)
                                                                                             aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc       >  2:294298/1‑100 (MQ=255)
                                                                                             aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc       <  2:65108/100‑1 (MQ=255)
                                                                                              aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGTTCa      >  2:1077256/1‑100 (MQ=255)
                                                                                              aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa      >  2:1940925/1‑100 (MQ=255)
                                                                                              aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa      <  2:821855/100‑1 (MQ=255)
                                                                                               gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa     >  2:649788/1‑100 (MQ=255)
                                                                                               gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa     <  1:10387/100‑1 (MQ=255)
                                                                                                 ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc   <  1:857949/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705952‑3706147

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: