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.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116
                                                                                                   |                                                                                                        
gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                                                          <  2:45238/100‑1 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCTGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTTCAAATTCAATACATTGCAGAAtt‑‑tt                                                                                          >  2:117613/1‑100 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATCACGGATGCAGCGTTAACTCCTTATCCGGCCTACAAAATCGTGCAAATTCATTACATTGCAGAAtt‑‑tt                                                                                          <  2:85273/100‑1 (MQ=17)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑g                                                                                      >  1:2857584/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑g                                                                                      >  1:1702126/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑g                                                                                      <  2:2293156/100‑1 (MQ=255)
                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑Gt                                                                                     >  1:701618/1‑100 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCATATTCAATACATTGCAGAATT‑‑TTC‑‑GTa                                                                                    <  1:1843893/100‑1 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTa                                                                                    <  2:371916/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAgg                                                                                  <  1:1201207/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAgg                                                                                  <  2:546365/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGTTTTATCCTGTTTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGcc                                                                                <  1:983872/100‑1 (MQ=18)
                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTATGCCt                                                                               >  2:1931382/1‑100 (MQ=255)
                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGAt                                                                            <  2:2811126/100‑1 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGc                                                                        >  1:2827266/1‑100 (MQ=255)
                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCgg                                                                      <  1:330353/100‑1 (MQ=255)
                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATCAGCGGa                                                                     >  2:1217667/1‑100 (MQ=255)
                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                    >  1:587145/1‑100 (MQ=255)
                                 cgcTTTGTTGATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                    <  2:2570162/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  <  1:3040692/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                   <  1:405163/99‑1 (MQ=255)
                                      tttgtaggCCGGATAAGTCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTACAGAGAT‑‑TGC‑‑GTAGGCCTGATAAGCGTAg                                                                    >  2:1462092/8‑95 (MQ=11)
                                      tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                               >  1:2836239/8‑100 (MQ=18)
                                      tGTTTCTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAAT‑‑TTC‑‑GTACGCCTGATAAGCGGAgcgc                                                                 >  2:1938451/1‑98 (MQ=255)
                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAATGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                               >  2:1671444/1‑100 (MQ=255)
                                       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATc                                                              >  2:2399942/7‑100 (MQ=18)
                                        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                             >  1:1798739/6‑100 (MQ=18)
                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                            >  2:3020255/5‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑CCT‑‑GTAGGCCGGATAAGCGTAGCGCATCAGGCa                                                         >  2:2050239/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:85464/2‑100 (MQ=14)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                         >  1:509184/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:1106044/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGAGAC‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                                                         >  1:401355/2‑100 (MQ=17)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                        >  1:1233024/1‑100 (MQ=255)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                        >  1:370551/1‑100 (MQ=255)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACACAAACATGCAAATTCAATAGATTGCAGAGAT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                        >  1:25540/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat                                                       >  1:494082/1‑97 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                        >  1:2262768/1‑99 (MQ=18)
                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                        >  1:449408/9‑99 (MQ=21)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                       >  1:146806/9‑98 (MQ=17)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                       >  1:2028258/1‑100 (MQ=255)
                                              ccGGATCCGGCGTTAACGCCTTATCCAGCCTACAAAATTGTGCAAAATCAATTGATTGCACAGCT‑‑AGC‑‑GTAGGCCTGATGCGCTACGCTCATCAGGc                                                          <  2:1483766/97‑1 (MQ=18)
                                               cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                       >  2:1357700/8‑97 (MQ=18)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1008095/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCGTCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGGTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                     >  1:2580452/7‑96 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1559778/7‑100 (MQ=14)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:131475/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                     >  2:756655/7‑96 (MQ=25)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                     >  2:1708710/7‑96 (MQ=25)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                      >  1:3034549/7‑99 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACT‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                     >  2:1239319/7‑96 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                    >  1:1066030/7‑98 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACG‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                     >  2:2269240/7‑96 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCATACAAATCATGCAAATCCAAAATATTGCAGGAGA‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACAttta                                                     >  1:2756742/7‑98 (MQ=18)
                                                 gataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                    >  2:2192871/6‑100 (MQ=17)
                                                  aTGGGGCGTTATCGCTTTATAGGGCGAAAAAAATCGTGGAAATTCACTACATTGCAGAATT‑‑TTC‑‑GGAGGGCTAATAAGTGGAGCGCATCAGGCATTtttg                                                   <  2:1974961/100‑1 (MQ=16)
                                                  aTGCGGGGTTAACTCCTTATCCGGCATACAAAATGGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTTATAAGCGGAGTGCATCAGGGATTtttg                                                   <  2:1270139/100‑1 (MQ=18)
                                                    gcggcgTTAACGCCTTATCCAGCCTACAAAATAGTACAAAATCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGAGCATCAGGCATTtttgcc                                                 >  2:704377/1‑100 (MQ=18)
                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                <  1:2913962/100‑1 (MQ=255)
                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                >  1:2748791/1‑100 (MQ=255)
                                                      ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                               >  1:1614734/1‑100 (MQ=255)
                                                      ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCATGATAAGCGGAGCGCATCAGGCATTtttgcctt                                               >  1:6083/1‑100 (MQ=255)
                                                        cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             >  2:636789/1‑100 (MQ=255)
                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGTAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             <  2:2555109/98‑1 (MQ=255)
                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                           <  2:475307/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                           >  2:1710077/1‑99 (MQ=255)
                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          <  2:1381428/100‑1 (MQ=255)
                                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                         >  1:2196911/1‑100 (MQ=255)
                                                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        <  1:2420009/100‑1 (MQ=255)
                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCa                                       >  1:2333888/1‑100 (MQ=255)
                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGATTT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                      <  1:773849/100‑1 (MQ=255)
                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                      <  1:1116001/100‑1 (MQ=255)
                                                                ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTCCTTTGCCATCatt                                     >  2:3045457/1‑98 (MQ=255)
                                                                ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                     >  2:1744233/1‑100 (MQ=255)
                                                                       cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              >  1:464485/1‑100 (MQ=255)
                                                                         gCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑TTAGGCCTGATAAGCGTAGCGCATCATGCATTTTTGCCTTTGCCATCAGTCtt                                  <  2:2656694/94‑1 (MQ=255)
                                                                         gCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              >  1:1982441/1‑98 (MQ=255)
                                                                          ccTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                           <  2:1831134/100‑1 (MQ=255)
                                                                           tctcAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                          <  2:2102119/97‑1 (MQ=255)
                                                                              cAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtaga                       >  1:2102142/1‑97 (MQ=255)
                                                                                 aaTCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                    <  2:496477/100‑1 (MQ=255)
                                                                                  aTCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                   >  2:919909/1‑100 (MQ=255)
                                                                                   tCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                  <  2:796912/100‑1 (MQ=255)
                                                                                   tCGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                  <  2:388500/100‑1 (MQ=255)
                                                                                    cGTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                 >  1:1187900/1‑100 (MQ=255)
                                                                                     gTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                >  2:2680400/1‑100 (MQ=255)
                                                                                     gTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                >  1:1933076/1‑100 (MQ=255)
                                                                                     gTGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                >  2:2082608/1‑100 (MQ=255)
                                                                                      tGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCTCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  1:997315/100‑1 (MQ=255)
                                                                                      tGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               >  1:2595544/1‑100 (MQ=255)
                                                                                      tGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  1:1549065/100‑1 (MQ=255)
                                                                                      tGCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  2:2148835/100‑1 (MQ=255)
                                                                                       gCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              >  2:2894226/1‑100 (MQ=255)
                                                                                       gCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              <  1:1476580/100‑1 (MQ=255)
                                                                                       gCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              >  2:1334593/1‑100 (MQ=255)
                                                                                       gCAAATTCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              >  2:2394884/1‑100 (MQ=255)
                                                                                         aaaTTCAATACATTGCAGAAGT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAGCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt            <  2:2741684/100‑1 (MQ=255)
                                                                                             tCAATACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag        >  1:962990/1‑100 (MQ=255)
                                                                                                aTACATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg     >  2:1666298/1‑100 (MQ=255)
                                                                                                   cATTGCAGAATT‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt  >  1:68886/1‑100 (MQ=255)
                                                                                                   |                                                                                                        
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑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: