Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 32.5% intergenic (‑195/+113) 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,705,9850TC32.5% 183.2 / 120.0 130intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (39/48);  new base C (39/3);  total (78/52)
Fisher's exact test for biased strand distribution p-value = 3.98e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-03

Only 100 of 133 total aligned reads displayed.

GCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705890‑3706084
                                                                                               |                                                                                                           
gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                         >  2:451302/1‑100 (MQ=255)
  cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                       <  2:1793185/100‑1 (MQ=255)
   ggtgggCTCTGCACGCAGTGGCACAGCGCGACGGCGTATACAGATTGTTGACAAAGTGCGCTTGGTTTATGCTGGATAAGGCGTTAACGCCCGATCCGGc                                                                                                      <  2:1485043/96‑1 (MQ=18)
   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                      <  1:1659655/100‑1 (MQ=255)
   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                      >  1:2031096/1‑100 (MQ=255)
   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                      >  1:533609/1‑100 (MQ=255)
   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                      <  1:2265040/100‑1 (MQ=37)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTTTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                     <  1:1680355/100‑1 (MQ=255)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                     <  2:1527873/100‑1 (MQ=255)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                     <  2:2246673/100‑1 (MQ=255)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                     <  2:257273/100‑1 (MQ=255)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                     <  1:1026546/100‑1 (MQ=255)
     aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                    <  2:1049687/100‑1 (MQ=255)
        gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                                >  2:1289056/1‑100 (MQ=255)
                 gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                       <  2:1625146/100‑1 (MQ=255)
                 gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                       <  1:1757996/100‑1 (MQ=255)
                  cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                      <  2:95124/100‑1 (MQ=255)
                   agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                       <  1:970509/98‑1 (MQ=255)
                    gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                                    >  2:367193/1‑100 (MQ=255)
                     agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                                   >  2:910870/1‑100 (MQ=255)
                       gCACAGCGCCACCGGGTATACAGATTGTTGACCAAGTGCGCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                                 <  2:1672339/100‑1 (MQ=255)
                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                                 <  1:531135/100‑1 (MQ=255)
                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                                 <  1:649986/100‑1 (MQ=255)
                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                                >  2:1583913/1‑100 (MQ=255)
                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                                <  2:2132931/100‑1 (MQ=255)
                         acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                               <  1:595832/100‑1 (MQ=255)
                         acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                               >  1:583378/1‑100 (MQ=255)
                            gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                            >  1:791972/1‑100 (MQ=255)
                               ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                         >  2:2019038/1‑100 (MQ=255)
                                   cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                       <  2:2011285/98‑1 (MQ=255)
                                   cGGGTATACAGATTGTTCACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                                      <  2:2168621/99‑1 (MQ=255)
                                     ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                                  >  2:1997370/1‑100 (MQ=255)
                                         taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                          >  2:1234753/1‑100 (MQ=255)
                                          aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTc                                                         >  2:336257/1‑100 (MQ=255)
                                           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑g                                                      <  1:1540601/100‑1 (MQ=255)
                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAgg                                                  >  1:1556835/7‑100 (MQ=18)
                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAgg                                                  >  1:224687/7‑100 (MQ=18)
                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAgg                                                  >  1:1159376/7‑100 (MQ=17)
                                                tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGc                                                 >  1:1266934/1‑100 (MQ=255)
                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTg                                              <  1:1234760/100‑1 (MQ=255)
                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTg                                              >  1:389865/1‑100 (MQ=255)
                                                       aaatggcgcTTTATATAGGCCGGATGCGGCGTTAACGCCCTATCAGGCCTGCAAAAGCGTGCATTTTCACTACATTGCATAATT‑‑‑TTC‑‑GTAGTCCTGATa                                           <  2:1548107/94‑1 (MQ=16)
                                                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTCCAAACTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                      <  2:1544585/100‑1 (MQ=255)
                                                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                      <  2:1472279/100‑1 (MQ=255)
                                                            gcgcTTTGTTTATGTCGGATGCGGCGTTAACGCCTTATCCGGGCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                     <  2:851541/100‑1 (MQ=255)
                                                               gattttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAgcg                                  >  1:1437372/11‑100 (MQ=18)
                                                               gaatttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAgcgc                                 >  1:1962390/11‑100 (MQ=18)
                                                                 ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCa                                >  1:1368646/9‑100 (MQ=18)
                                                                 ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                <  1:1472132/100‑1 (MQ=255)
                                                                   gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                              >  2:797622/1‑100 (MQ=255)
                                                                   gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                              <  2:1065081/100‑1 (MQ=255)
                                                                    tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                             <  1:1806524/100‑1 (MQ=255)
                                                                     gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  1:61560/5‑100 (MQ=21)
                                                                     gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:1017640/5‑100 (MQ=21)
                                                                      taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:1923143/4‑100 (MQ=17)
                                                                      taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAGGCGTAGCGCATCAgg                           >  1:2351496/4‑100 (MQ=18)
                                                                      taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:2073405/4‑100 (MQ=21)
                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                         >  1:1094191/3‑99 (MQ=21)
                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:1311126/3‑100 (MQ=18)
                                                                       aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                          >  1:2300622/3‑100 (MQ=18)
                                                                       aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                          >  1:1154672/1‑100 (MQ=255)
                                                                        ggcctgataaggcgTTTACGCCGCATCAGGCCTACAAGATCGTGCAAATTCAACATATTGCCACTC‑‑‑ACCCAGTAGGCCTGATAAGCGCAGCGCATCAGGc                          >  2:456011/11‑100 (MQ=21)
                                                                        ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:181213/2‑100 (MQ=17)
                                                                        ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                         >  2:2323033/2‑100 (MQ=18)
                                                                         gCCGGATGCGGCGTGAACGCCTTATCCGGCCTATAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATGCGCTGCGCTTATCAGGcct                        <  1:464619/100‑3 (MQ=21)
                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:1884715/1‑99 (MQ=18)
                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:2153720/1‑99 (MQ=18)
                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:2360505/1‑99 (MQ=18)
                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:1336598/1‑99 (MQ=18)
                                                                          ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAA‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCagt                       >  1:2178199/9‑98 (MQ=18)
                                                                          ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       <  1:802835/100‑1 (MQ=255)
                                                                          ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  1:150657/1‑100 (MQ=255)
                                                                          ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  1:2233425/1‑99 (MQ=255)
                                                                          caggataaggcgTTCACGCCGCATCCAGCCTACAAAATTGTGCAAAATCAATGGATTGCACAGCT‑‑‑AGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaga                       >  1:1611143/9‑98 (MQ=25)
                                                                           cggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatc                      >  1:1920055/8‑97 (MQ=25)
                                                                           cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                      >  1:2206527/8‑97 (MQ=21)
                                                                           cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                      >  1:1136997/8‑97 (MQ=21)
                                                                           cggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTC‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                      >  1:1713102/8‑97 (MQ=18)
                                                                           cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                      >  1:272095/1‑100 (MQ=255)
                                                                            ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:1233241/7‑100 (MQ=17)
                                                                            ggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCAatt                      >  1:309886/7‑100 (MQ=21)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:578233/7‑100 (MQ=18)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:834119/7‑100 (MQ=18)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:1329030/7‑100 (MQ=18)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:1329654/7‑100 (MQ=18)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                     >  1:96255/7‑96 (MQ=25)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                     >  2:70859/7‑96 (MQ=25)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:36153/7‑100 (MQ=17)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                     >  2:542949/7‑96 (MQ=21)
                                                                            ggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAAGCCTGATAAGCGGAGCGCATCAGGCAtttt                     >  1:992215/7‑100 (MQ=18)
                                                                             gataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                    >  2:1091493/6‑100 (MQ=17)
                                                                             gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                    <  1:1981795/100‑1 (MQ=255)
                                                                               tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                  >  1:1157620/1‑100 (MQ=255)
                                                                                gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                 >  1:879450/1‑100 (MQ=255)
                                                                                   gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt              >  1:833794/1‑100 (MQ=255)
                                                                                         aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  2:693720/100‑1 (MQ=255)
                                                                                            ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  1:1286161/1‑100 (MQ=255)
                                                                                              ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  2:56540/100‑1 (MQ=255)
                                                                                               tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  2:252956/1‑100 (MQ=255)
                                                                                               tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  1:417887/1‑100 (MQ=255)
                                                                                               |                                                                                                           
GCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705890‑3706084

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