Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,970 G→A 23.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,9700GA23.1% 393.0 / 107.6 182intergenic (‑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 (69/71);  new base A (33/9);  total (102/80)
Fisher's exact test for biased strand distribution p-value = 7.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-02

Only 100 of 189 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGC  >  NC_000913/3705873‑3706068
                                                                                                 |                                                                                                         
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  1:2747145/1‑100 (MQ=255)
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  1:3035289/1‑100 (MQ=255)
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  2:3461548/1‑100 (MQ=255)
   aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                      <  1:856089/100‑1 (MQ=255)
   aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                      <  1:947869/100‑1 (MQ=255)
     tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                    <  2:71224/100‑1 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                   >  1:1229326/1‑100 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                  >  1:85964/1‑100 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                 >  2:2484952/1‑100 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                 <  2:1888034/99‑1 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                >  1:3673142/1‑100 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                >  2:2011281/1‑100 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              <  1:3240886/100‑1 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                            <  1:564723/99‑1 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                           >  1:2170085/1‑100 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                          >  1:1952097/1‑100 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                          >  1:1994346/1‑100 (MQ=255)
                tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                         <  1:3196072/100‑1 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  2:590897/1‑100 (MQ=255)
                  ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                       <  1:3371274/100‑1 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                      <  2:2199478/100‑1 (MQ=37)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                      <  2:2198857/100‑1 (MQ=37)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                     <  1:1830443/100‑1 (MQ=255)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                     <  2:654093/100‑1 (MQ=37)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                    <  1:1439088/100‑2 (MQ=37)
                     gATGGCGCGGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTGTGGTTATGCCGGATGCGGCGTTATCGCCTTATCCGGcc                                                                                    <  1:1279210/100‑1 (MQ=25)
                       tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                  <  1:1804698/100‑1 (MQ=37)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATCCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                               <  1:2011271/100‑1 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCAGCCTAC‑a                                                                               >  1:363104/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCAGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                              >  1:2173830/1‑100 (MQ=255)
                              gCACGCAGAGGCACGGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                          >  1:2716316/1‑100 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGTATATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                           >  1:1924241/1‑99 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                          >  1:3744004/1‑100 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                         <  2:2170651/100‑1 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                         <  2:573253/100‑1 (MQ=255)
                                aCGCAGAGGCACAGCGCCACCGGGTATACAGACTGTTGAGAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑agaa                                                                            >  1:3529737/1‑96 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AACATCGTg                                                                       >  2:1401492/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                       <  2:2630196/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                       >  1:2084578/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTg                                                                       >  2:3473030/1‑100 (MQ=21)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                      <  2:749898/100‑1 (MQ=255)
                                     gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑caa                                                                  >  1:798010/1‑100 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTTCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAAtt                                                                >  1:2348052/1‑100 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAAtt                                                                >  2:3407449/1‑100 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCa                                                              <  2:857784/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCa                                                              <  2:3103257/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATAATAATTGTTGCCAAAGTGCGCATTGTTTATGCCGGATGCGGCGTTAACGCCCTATCCGGCCTAC‑AAAATCGTGccattt                                                                 >  1:2654461/1‑93 (MQ=17)
                                           caacgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTGTGCCGGATGCGGCGTTAACACCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCa                                                              >  1:2379309/4‑98 (MQ=255)
                                           caGCGTCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAAt                                                            <  2:1266067/100‑1 (MQ=255)
                                             gagcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATAc                                                          <  2:2470253/98‑1 (MQ=255)
                                               gcCACTGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACAt                                                        >  2:1213499/1‑100 (MQ=255)
                                                  aCCGGGTATACACCTTGTTGACAAAGTCCGCTTTGTTTATCTCGGATGCGTCCTTAACGCCTTATCCGTCCTAC‑AAAATCGTGC‑AAATTCATTACATTGc                                                     <  1:313699/100‑1 (MQ=17)
                                                     gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGa                                                >  2:1992425/1‑100 (MQ=255)
                                                       gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑CAATTCAATACATTGC‑‑AGAAt                                              <  1:2326995/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑g                                      <  2:492490/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑g                                      >  1:1904156/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑g                                      >  2:3218815/1‑100 (MQ=255)
                                                                 tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGc                                 >  1:3800681/1‑100 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGcc                                <  2:2951956/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGcc                                >  2:161064/1‑100 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCt                               >  2:1823407/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTg                              >  1:3438451/1‑100 (MQ=255)
                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGAt                            <  2:1501042/100‑1 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATa                           >  1:3196531/1‑100 (MQ=255)
                                                                            tGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCgg                      <  1:1307701/100‑1 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                    <  1:2709829/100‑1 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                   >  1:2611423/1‑100 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGTTAAGCg                       >  2:2423031/1‑97 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                   >  2:2063265/1‑100 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                  >  2:2643108/1‑100 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                  <  1:1419714/100‑1 (MQ=255)
                                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                >  2:301389/1‑100 (MQ=255)
                                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                <  1:373835/100‑1 (MQ=255)
                                                                                   tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt               >  1:1198358/1‑100 (MQ=255)
                                                                                    ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCC‑AAATTCAATATATTGC‑‑AGAAA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATc              >  2:2307935/7‑100 (MQ=18)
                                                                                    ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGC‑AAATTCAATATATTGC‑‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATc              >  2:1720692/7‑100 (MQ=18)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc              >  1:2451771/1‑100 (MQ=255)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc              <  2:1256428/100‑1 (MQ=255)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg            <  2:2716330/100‑1 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAATTC‑TGC‑AAATTCAACAAATTGC‑‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg           >  1:1548395/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGA‑AAATTCAATATATTGC‑‑AGGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg            >  2:2838790/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTAT‑AAATTCAATGTATTGC‑‑AGTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg           >  2:2889443/4‑100 (MQ=25)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGA‑AAATTCAATATATTGC‑‑AGGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg           >  2:1252712/4‑100 (MQ=17)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCTTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg           <  1:3496329/100‑1 (MQ=255)
                                                                                        nggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑ATACCGCGC‑AAATTCAATATATTGCGGAGAAA‑AT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCGTCAGGct         >  1:2043341/3‑99 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGC‑TAATTCAATATATTGC‑‑AGAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc          >  1:431207/3‑100 (MQ=17)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGC‑AAATTCAATATATTGC‑‑AGGAA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc          >  2:172828/3‑100 (MQ=18)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa         >  2:1530808/1‑100 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa         <  2:406035/100‑1 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa         >  1:1379054/1‑100 (MQ=255)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGC‑CAATTCAATATATTGC‑‑AGGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa         >  1:2265583/2‑100 (MQ=18)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt        >  1:2885851/1‑100 (MQ=255)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCC‑AAATTCAATATATTGC‑‑AGAAA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa        >  2:2975010/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGC‑AAATTCAATAAATTGC‑‑A‑AAT‑CAACATGTAGGCCGGATAAGCGCAGCGCATCAGGCa         >  1:2571710/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGC‑TAATTCAATATATTGT‑‑GGGCT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCTCATCAGGCag        >  1:2497197/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGC‑TAATTCAATATATTGC‑‑AGGGA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag        >  1:1217899/1‑99 (MQ=18)
                                                                                            cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGATCGCATCAGGCAtt       >  2:1094665/1‑99 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg   <  1:3025232/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg   <  2:2232917/100‑1 (MQ=255)
                                                                                                tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc  <  2:1952107/100‑1 (MQ=255)
                                                                                                 |                                                                                                         
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGC  >  NC_000913/3705873‑3706068

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