Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 19.9% intergenic (‑195/+113) 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,9850TC19.9% 423.0 / 93.3 182intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (63/82);  new base C (34/2);  total (98/84)
Fisher's exact test for biased strand distribution p-value = 9.69e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

Only 100 of 187 total aligned reads displayed.

TGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705889‑3706084
                                                                                                |                                                                                                          
tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCAGATGCGGCGTTAACGCCTTATc                                                                                                         <  1:1338895/100‑1 (MQ=255)
 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        >  2:3239895/1‑100 (MQ=255)
 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        <  2:2684741/100‑1 (MQ=255)
     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  2:3080071/100‑1 (MQ=255)
     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  1:142860/100‑1 (MQ=255)
      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                   <  1:1020673/100‑1 (MQ=255)
      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                   <  2:2946884/100‑1 (MQ=255)
      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                   <  1:1487344/100‑3 (MQ=37)
         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                               <  1:984224/100‑1 (MQ=255)
              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                          >  2:2689433/1‑100 (MQ=255)
               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                         <  2:1273863/100‑1 (MQ=255)
                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                      <  2:1839423/100‑1 (MQ=255)
                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                      >  1:2278993/1‑100 (MQ=255)
                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                      <  1:2258861/100‑1 (MQ=255)
                     gagGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAATGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                                   >  1:2007920/1‑100 (MQ=255)
                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                                <  1:3087489/100‑1 (MQ=255)
                         cacaGCGCCACCGGGTATACAGATTGTTGACATAGTGCGGTTCGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                               <  1:2346787/100‑1 (MQ=255)
                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                               <  2:1583771/100‑1 (MQ=255)
                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                               <  1:3320623/100‑1 (MQ=255)
                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                               >  1:3197023/1‑100 (MQ=255)
                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                             >  1:956647/1‑100 (MQ=255)
                               gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                         >  2:139067/1‑100 (MQ=255)
                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                        >  1:644144/1‑100 (MQ=255)
                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                      <  2:1691793/100‑1 (MQ=255)
                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                      <  2:1203671/100‑1 (MQ=255)
                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                      <  2:1153347/100‑1 (MQ=255)
                                       gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAt                                                             <  2:3143498/100‑1 (MQ=255)
                                        tatataGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAATTCGTGCAAATTCAATACATTGC‑‑AGAAtt                                                           <  2:1898055/100‑1 (MQ=255)
                                           aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTc                                                       <  2:535277/100‑1 (MQ=255)
                                           aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTc                                                       <  1:1219074/100‑1 (MQ=255)
                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGTCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCg                                                      <  1:1827005/100‑1 (MQ=255)
                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCg                                                      <  2:196868/100‑1 (MQ=255)
                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCg                                                      <  2:1716538/100‑1 (MQ=255)
                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCg                                                      >  2:581642/1‑100 (MQ=255)
                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTa                                                    <  1:2238269/100‑1 (MQ=255)
                                               aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAg                                                   <  2:1527489/100‑1 (MQ=255)
                                                ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAgg                                                  <  2:1027915/100‑1 (MQ=255)
                                                gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATAGC‑‑AGAATTTTCGTAgg                                                  >  1:1956472/7‑100 (MQ=18)
                                                 cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGc                                                 >  1:3043652/6‑100 (MQ=18)
                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGcc                                                >  1:2320545/1‑100 (MQ=255)
                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGcc                                                >  1:1108108/1‑100 (MQ=255)
                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCt                                               <  2:1467663/100‑1 (MQ=255)
                                                    tgACAGAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAGAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTg                                              <  1:1509216/100‑1 (MQ=255)
                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGa                                             >  1:2210422/1‑100 (MQ=255)
                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGAt                                            <  2:160962/100‑1 (MQ=255)
                                                         aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAg                                         >  2:2904942/1‑100 (MQ=255)
                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCg                                       <  2:1547524/100‑1 (MQ=255)
                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCgg                                      <  1:501890/100‑1 (MQ=255)
                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCgg                                      >  1:1679104/1‑100 (MQ=255)
                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAgcg                                  <  1:1908111/100‑1 (MQ=255)
                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAgcg                                  <  1:581642/100‑1 (MQ=255)
                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCa                                >  1:467457/1‑100 (MQ=255)
                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCa                                >  1:101451/1‑100 (MQ=255)
                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCAt                               >  1:2916067/8‑100 (MQ=18)
                                                                   cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑‑ATGAGCCACGTAGGCCTGATAAGCGTAGCGCAt                               >  2:2597289/8‑100 (MQ=17)
                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATc                              >  1:395965/1‑100 (MQ=255)
                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  2:1108111/100‑1 (MQ=255)
                                                                     tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATACGTGTAGGCCTGATAAGCGTAGCGCATCa                             >  2:3311802/6‑100 (MQ=18)
                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg                            >  1:2736833/1‑100 (MQ=255)
                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg                            >  1:2145762/1‑100 (MQ=255)
                                                                      gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAg                            >  1:815794/5‑100 (MQ=17)
                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAACCGCGCAAATTCAATATATTGCGGAGAAAAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:2491171/5‑100 (MQ=18)
                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:2566539/4‑100 (MQ=18)
                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2545684/3‑100 (MQ=18)
                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑‑AGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGc                          >  1:2687291/3‑100 (MQ=18)
                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGC‑‑AGAGATCATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2334551/2‑100 (MQ=17)
                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:2164984/1‑99 (MQ=18)
                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:182183/1‑99 (MQ=18)
                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCTGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  1:1802484/1‑100 (MQ=17)
                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGA‑‑TCGATCCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:1289661/1‑99 (MQ=21)
                                                                           ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGGACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  1:2438299/9‑98 (MQ=18)
                                                                           ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  1:1042575/1‑100 (MQ=255)
                                                                             ngataaggcgTTCACGCCGCATCCAGCCTACAAAATTGTGCAAAATCAATGGATTGC‑‑ACAGCTAGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagatt                     >  1:124068/7‑96 (MQ=37)
                                                                             ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                    >  1:1090501/7‑98 (MQ=21)
                                                                             ggataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑‑AGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                     >  2:2270858/7‑100 (MQ=18)
                                                                              gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                    <  2:382968/100‑1 (MQ=255)
                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttttn                   <  1:869271/100‑2 (MQ=255)
                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:3018944/100‑1 (MQ=255)
                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:139066/100‑1 (MQ=255)
                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:2007934/100‑1 (MQ=255)
                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:2145778/100‑1 (MQ=255)
                                                                                  cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                <  1:3083147/100‑1 (MQ=255)
                                                                                      gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            >  2:1038959/1‑100 (MQ=255)
                                                                                       ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:1630361/100‑1 (MQ=255)
                                                                                        tAACGCCTTATCCGGCCTACAACATCGTGCACATTCATTACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          <  1:2253522/100‑1 (MQ=255)
                                                                                         aaCGCCTTATCCGGCCTACAAAATCGTGCAACTTCAATACATTGC‑‑AGAATTTTCATAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  1:2800636/1‑99 (MQ=255)
                                                                                          aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATTCATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:1531702/100‑1 (MQ=255)
                                                                                          aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:700750/100‑1 (MQ=255)
                                                                                          aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:3057161/100‑1 (MQ=255)
                                                                                           cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       >  2:1990075/1‑100 (MQ=255)
                                                                                           cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAGCa       >  2:511778/1‑100 (MQ=255)
                                                                                            gCCTTATCCGGCCTACAAAATCGTGTAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  2:1634866/100‑1 (MQ=255)
                                                                                            gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  1:2376126/100‑1 (MQ=255)
                                                                                             ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  2:766982/1‑100 (MQ=255)
                                                                                              cTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc    >  2:2402876/1‑100 (MQ=255)
                                                                                               ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  2:2114964/100‑1 (MQ=255)
                                                                                               ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  2:3262999/100‑1 (MQ=255)
                                                                                                tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  2:1939554/1‑100 (MQ=255)
                                                                                                tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  1:2429667/1‑100 (MQ=255)
                                                                                                tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   >  2:656396/1‑99 (MQ=255)
                                                                                                |                                                                                                          
TGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705889‑3706084

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: