Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,028 T→C 23.3% intergenic (‑238/+70) 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,0280TC23.3% 609.4 / 192.9 324intergenic (‑238/+70)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (139/103);  new base C (71/3);  total (217/107)
Fisher's exact test for biased strand distribution p-value = 1.53e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 326 total aligned reads displayed.

TACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTG  >  NC_000913/3705931‑3706124
                                                                                                    |                                                                                                      
taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                                                                  <  1:2567488/100‑1 (MQ=255)
  cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑g                                                                                              <  2:799482/100‑1 (MQ=255)
   aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑Gt                                                                                             >  1:1320482/1‑100 (MQ=255)
    gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTa                                                                                            >  1:3090983/1‑100 (MQ=255)
     aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAg                                                                                           <  1:2286545/100‑1 (MQ=255)
       tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGc                                                                                         <  1:860543/100‑1 (MQ=255)
        cttccgatcttGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGcc                                                                                        <  2:1856683/90‑1 (MQ=255)
           gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGa                                                                                     <  1:1805404/100‑1 (MQ=255)
                aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGc                                                                                >  2:3349954/1‑100 (MQ=255)
                  tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                                                              <  1:2174010/100‑1 (MQ=255)
                     gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                            <  1:1042068/99‑1 (MQ=255)
                      cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGACTT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                           >  2:1645393/1‑99 (MQ=255)
                        ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                        <  2:36542/100‑1 (MQ=255)
                         tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                       >  1:2618801/8‑100 (MQ=17)
                           tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                     <  2:1601408/100‑1 (MQ=255)
                           tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                    >  1:1102970/6‑100 (MQ=18)
                            gaaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTT‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAg                                                                    >  1:2830811/5‑100 (MQ=18)
                             taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                   >  2:2482918/4‑100 (MQ=17)
                             tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                    <  2:1484001/99‑1 (MQ=255)
                              aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                  >  2:1979436/3‑100 (MQ=18)
                              aggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                  >  1:3366364/3‑100 (MQ=17)
                              aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                  >  1:1359021/3‑100 (MQ=21)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                >  1:3628725/2‑100 (MQ=18)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATC‑‑‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                >  1:864509/2‑100 (MQ=18)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  2:1747157/2‑100 (MQ=18)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  2:3765236/2‑100 (MQ=17)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  1:75441/2‑100 (MQ=17)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                  >  1:3585531/2‑100 (MQ=21)
                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  2:3096200/2‑100 (MQ=21)
                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  1:472568/2‑100 (MQ=21)
                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  2:3855021/2‑100 (MQ=18)
                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                                                                >  1:3519790/2‑98 (MQ=21)
                                gCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                >  2:4136408/1‑99 (MQ=18)
                                gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                                >  2:1423190/1‑98 (MQ=18)
                                gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                  >  2:3690750/1‑98 (MQ=17)
                                 cctgataaggcgTTTACGCCTTATCAGGCCTACAAGATCGTGCAAATTCAACATATTGCCACTC‑‑‑ACCCAGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                >  2:2651227/9‑99 (MQ=21)
                                 ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                               >  1:1072803/9‑98 (MQ=17)
                                 ccggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCATAACGCATCAGGCa                                                                 >  1:2668686/9‑98 (MQ=17)
                                 ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                               >  2:914270/9‑98 (MQ=21)
                                 ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                                                               >  1:4028903/9‑98 (MQ=17)
                                 ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                               >  1:4144737/1‑100 (MQ=255)
                                 ccGGATAAGGCGGTCACGCCGCATCCGGCCTACGAAACCGCAAGAATTCAATATATTGCAGGAGC‑‑‑GGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                               >  2:3085/1‑98 (MQ=18)
                                  cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                              <  2:3259963/100‑1 (MQ=255)
                                   ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                            >  1:2464264/7‑98 (MQ=21)
                                   ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  2:1209954/7‑100 (MQ=18)
                                   ggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                                                              >  1:3842269/7‑97 (MQ=21)
                                   ggataaggcgTTCACGCCTTATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  1:621797/7‑100 (MQ=21)
                                   ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  1:3784957/7‑100 (MQ=18)
                                   ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  1:2400924/7‑100 (MQ=18)
                                   ggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  1:2149184/7‑100 (MQ=21)
                                   ggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  2:1358885/7‑100 (MQ=18)
                                   ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  1:1468528/7‑100 (MQ=17)
                                   ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                             >  1:1899437/7‑100 (MQ=17)
                                    gataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                                                           >  1:203771/6‑97 (MQ=21)
                                    gataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                            <  2:2100866/95‑1 (MQ=17)
                                    gataaggcgTTCACGCCGCATCCTGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGAT‑‑‑GGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttttt                                                            >  1:3250564/6‑100 (MQ=14)
                                     aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                           <  2:396156/100‑1 (MQ=255)
                                            gTTAACGCCTTATCCGGCCTACAAAATCGTGCATATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                    >  2:2517834/1‑100 (MQ=255)
                                            gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                    <  1:1497072/100‑1 (MQ=255)
                                             ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATCCATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                   <  2:2698226/100‑1 (MQ=255)
                                              tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                  >  1:1182755/1‑100 (MQ=255)
                                               aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                 >  2:3610401/1‑100 (MQ=255)
                                                aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                <  2:1332259/100‑1 (MQ=255)
                                                aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                <  2:1706227/100‑1 (MQ=255)
                                                       aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                         <  2:3005166/100‑1 (MQ=255)
                                                        tCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                        >  1:710336/1‑100 (MQ=255)
                                                         ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                       <  1:1698568/100‑1 (MQ=255)
                                                           ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                     <  2:2185394/100‑1 (MQ=255)
                                                           ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                     <  2:257182/100‑1 (MQ=255)
                                                            gCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                    >  1:3936350/1‑100 (MQ=255)
                                                                aCAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                <  1:2935653/100‑1 (MQ=255)
                                                                aCAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                <  1:3149966/100‑1 (MQ=255)
                                                                   aaaaTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                              <  2:4046120/100‑1 (MQ=255)
                                                                        cGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                         <  1:3955061/100‑1 (MQ=255)
                                                                         gTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                        <  2:549819/100‑1 (MQ=255)
                                                                         gTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                        >  1:496008/1‑100 (MQ=255)
                                                                         gTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                        >  2:1069108/1‑100 (MQ=255)
                                                                          tGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                       <  1:940906/100‑1 (MQ=255)
                                                                          tGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                       >  2:3494064/1‑100 (MQ=255)
                                                                          tGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                       <  1:2629258/100‑1 (MQ=255)
                                                                          tGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                       <  1:1606358/100‑1 (MQ=255)
                                                                           gCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                      <  2:3975158/100‑1 (MQ=255)
                                                                           ccgatctCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                      <  1:4070815/94‑1 (MQ=255)
                                                                            cAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATGTTTGCCTTTGCAATCAGTCTTGTATGGCTTTTAattatt                        >  2:2951599/1‑97 (MQ=255)
                                                                              aaTTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                   <  1:660167/100‑1 (MQ=255)
                                                                                ttCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                 >  2:829143/1‑100 (MQ=255)
                                                                                ttCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                 <  1:2967100/100‑1 (MQ=255)
                                                                                   aaTACATTGCAGAATT‑‑‑TCC‑‑GTAGGCCTGATAAGCGGAGCGCATAAGGCATTTTTGACTTTGCCATCAGTCACGTATGGCTTTTACTTAGTCGTTAGAGAc              >  2:1139590/1‑100 (MQ=25)
                                                                                    aTACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg             >  2:215955/1‑100 (MQ=255)
                                                                                         ttGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc        <  1:3683144/100‑1 (MQ=255)
                                                                                         ttGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc        >  1:3895367/1‑100 (MQ=255)
                                                                                         ttGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc        >  2:223491/1‑100 (MQ=255)
                                                                                          tGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt         <  2:3441592/98‑1 (MQ=255)
                                                                                           gCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa      >  1:3726310/1‑100 (MQ=255)
                                                                                           gCAGAATG‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTGGCCATCAGTCTTGTATGGCGTTTAATTATTCGATAGAGACGTTTTCGa      <  2:2771372/100‑1 (MQ=255)
                                                                                           gCAGAATC‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa      <  1:1580369/100‑1 (MQ=255)
                                                                                               gAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt   <  2:4184273/100‑1 (MQ=255)
                                                                                               gAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt   >  1:3326233/1‑100 (MQ=255)
                                                                                                aaTT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg  >  1:2241978/1‑100 (MQ=255)
                                                                                                 aTT‑‑‑TTC‑‑GTAGGCCTGATAAGCCGAGCACATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg  <  2:505188/99‑1 (MQ=255)
                                                                                                    |                                                                                                      
TACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTG  >  NC_000913/3705931‑3706124

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: