Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 27.0% intergenic (‑258/+50) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,706,0480GT27.0% 613.8 / 217.5 295intergenic (‑258/+50)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (103/107);  new base T (74/2);  total (186/109)
Fisher's exact test for biased strand distribution p-value = 1.64e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 305 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147
                                                                                                        |                                                                                                   
    gattttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAgcg                                                                                                 >  2:1553592/11‑100 (MQ=18)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 <  1:1531116/100‑1 (MQ=255)
      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                                               <  2:3759090/100‑1 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                             >  2:2683847/1‑100 (MQ=255)
         tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  1:639892/6‑100 (MQ=17)
          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  2:2557763/5‑100 (MQ=17)
           taggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  2:289032/4‑100 (MQ=18)
           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATAAATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  1:3001215/4‑100 (MQ=17)
            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:2839328/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  2:396916/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                                         >  1:955708/3‑100 (MQ=18)
            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                         <  1:1191859/100‑1 (MQ=255)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                        >  2:1910432/1‑100 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTTCAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑ATAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:2569422/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGct                                                                                        >  2:598624/2‑99 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  2:1153356/2‑100 (MQ=25)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                        >  1:3434330/1‑99 (MQ=255)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                      >  1:1646114/1‑100 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  1:438969/1‑99 (MQ=17)
               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                      >  2:4237665/9‑98 (MQ=21)
               ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGTGCATCAGGCAtt                                                                                      <  2:2811865/100‑1 (MQ=255)
                cggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                                                     >  1:1786424/8‑100 (MQ=17)
                cggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatc                                                                                     >  2:1626921/8‑97 (MQ=25)
                cGGATAAGGCGGTCACGCCGCATCCGGCCTACGAAACCGCAAGAATTCAATATCTTGCAGGAGCGGT‑‑GTAGGCCTGATAAGCGTAGCGAATCAGGCa                                                                                        >  1:1565797/1‑97 (MQ=18)
                 tgataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                                    <  1:3466821/94‑5 (MQ=21)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:785283/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                                                                                    >  1:280916/7‑96 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattc                                                                                    >  2:1327207/7‑99 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                    >  1:101432/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                    >  1:1320148/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                                                    >  1:224714/7‑96 (MQ=25)
                  gataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctta                                                                                   >  2:3466834/6‑95 (MQ=21)
                  gataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaactac                                                                                   >  2:3734328/6‑95 (MQ=25)
                  gataacgcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCCGAAACCCT‑‑GTAGGCCTCAACAGCGGCGCGGATCACGc                                                                                         >  1:2085452/7‑94 (MQ=9)
                  gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                                                   <  2:2747352/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  2:4112468/100‑1 (MQ=255)
                    taaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAAGCAGTtttgc                                                                                 >  2:1971434/4‑100 (MQ=17)
                    tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                 >  2:2425261/1‑100 (MQ=255)
                    tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                 >  1:1424529/1‑100 (MQ=255)
                     gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                                <  2:224487/100‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATACGCGGAGCGCATCAGGCATTtttgcct                                                                               <  2:2605978/100‑1 (MQ=255)
                       ggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgcgtt                                                                              <  1:1948023/100‑1 (MQ=21)
                         cgTTCACGCCTTATCCGGCCTACAAAATCATGCAAACTCAATACATTGCCGAATTTTC‑‑GTAGGCCTGATAAGCGGAACGCATCAGGCATTtttgcct                                                                               >  1:228780/1‑97 (MQ=17)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTttttccttt                                                                             >  2:429108/1‑99 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                           <  2:871169/99‑1 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                         >  2:348360/1‑100 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATGGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                         <  1:2116188/100‑1 (MQ=255)
                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       >  2:2132203/1‑100 (MQ=255)
                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       >  1:3835576/1‑100 (MQ=255)
                                     aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  1:3765750/100‑1 (MQ=255)
                                     aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  1:1683638/100‑1 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  1:692668/100‑1 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  1:1910428/100‑1 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              >  2:4073396/1‑100 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              >  1:3990390/1‑100 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  2:2507916/100‑1 (MQ=255)
                                          gCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                                           <  2:1841856/100‑1 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                          <  2:3862681/100‑1 (MQ=255)
                                            cTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                                         <  1:3177516/100‑1 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       >  1:4147880/1‑100 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       >  1:3591117/1‑100 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       >  2:532106/1‑100 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCAGAATTTTA‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAGTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       >  2:3016187/1‑100 (MQ=255)
                                               cAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCGGGCATTTTTGCCTTTGCCATCAGTATTGTATGGCtttt                                                      <  1:2630123/100‑1 (MQ=255)
                                                   aTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTATGCCTGATAAGCGGAGCGCACCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                   <  1:1465236/99‑1 (MQ=255)
                                                      gTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                               >  2:708234/1‑100 (MQ=255)
                                                       tGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  2:2028906/100‑1 (MQ=255)
                                                       tGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              >  2:1511471/1‑100 (MQ=255)
                                                       tGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:3497755/100‑1 (MQ=255)
                                                        gCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                             <  1:2911919/100‑1 (MQ=255)
                                                          aaaTTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCTGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:2739882/97‑1 (MQ=255)
                                                              tCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                       <  2:3150924/100‑1 (MQ=255)
                                                               cAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                      <  2:3896522/100‑1 (MQ=255)
                                                                aaTACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCCATagag                                       <  1:3725661/98‑1 (MQ=255)
                                                                      ttGCAGAATTTTC‑‑GTAGGCCTGATAATCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  2:2597930/1‑100 (MQ=255)
                                                                       tGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                              >  2:2909039/1‑100 (MQ=255)
                                                                        gCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                              >  1:2529320/1‑99 (MQ=255)
                                                                        gCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             <  1:2511206/100‑1 (MQ=255)
                                                                        gCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  2:703505/1‑100 (MQ=255)
                                                                         cAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                            >  1:1594469/1‑100 (MQ=255)
                                                                          aGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  2:3245726/100‑1 (MQ=255)
                                                                            gAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          <  2:2062951/100‑1 (MQ=255)
                                                                             aaTTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                         <  1:814430/100‑1 (MQ=255)
                                                                               aTTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGa                        <  2:1268009/100‑1 (MQ=255)
                                                                                  tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGTCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                       >  1:2555370/1‑99 (MQ=255)
                                                                                  tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATg                      <  2:3331996/100‑1 (MQ=255)
                                                                                   ttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGt                     <  1:2423206/100‑1 (MQ=255)
                                                                                     c‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGttt                   >  2:1471837/1‑100 (MQ=255)
                                                                                     c‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGttt                   >  2:3634255/1‑100 (MQ=255)
                                                                                            gCCTGATAAGCGGAGCGCATCAGGCATTTTTGTCTTTGCAATCAGTATTGTATGGCTTTTACTTATTCGATAGAGACGTTATAGAAGTGATGTTTGCCTa              <  1:4066723/100‑1 (MQ=37)
                                                                                             ccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa             >  2:1670034/1‑100 (MQ=255)
                                                                                             ccTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa             <  2:2405385/100‑1 (MQ=255)
                                                                                              cTGATAAGCGGAGCGCATAAGGCATTTTTGCCTTTTCCATCAGTCTTGTATCGCTTTTAATTATTCGATAGAGCCGTTTACGAAGTGATGTTTGCCTAAt            >  1:3948043/1‑100 (MQ=255)
                                                                                               tGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAGTCTTGTATGGCCTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTCAt            >  1:3879287/1‑99 (MQ=255)
                                                                                               tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           >  2:1764514/1‑100 (MQ=255)
                                                                                                   aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc       <  1:4083380/100‑1 (MQ=255)
                                                                                                   aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc       <  1:2039935/100‑1 (MQ=255)
                                                                                                    aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa      >  1:2189407/1‑100 (MQ=255)
                                                                                                      cGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAc    <  2:3409514/100‑1 (MQ=255)
                                                                                                       ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc   <  1:3083738/100‑1 (MQ=255)
                                                                                                        |                                                                                                   
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

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