Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 28.8% intergenic (‑194/+114) 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,9840TG28.8% 506.2 / 248.0 257intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (92/91);  new base G (69/5);  total (161/96)
Fisher's exact test for biased strand distribution p-value = 5.13e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-02

Only 100 of 266 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                         
ttGCTGCCGGATGTCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                           <  1:602182/100‑1 (MQ=255)
ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                           >  1:1525855/1‑100 (MQ=255)
ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                           >  1:2651115/1‑100 (MQ=255)
 tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                          >  2:1664284/1‑100 (MQ=255)
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                         >  2:2561829/1‑100 (MQ=255)
   ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                        >  2:299896/1‑100 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                       <  1:2080456/100‑1 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                      <  2:2095945/100‑1 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                      <  2:2834720/100‑1 (MQ=255)
      ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                     >  2:3311636/1‑100 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                   <  2:677890/100‑1 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                  <  1:3599058/100‑1 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                  <  2:1114287/100‑1 (MQ=255)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                >  1:1392216/1‑100 (MQ=255)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                             <  2:2210724/100‑7 (MQ=25)
               gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                           >  2:1995120/1‑100 (MQ=255)
               cctgtaggaagagGCGCAGTGCCACCGGGTATACAGATTGCTGACAAAGTGAGCTCTGATTATGCCCGATGCGTCGTTGACGCCTTATCCGGCCTAc                                                                                               >  1:1775629/10‑97 (MQ=17)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                       <  1:3309763/100‑1 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                       >  1:1604397/1‑100 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                     >  1:2985076/1‑100 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                     >  2:2055967/1‑100 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                     >  1:315175/1‑100 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                     <  1:3120527/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCATg                                                                                     >  2:1064694/1‑100 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                   >  2:2373106/1‑100 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                   >  2:2339947/1‑100 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                   >  1:3561829/1‑100 (MQ=255)
                         gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                                 <  1:1664276/100‑1 (MQ=255)
                          agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                                <  2:1633978/100‑1 (MQ=255)
                           ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                               <  2:3226075/100‑1 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                             >  2:3069004/1‑100 (MQ=255)
                              acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                            >  2:1905814/1‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                           <  2:2593043/100‑1 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                         >  2:504571/1‑100 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                      >  1:2629182/1‑100 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCCGATGCGGCGTTAACGCCTTATTCGGCCTACAAAATCGTGCAAATTCAATAcct                                                                       >  1:3166936/1‑97 (MQ=255)
                                        cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                                  <  2:2073775/100‑1 (MQ=255)
                                         gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa                                                                 <  2:3612945/100‑1 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                                         <  2:277542/100‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑g                                                     <  2:410958/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑Gt                                                    <  1:338780/100‑1 (MQ=255)
                                                 aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑Gt                                                    <  1:98090/100‑1 (MQ=255)
                                                   aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAg                                                  <  2:1525862/100‑1 (MQ=255)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGcc                                               >  2:1630060/1‑100 (MQ=255)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGcc                                               <  1:1002677/100‑1 (MQ=255)
                                                       ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCt                                              <  1:2213313/100‑1 (MQ=255)
                                                       ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCt                                              >  2:3422694/1‑100 (MQ=255)
                                                        tgACAAAGTGAACTTTGTTTATGCCGGATGCGGCATTAACGCCTTATCCGGCCTACAAAATCGTGAAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAgcc                                                >  1:1433280/1‑95 (MQ=17)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGa                                            >  2:2610607/1‑100 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                          >  1:1070408/1‑100 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                          <  1:2688472/100‑1 (MQ=255)
                                                           cAAAGTGCGCTTGGTTTATGCCGGATGCGGCGTTCACGTCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATa                                          <  2:74374/100‑1 (MQ=255)
                                                                  tccgatctttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCATAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                   <  2:2254921/93‑1 (MQ=255)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                              >  1:3497688/1‑100 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            >  2:1459496/1‑100 (MQ=255)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:2676870/5‑100 (MQ=21)
                                                                           taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:3142398/4‑100 (MQ=18)
                                                                           ctggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAg                           <  1:1322762/97‑1 (MQ=21)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:1558833/3‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:1627716/3‑100 (MQ=18)
                                                                            aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                         >  2:2466547/1‑100 (MQ=255)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1375208/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2945391/2‑100 (MQ=14)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTC‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:3254941/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                        >  1:1164148/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:3187227/2‑99 (MQ=18)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                       >  1:851888/1‑100 (MQ=255)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATTTTGCAGAAAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:928360/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                       >  1:2134425/1‑98 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                       >  2:350184/1‑98 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGC‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  2:913008/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTT‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaa                       >  1:637443/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCac                       >  2:1199817/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  2:1328885/1‑99 (MQ=17)
                                                                               ccggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGTCTGATAAGCGTAGCGCATCAGGCagt                      >  1:710310/9‑98 (MQ=18)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  1:2334547/9‑98 (MQ=21)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTC‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:1290245/9‑98 (MQ=18)
                                                                                cggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgc                     >  2:1224423/8‑97 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:774148/7‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:478513/7‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:1898626/7‑100 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                    >  2:112746/7‑96 (MQ=25)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:3620743/7‑100 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                    >  2:64034/7‑100 (MQ=17)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  2:2706498/100‑1 (MQ=255)
                                                                                    tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 <  1:2659298/100‑1 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAATTTAAATACATTGCAGAATT‑‑‑TTC‑‑GTAATCATGATAAGCGGAGTGAATCAGGCATTcttgcct               >  1:2553076/1‑100 (MQ=18)
                                                                                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  1:3470014/1‑100 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  2:1755694/1‑100 (MQ=255)
                                                                                       ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt              >  1:2918787/1‑100 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  1:3033038/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            >  1:2974590/1‑100 (MQ=255)
                                                                                          gTTAACGCGTTATCCGGCCGACCGAATCGTGCAAATCCAAGACATTGCAGAATT‑‑‑TTC‑‑CTAGGCCTGATAAGCAGACCGCATCAGGCATTtttgcct               <  1:2812724/96‑1 (MQ=18)
                                                                                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          <  2:359905/100‑1 (MQ=255)
                                                                                           ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          <  2:851893/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  1:2484278/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATACGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  1:1947466/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  1:1243904/1‑100 (MQ=255)
                                                                                               cGCCGTATCCGGCCTACAAAATCGTGCAAATTCAAGACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  1:1051884/99‑1 (MQ=255)
                                                                                                  cTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc   >  2:230982/1‑100 (MQ=255)
                                                                                                   |                                                                                                         
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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: