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.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147
                                                                                                        |                                                                                                   
    gattttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAgcg                                                                                                 >  2:1553592/11‑100 (MQ=18)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 <  1:1531116/100‑1 (MQ=255)
      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                                               <  2:3759090/100‑1 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                             >  2:2683847/1‑100 (MQ=255)
         tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCA‑GA‑GA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  1:639892/6‑100 (MQ=17)
          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  2:2557763/5‑100 (MQ=17)
           taggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  2:289032/4‑100 (MQ=18)
           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATAAATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                          >  1:3001215/4‑100 (MQ=17)
            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑‑A‑AT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:2839328/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  2:396916/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                                         >  1:955708/3‑100 (MQ=18)
            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                                         <  1:1191859/100‑1 (MQ=255)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                        >  2:1910432/1‑100 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTTCAAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑ATAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:2569422/2‑100 (MQ=17)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑‑A‑AT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGct                                                                                        >  2:598624/2‑99 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  2:1153356/2‑100 (MQ=25)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                        >  1:3434330/1‑99 (MQ=255)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                      >  1:1646114/1‑100 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                       >  1:438969/1‑99 (MQ=17)
               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                      >  2:4237665/9‑98 (MQ=21)
               ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGTGCATCAGGCAtt                                                                                      <  2:2811865/100‑1 (MQ=255)
                cggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                                                     >  1:1786424/8‑100 (MQ=17)
                cggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatc                                                                                     >  2:1626921/8‑97 (MQ=25)
                cGGATAAGGCGGTCACGCCGCATCCGGCCTACGAAACCGCAAGAATTCAATATCTTGCA‑GG‑AG‑CGGT‑‑GTAGGCCTGATAAGCGTAGCGAATCAGGCa                                                                                        >  1:1565797/1‑97 (MQ=18)
                 tgataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                                                    <  1:3466821/94‑5 (MQ=21)
                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:785283/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                                                                                    >  1:280916/7‑96 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattc                                                                                    >  2:1327207/7‑99 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                    >  1:101432/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCA‑GA‑GA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                                                    >  1:1320148/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGA‑AC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                                                    >  1:224714/7‑96 (MQ=25)
                  gataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctta                                                                                   >  2:3466834/6‑95 (MQ=21)
                  gataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGA‑AC‑GAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaactac                                                                                   >  2:3734328/6‑95 (MQ=25)
                  gataacgcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCC‑GA‑AA‑CCCT‑‑GTAGGCCTCAACAGCGGCGCGGATCACGc                                                                                         >  1:2085452/7‑94 (MQ=9)
                  gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                                                   <  2:2747352/100‑1 (MQ=255)
                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  2:4112468/100‑1 (MQ=255)
                    taaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAAGCAGTtttgc                                                                                 >  2:1971434/4‑100 (MQ=17)
                    tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                 >  2:2425261/1‑100 (MQ=255)
                    tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                                                 >  1:1424529/1‑100 (MQ=255)
                     gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                                <  2:224487/100‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATACGCGGAGCGCATCAGGCATTtttgcct                                                                               <  2:2605978/100‑1 (MQ=255)
                       ggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgcgtt                                                                              <  1:1948023/100‑1 (MQ=21)
                         cgTTCACGCCTTATCCGGCCTACAAAATCATGCAAACTCAATACATTGCC‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAACGCATCAGGCATTtttgcct                                                                               >  1:228780/1‑97 (MQ=17)
                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTttttccttt                                                                             >  2:429108/1‑99 (MQ=255)
                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                           <  2:871169/99‑1 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                         >  2:348360/1‑100 (MQ=255)
                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATGGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                         <  1:2116188/100‑1 (MQ=255)
                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       >  2:2132203/1‑100 (MQ=255)
                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                       >  1:3835576/1‑100 (MQ=255)
                                     aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  1:3765750/100‑1 (MQ=255)
                                     aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  1:1683638/100‑1 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  1:692668/100‑1 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  1:1910428/100‑1 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              >  2:4073396/1‑100 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              >  1:3990390/1‑100 (MQ=255)
                                       ccGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                              <  2:2507916/100‑1 (MQ=255)
                                          gCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                                           <  2:1841856/100‑1 (MQ=255)
                                           ccTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                          <  2:3862681/100‑1 (MQ=255)
                                            cTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                                                         <  1:3177516/100‑1 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       >  1:4147880/1‑100 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       >  1:3591117/1‑100 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       >  2:532106/1‑100 (MQ=255)
                                              aCAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTA‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAGTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                       >  2:3016187/1‑100 (MQ=255)
                                               cAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCGGGCATTTTTGCCTTTGCCATCAGTATTGTATGGCtttt                                                      <  1:2630123/100‑1 (MQ=255)
                                                   aTCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTATGCCTGATAAGCGGAGCGCACCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                   <  1:1465236/99‑1 (MQ=255)
                                                      gTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                               >  2:708234/1‑100 (MQ=255)
                                                       tGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  2:2028906/100‑1 (MQ=255)
                                                       tGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              >  2:1511471/1‑100 (MQ=255)
                                                       tGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:3497755/100‑1 (MQ=255)
                                                        gCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                             <  1:2911919/100‑1 (MQ=255)
                                                          aaaTTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCTGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              <  1:2739882/97‑1 (MQ=255)
                                                              tCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag                                       <  2:3150924/100‑1 (MQ=255)
                                                               cAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga                                      <  2:3896522/100‑1 (MQ=255)
                                                                aaTACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCCATagag                                       <  1:3725661/98‑1 (MQ=255)
                                                                      ttGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAATCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  2:2597930/1‑100 (MQ=255)
                                                                       tGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                              >  2:2909039/1‑100 (MQ=255)
                                                                        gCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                              >  1:2529320/1‑99 (MQ=255)
                                                                        gCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             <  1:2511206/100‑1 (MQ=255)
                                                                        gCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                             >  2:703505/1‑100 (MQ=255)
                                                                         cA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                            >  1:1594469/1‑100 (MQ=255)
                                                                          a‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  2:3245726/100‑1 (MQ=255)
                                                                            gA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          <  2:2062951/100‑1 (MQ=255)
                                                                             a‑aT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                         <  1:814430/100‑1 (MQ=255)
                                                                               aT‑TTTC‑‑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)
                                                                                                        |                                                                                                   
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑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: ‑