Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 37.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,0480GT37.0% 489.8 / inf 328intergenic (‑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 (90/111);  new base T (113/3);  total (214/114)
Fisher's exact test for biased strand distribution p-value = 4.31e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

Only 100 of 338 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147
                                                                                                            |                                                                                                   
tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                                                                                     <  1:2174010/100‑1 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  >  1:3912159/1‑100 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                   <  1:1042068/99‑1 (MQ=255)
       tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:2618801/8‑100 (MQ=17)
       tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:2389468/8‑100 (MQ=17)
       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                              >  2:1606354/1‑100 (MQ=255)
       tGTCTATGCCAGAGGGGGCCTTAACGGCGTATCCGGCCTACAAGAGCGTGCAAATTCAATTCGTTGCAGAAGT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGTgc                                                                                                  <  2:202789/96‑1 (MQ=12)
        ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  1:3675378/7‑100 (MQ=18)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                             >  1:2118781/1‑100 (MQ=255)
         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  1:1495145/6‑100 (MQ=18)
         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  2:2963406/6‑100 (MQ=17)
          gaaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTT‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAg                                                                                           >  1:2830811/5‑100 (MQ=18)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           <  2:1484001/99‑1 (MQ=255)
            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTT‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc                                                                                         >  2:2507734/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:1359021/3‑100 (MQ=21)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:2197207/3‑100 (MQ=17)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:366725/3‑100 (MQ=21)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                        <  2:3006328/100‑1 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:3237597/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                                                                                        >  1:3112745/2‑99 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATGCATTGCACGGGC‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:1137533/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:75441/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                         >  1:3585531/2‑100 (MQ=21)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTC‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:1888515/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:3855021/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                                                                                       >  1:3519790/2‑98 (MQ=21)
             gcccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:3893438/11‑100 (MQ=17)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAggaa                                                                                        >  2:2838325/1‑97 (MQ=255)
              gCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  1:1816695/1‑99 (MQ=17)
              gCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  2:4136408/1‑99 (MQ=18)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                      >  1:3365744/1‑100 (MQ=18)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                                                       >  2:1423190/1‑98 (MQ=18)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  1:2267381/1‑100 (MQ=21)
               cctgataaggcgTTTACGCCTTATCAGGCCTACAAGATCGTGCAAATTCAACATATTGCCACTC‑‑‑ACCCAGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                       >  2:2651227/9‑99 (MQ=21)
               ccggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                      >  1:2446499/9‑98 (MQ=18)
               ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                      >  1:1072803/9‑98 (MQ=17)
               ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                       >  1:92542/9‑99 (MQ=21)
               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                      >  2:914270/9‑98 (MQ=21)
               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                                                                                      >  1:4028903/9‑98 (MQ=17)
               ccggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                                                     >  1:2554907/9‑100 (MQ=21)
               ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                                      >  2:1190945/1‑100 (MQ=255)
                cggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                                                     >  2:2109531/8‑97 (MQ=17)
                cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                                                     <  2:3259963/100‑1 (MQ=255)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:3278712/7‑100 (MQ=17)
                 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:860604/7‑100 (MQ=18)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  1:116989/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  2:396156/100‑1 (MQ=255)
                        gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                                                             >  1:1514490/1‑100 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  1:2731698/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  2:3091136/100‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATCCATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                          <  2:2698226/100‑1 (MQ=255)
                             aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                                        >  2:3610401/1‑100 (MQ=255)
                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       <  1:448589/100‑1 (MQ=255)
                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       <  2:1332259/100‑1 (MQ=255)
                                  cTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                                                   >  1:131155/1‑100 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCaagatcgga                                                              >  1:2638605/1‑92 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  1:672771/100‑1 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  1:1698568/100‑1 (MQ=255)
                                         ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                            <  2:257182/100‑1 (MQ=255)
                                            cTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                                         >  2:3146875/1‑100 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       <  2:3782099/100‑1 (MQ=255)
                                                 aaaaTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                     <  2:4046120/100‑1 (MQ=255)
                                                  aaaTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                    <  1:4120659/100‑1 (MQ=255)
                                                   aaTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                   <  1:2975147/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                                               >  2:1069108/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              >  2:3494064/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:1606358/100‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:732078/100‑1 (MQ=255)
                                                         gCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                             <  2:3975158/100‑1 (MQ=255)
                                                            aaTTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                                          <  1:1903744/100‑1 (MQ=255)
                                                              ttCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                                        >  2:829143/1‑100 (MQ=255)
                                                                       ttGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTGTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  1:2940773/1‑100 (MQ=255)
                                                                       ttGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  2:223491/1‑100 (MQ=255)
                                                                        tGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt                                <  2:3441592/98‑1 (MQ=255)
                                                                         gCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGTATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             <  2:1772733/100‑1 (MQ=255)
                                                                         gCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             <  2:1719885/100‑1 (MQ=255)
                                                                         gCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  1:3726310/1‑100 (MQ=255)
                                                                         gCAGAATG‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTGGCCATCAGTCTTGTATGGCGTTTAATTATTCGATAGAGACGTTTTCGa                             <  2:2771372/100‑1 (MQ=255)
                                                                           aGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  1:506221/100‑1 (MQ=255)
                                                                             gAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          <  2:4184273/100‑1 (MQ=255)
                                                                                aTT‑‑‑TTC‑‑GTAGGCCTGATAAGCCGAGCACATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                         <  2:505188/99‑1 (MQ=255)
                                                                                             tAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGc                 <  1:3610598/100‑1 (MQ=255)
                                                                                              aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc                <  2:1681876/100‑1 (MQ=255)
                                                                                               ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               <  1:4058084/100‑1 (MQ=255)
                                                                                                gCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTa              <  2:2452789/100‑1 (MQ=255)
                                                                                                   tGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           >  2:1278312/1‑100 (MQ=255)
                                                                                                   tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  1:490018/100‑1 (MQ=255)
                                                                                                    gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCCGTCTTGTATGGCTTTTAATTATTCCATAGAGACGTTTTCGAAGTGATGTTTGCGTAATgg          >  2:1359037/1‑100 (MQ=255)
                                                                                                     aTAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGa         <  1:3417110/100‑1 (MQ=255)
                                                                                                       aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc       <  2:664561/100‑1 (MQ=255)
                                                                                                        aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt        >  2:3758888/1‑98 (MQ=255)
                                                                                                        aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa      <  2:2165512/100‑1 (MQ=255)
                                                                                                         gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa     >  1:1206760/1‑100 (MQ=255)
                                                                                                         gCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCaa     <  2:3369003/100‑1 (MQ=255)
                                                                                                            gAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACca  <  1:1937218/100‑1 (MQ=255)
                                                                                                            |                                                                                                   
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑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: