Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,017 C→T 23.1% intergenic (‑227/+81) 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,0170CT23.1% 286.8 / 75.7 148intergenic (‑227/+81)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (61/47);  new base T (30/2);  total (98/50)
Fisher's exact test for biased strand distribution p-value = 4.28e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.96e-01

Only 100 of 151 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116
                                                                                                   |                                                                                                        
gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                                                          <  2:45238/100‑1 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCTGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTTCAAATTCAATACATTGCA‑GAAtt‑‑tt                                                                                          >  2:117613/1‑100 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATCACGGATGCAGCGTTAACTCCTTATCCGGCCTACAAAATCGTGCAAATTCATTACATTGCA‑GAAtt‑‑tt                                                                                          <  2:85273/100‑1 (MQ=17)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑g                                                                                      >  1:2857584/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑g                                                                                      >  1:1702126/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑g                                                                                      <  2:2293156/100‑1 (MQ=255)
                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑Gt                                                                                     >  1:701618/1‑100 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCATATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTa                                                                                    <  1:1843893/100‑1 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTa                                                                                    <  2:371916/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAgg                                                                                  <  1:1201207/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAgg                                                                                  <  2:546365/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGTTTTATCCTGTTTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGcc                                                                                <  1:983872/100‑1 (MQ=18)
                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTATGCCt                                                                               >  2:1931382/1‑100 (MQ=255)
                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGAt                                                                            <  2:2811126/100‑1 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGc                                                                        >  1:2827266/1‑100 (MQ=255)
                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCgg                                                                      <  1:330353/100‑1 (MQ=255)
                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATCAGCGGa                                                                     >  2:1217667/1‑100 (MQ=255)
                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                    >  1:587145/1‑100 (MQ=255)
                                 cgcTTTGTTGATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                    <  2:2570162/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  <  1:3040692/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                   <  1:405163/99‑1 (MQ=255)
                                      tttgtaggCCGGATAAGTCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTACA‑GAGAT‑‑TGC‑‑GTAGGCCTGATAAGCGTAg                                                                    >  2:1462092/8‑95 (MQ=11)
                                      tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCA‑ATTCT‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                               >  1:2836239/8‑100 (MQ=18)
                                      tGTTTCTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAAAT‑‑TTC‑‑GTACGCCTGATAAGCGGAgcgc                                                                 >  2:1938451/1‑98 (MQ=255)
                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAATGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                               >  2:1671444/1‑100 (MQ=255)
                                       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AATC‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATc                                                              >  2:2399942/7‑100 (MQ=18)
                                        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCA‑ATTCT‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                             >  1:1798739/6‑100 (MQ=18)
                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GGAGC‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                            >  2:3020255/5‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCA‑GAAAC‑‑CCT‑‑GTAGGCCGGATAAGCGTAGCGCATCAGGCa                                                         >  2:2050239/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GAGTT‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:85464/2‑100 (MQ=14)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                         >  1:509184/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGAT‑CGATC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:1106044/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCA‑GAGAC‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                                                         >  1:401355/2‑100 (MQ=17)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                        >  1:1233024/1‑100 (MQ=255)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                        >  1:370551/1‑100 (MQ=255)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACACAAACATGCAAATTCAATAGATTGCA‑GAGAT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                        >  1:25540/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCA‑‑ATAT‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat                                                       >  1:494082/1‑97 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AATC‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                        >  1:2262768/1‑99 (MQ=18)
                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCA‑‑ATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                        >  1:449408/9‑99 (MQ=21)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TGAGC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                       >  1:146806/9‑98 (MQ=17)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                       >  1:2028258/1‑100 (MQ=255)
                                              ccGGATCCGGCGTTAACGCCTTATCCAGCCTACAAAATTGTGCAAAATCAATTGATTGCA‑CAGCT‑‑AGC‑‑GTAGGCCTGATGCGCTACGCTCATCAGGc                                                          <  2:1483766/97‑1 (MQ=18)
                                               cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCA‑CCATC‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                       >  2:1357700/8‑97 (MQ=18)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GAGTT‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1008095/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCGTCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GGTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                     >  1:2580452/7‑96 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TGAGC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1559778/7‑100 (MQ=14)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCA‑TGAGC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:131475/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                     >  2:756655/7‑96 (MQ=25)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCA‑GTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                     >  2:1708710/7‑96 (MQ=25)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCA‑GGAGC‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                      >  1:3034549/7‑99 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACT‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                     >  2:1239319/7‑96 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAGT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                    >  1:1066030/7‑98 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACG‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                     >  2:2269240/7‑96 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCATAC‑AAATCATGCAAATCCAAAATATTGCA‑GGAGA‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACAttta                                                     >  1:2756742/7‑98 (MQ=18)
                                                 gataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCA‑GAGTT‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                    >  2:2192871/6‑100 (MQ=17)
                                                  aTGGGGCGTTATCGCTTTATAGGGCGAAAAAAATCGTGGAAATTCACTACATTGCA‑GAATT‑‑TTC‑‑GGAGGGCTAATAAGTGGAGCGCATCAGGCATTtttg                                                   <  2:1974961/100‑1 (MQ=16)
                                                  aTGCGGGGTTAACTCCTTATCCGGCATACAAAATGGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTTATAAGCGGAGTGCATCAGGGATTtttg                                                   <  2:1270139/100‑1 (MQ=18)
                                                    gcggcgTTAACGCCTTATCCAGCCTACAAAATAGTACAAAATCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGAGCATCAGGCATTtttgcc                                                 >  2:704377/1‑100 (MQ=18)
                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                <  1:2913962/100‑1 (MQ=255)
                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                >  1:2748791/1‑100 (MQ=255)
                                                      ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                               >  1:1614734/1‑100 (MQ=255)
                                                      ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCATGATAAGCGGAGCGCATCAGGCATTtttgcctt                                               >  1:6083/1‑100 (MQ=255)
                                                        cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             >  2:636789/1‑100 (MQ=255)
                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGTAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             <  2:2555109/98‑1 (MQ=255)
                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                           <  2:475307/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                           >  2:1710077/1‑99 (MQ=255)
                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          <  2:1381428/100‑1 (MQ=255)
                                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                         >  1:2196911/1‑100 (MQ=255)
                                                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        <  1:2420009/100‑1 (MQ=255)
                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCa                                       >  1:2333888/1‑100 (MQ=255)
                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GATTT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                      <  1:773849/100‑1 (MQ=255)
                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                      <  1:1116001/100‑1 (MQ=255)
                                                                ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTCCTTTGCCATCatt                                     >  2:3045457/1‑98 (MQ=255)
                                                                ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                     >  2:1744233/1‑100 (MQ=255)
                                                                       cGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              >  1:464485/1‑100 (MQ=255)
                                                                         gCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑TTAGGCCTGATAAGCGTAGCGCATCATGCATTTTTGCCTTTGCCATCAGTCtt                                  <  2:2656694/94‑1 (MQ=255)
                                                                         gCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              >  1:1982441/1‑98 (MQ=255)
                                                                          ccTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                           <  2:1831134/100‑1 (MQ=255)
                                                                           tctcAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                          <  2:2102119/97‑1 (MQ=255)
                                                                              cAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtaga                       >  1:2102142/1‑97 (MQ=255)
                                                                                 aaTCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                    <  2:496477/100‑1 (MQ=255)
                                                                                  aTCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                   >  2:919909/1‑100 (MQ=255)
                                                                                   tCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                  <  2:796912/100‑1 (MQ=255)
                                                                                   tCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                  <  2:388500/100‑1 (MQ=255)
                                                                                    cGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                 >  1:1187900/1‑100 (MQ=255)
                                                                                     gTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                >  2:2680400/1‑100 (MQ=255)
                                                                                     gTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                >  1:1933076/1‑100 (MQ=255)
                                                                                     gTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                >  2:2082608/1‑100 (MQ=255)
                                                                                      tGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCTCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  1:997315/100‑1 (MQ=255)
                                                                                      tGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               >  1:2595544/1‑100 (MQ=255)
                                                                                      tGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  1:1549065/100‑1 (MQ=255)
                                                                                      tGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  2:2148835/100‑1 (MQ=255)
                                                                                       gCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              >  2:2894226/1‑100 (MQ=255)
                                                                                       gCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              <  1:1476580/100‑1 (MQ=255)
                                                                                       gCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              >  2:1334593/1‑100 (MQ=255)
                                                                                       gCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              >  2:2394884/1‑100 (MQ=255)
                                                                                         aaaTTCAATACATTGCA‑GAAGT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAGCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt            <  2:2741684/100‑1 (MQ=255)
                                                                                             tCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag        >  1:962990/1‑100 (MQ=255)
                                                                                                aTACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg     >  2:1666298/1‑100 (MQ=255)
                                                                                                   cATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt  >  1:68886/1‑100 (MQ=255)
                                                                                                   |                                                                                                        
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCA‑GAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116

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