Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,970 G→A 24.7% 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,9700GA24.7% 508.4 / 153.5 247intergenic (‑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 (85/101);  new base A (54/7);  total (139/108)
Fisher's exact test for biased strand distribution p-value = 1.03e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

Only 100 of 260 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069
                                                                                                   |                                                                                                        
cAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                                                          <  1:922700/100‑1 (MQ=255)
cAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                                                          <  1:2643282/100‑1 (MQ=255)
 aGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                                                                         <  1:2975533/100‑1 (MQ=255)
 aGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgc                                                                                                         <  2:3252070/100‑1 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                        >  1:1904255/1‑100 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGAGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                    >  2:2656873/1‑100 (MQ=255)
       tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                   <  1:1181234/100‑1 (MQ=255)
       tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                   <  2:3080685/100‑1 (MQ=255)
       tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                   <  2:3658811/100‑1 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                  <  1:3152769/100‑1 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                  <  1:2451423/100‑1 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                  <  2:479171/100‑1 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                             >  2:238128/1‑100 (MQ=255)
               tgctgcCGGATGGCGCTGCAGGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                           >  2:2327405/1‑100 (MQ=255)
               tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                           >  2:3728242/1‑100 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                          >  2:1793100/1‑100 (MQ=255)
                gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                          >  1:2105192/1‑100 (MQ=255)
                  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                        <  1:3123439/100‑1 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                       >  1:3082497/1‑100 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                       >  2:241716/1‑100 (MQ=255)
                     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGTTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                     <  1:796663/100‑1 (MQ=37)
                     cGGATGGCGCTGCACGCAGAGGCACAGAGCCACCGGGTACACAGATTGTTGACAAAGTGCACTTTGTTTATTCCGAATGCGGCGTTAACGCCTTATCCg                                                                                      >  2:3456794/1‑99 (MQ=25)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                    <  2:3401360/100‑1 (MQ=37)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                   <  1:3307969/100‑1 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                   <  2:3691514/100‑1 (MQ=255)
                         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                 >  2:376874/1‑100 (MQ=255)
                          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGCTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                <  2:1157721/100‑2 (MQ=25)
                            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                              <  1:3450998/100‑7 (MQ=25)
                            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCGTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaac                                                                              <  1:2034563/100‑4 (MQ=25)
                            cgcTGCACCCAGAGGCACAGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCAGCCGGCATaaa                                                                               <  2:3592721/99‑3 (MQ=21)
                              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                            <  2:5648/100‑1 (MQ=255)
                                gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                          >  1:1083245/1‑100 (MQ=255)
                                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                        <  2:430322/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                       >  1:3435862/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                       >  2:967951/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                       <  1:3714160/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                       <  1:3817292/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACACATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCAGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                       <  2:2290059/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                     <  2:120833/100‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                >  2:1401931/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                               >  1:1099019/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                               <  1:2390026/100‑1 (MQ=255)
                                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                         <  2:110714/100‑1 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                        >  1:1150963/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                        >  1:2515687/1‑100 (MQ=255)
                                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                       <  1:1439973/100‑1 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                      <  2:1616706/100‑1 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                      >  1:1962631/1‑100 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                      <  1:976818/100‑1 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                      <  2:1214979/100‑1 (MQ=255)
                                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑a                                                   <  2:2807935/100‑1 (MQ=255)
                                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGa                                                >  1:446179/1‑100 (MQ=255)
                                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAt                                              <  2:1882066/100‑1 (MQ=255)
                                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAttt                                            >  2:2846995/1‑100 (MQ=255)
                                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAttt                                            >  2:3339301/1‑100 (MQ=255)
                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtttt                                           >  1:2234311/1‑100 (MQ=255)
                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTACGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtttt                                           >  2:180599/1‑100 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑g                                       >  2:2249496/1‑100 (MQ=255)
                                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGc                                  >  2:1122334/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGcc                                 >  1:912497/1‑100 (MQ=255)
                                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGAt                             >  1:1678092/1‑100 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATa                            >  1:1760043/1‑100 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATa                            <  2:2198370/100‑1 (MQ=255)
                                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGc                         >  1:656338/1‑100 (MQ=255)
                                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCgg                       <  1:2170767/100‑1 (MQ=255)
                                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAg                     <  1:3339326/100‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgc                    <  1:2677106/100‑1 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                   >  1:1598521/1‑100 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                   >  1:61407/1‑100 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAACCGGAgcg                   >  2:593317/1‑100 (MQ=255)
                                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcgc                  <  1:734008/100‑1 (MQ=255)
                                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                   >  1:3771752/1‑99 (MQ=255)
                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                 <  1:3634411/100‑1 (MQ=255)
                                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                 <  2:1312236/100‑1 (MQ=255)
                                                                                      ttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGCGCGCATc               >  2:1453080/7‑100 (MQ=17)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg             <  1:1316302/100‑1 (MQ=255)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAAGCGTGCAAATTCAATACATTGC‑‑ATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:1926802/5‑100 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGC‑‑AATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc           >  2:2944901/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGC‑‑AATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg          >  1:1479990/3‑99 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:2557204/3‑100 (MQ=17)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAAT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct          >  2:3467828/3‑99 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc           >  2:2116630/3‑100 (MQ=17)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  1:987211/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa          >  1:3258235/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGC‑‑AGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:2461991/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  1:3785975/2‑100 (MQ=21)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  1:3595745/1‑100 (MQ=255)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTTCAAAATCGTGCAAATTCAACATATTGC‑‑AATTCTCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctg         >  2:1675857/1‑98 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGC‑‑AATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat        >  1:1139482/1‑97 (MQ=21)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:1234710/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  2:3191290/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag         >  1:50490/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  1:858963/1‑100 (MQ=18)
                                                                                            gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  1:1383822/1‑100 (MQ=18)
                                                                                             ccGGATAAGGCGTTCGCGCCGCATCCGGCCTACAAAAGTTTGCAAATTCAATAAATTGC‑‑AGAATTCAT‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAtt        >  1:1516480/1‑100 (MQ=18)
                                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt       >  1:2923840/1‑100 (MQ=255)
                                                                                              cGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGC‑‑AGTTTATTG‑‑GAAGGCCTGATAAGCGCAGCGCATGAGGCaatc       >  1:1367584/1‑97 (MQ=21)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  1:1863991/100‑1 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  2:3802578/100‑1 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  2:1574290/100‑1 (MQ=255)
                                                                                                   |                                                                                                        
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069

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