Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 24.3% 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,9850TC24.3% 517.8 / 168.6 259intergenic (‑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 (88/108);  new base C (60/3);  total (148/111)
Fisher's exact test for biased strand distribution p-value = 3.94e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-02

Only 100 of 264 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084
                                                                                                  |                                                                                                              
 ctgagggATGGCGCTGCACGCAGAGGCAGCGCGTAAACGGGTGTACAGATTGTTGTCAAAGTGCGATTTGTGTATGGCGTAGCCCGCGTTAACACCttc                                                                                                               <  2:362567/94‑2 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                          <  2:2211313/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                        <  1:3776470/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCGCATCCGGc                                                                                                         <  2:1624355/99‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                        <  1:4138406/100‑2 (MQ=37)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGCGCGCTTCGTTTATGCCCGATGCGGCGCTAACGCCTTCCCCGGCCTa                                                                                                      >  2:1857371/1‑100 (MQ=25)
          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                                     >  2:348834/1‑100 (MQ=255)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGGTAACGCCTTATCCGGCCTACaa                                                                                                  >  1:105899/1‑100 (MQ=255)
               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                               <  1:3178730/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                             <  1:3737482/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                           >  1:2927163/1‑100 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                          <  2:3907159/100‑1 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                         <  1:3900806/100‑1 (MQ=255)
                         ggCCCAGCCCCACCGAGTATACATATCGTTGACAAAGTGCGCTTTGTTTGTGCCGGAGGCGGCGTTACCGACTTATCCGGCCTACACAACCGTGCa                                                                                         <  2:3158118/96‑1 (MQ=17)
                         agCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                                     >  2:1754591/2‑100 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                                    <  2:1387920/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                                 <  2:3087705/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                                 <  1:144595/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                                <  1:3314644/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCAGCCTACAAAATCGTGCAAATTCAATa                                                                                <  1:2713647/100‑1 (MQ=255)
                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                               >  1:2310467/1‑100 (MQ=255)
                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                               <  1:3291938/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCGACAAAATCGTGCAAATTCAATACATTGCAg                                                                      <  2:1675371/100‑1 (MQ=255)
                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGAATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                                           <  2:2903065/97‑1 (MQ=255)
                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa                                                                     <  1:2768258/100‑1 (MQ=255)
                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTCCAAAATCGTGCAAATTCAATACATTGc                                                                          >  2:2491604/1‑96 (MQ=255)
                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                                    >  1:919460/1‑100 (MQ=255)
                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑t‑‑‑t                                                            <  1:2488960/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGTATTCGGCGTTAACTCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑t‑‑‑tt                                                          >  1:2268689/1‑100 (MQ=255)
                                              cgatctgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑g                                                      <  2:2288436/95‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTa                                                    >  1:2568812/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTa                                                    <  2:4151717/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTa                                                    >  1:3397320/1‑100 (MQ=255)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGc                                                 >  2:2981095/1‑100 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGcc                                                <  1:2124849/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGcc                                                >  1:4052532/1‑100 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTg                                              >  2:2378451/1‑100 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTg                                              >  2:857927/1‑100 (MQ=255)
                                                      cgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTg                                              >  2:281776/2‑100 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATa                                           <  1:871671/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATa                                           <  1:4047399/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATa                                           <  2:2646751/100‑1 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATaa                                          >  2:3743007/1‑100 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCg                                       >  2:2916819/1‑100 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                      >  2:75980/1‑100 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                     >  1:233793/1‑99 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                    <  2:3729640/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  1:281779/100‑1 (MQ=255)
                                                                   attttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAATCGTGCTAATTCAATATATTGCAGAAA‑‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAgcgc                                 >  1:1406379/10‑100 (MQ=17)
                                                                    ttGTTTATGCCGGATGCGCCGTTAACGCCTTATCCGGCCTACAAACTCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                <  1:1157764/100‑1 (MQ=255)
                                                                    atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                               >  1:2671933/9‑100 (MQ=18)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                               <  2:1924994/100‑1 (MQ=255)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAG‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCa                                >  1:1294540/7‑98 (MQ=17)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                             <  2:2395782/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCTACAAGAACGTGCAAATTCAACAAATTGTGATGCAT‑‑C‑‑‑CT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                             >  1:1023285/5‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATA‑‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:3794943/4‑100 (MQ=18)
                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTCCAGAAT‑‑T‑‑‑TTC‑‑GTAGACCTGATAAGAGGAGCGCATCAgg                           <  2:2567147/100‑1 (MQ=255)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAA‑‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:3633099/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAT‑‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                           >  1:3687641/3‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGC‑‑C‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:3332077/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCT‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:1516650/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:3830939/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2057983/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCT‑‑C‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:3660336/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                          >  1:1366761/2‑99 (MQ=18)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  2:933234/1‑100 (MQ=255)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  1:1754180/1‑100 (MQ=255)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAA‑‑T‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  1:1641078/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  1:743393/1‑100 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGAT‑‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:3517283/1‑99 (MQ=21)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAAT‑‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:3561992/9‑98 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑T‑‑‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagc                       >  1:2655868/9‑98 (MQ=18)
                                                                             ccGGATGCGGAGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGACCTGATAAGCTGAGCGCATAACGCAtt                       >  2:2547629/1‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAAT‑‑C‑‑‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                    >  1:4005093/7‑98 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAA‑‑‑CTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                      >  1:2392307/7‑97 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                     >  1:2524187/7‑96 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:535835/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAAC‑‑G‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                     >  1:3390970/7‑96 (MQ=25)
                                                                               ggacaaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                     >  2:1197906/7‑96 (MQ=17)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑‑C‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                    >  2:2828621/6‑100 (MQ=18)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:3438729/100‑1 (MQ=255)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:2947696/100‑1 (MQ=255)
                                                                                    aggcgTTAACGCCTTTTTTGGCCTACTAAATCATGCAAATTCAATTCATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTAATAAGCGGAGCGCATCAGGCATTtttgc                  >  2:3554594/2‑98 (MQ=18)
                                                                                     ngcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATCtttgcctt               >  2:1041620/2‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               <  2:11452/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  2:320447/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  2:481015/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGACCAACAAATTATGCAACTCCACAATATTGCAGAAA‑‑T‑‑‑TAA‑‑GTAGGACTTATAAGCATAACGCATAAAACAtttattcgttt              >  2:329578/1‑90 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTT‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  2:883385/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  1:3215199/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  2:1096380/100‑1 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt         >  1:3365096/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:1718252/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAACGCATCAGGCATTTTTGCCTTTGCCATCa       >  2:3967615/1‑100 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGGCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  1:3960608/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCan      <  2:535133/100‑2 (MQ=255)
                                                                                               ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  1:3105986/1‑100 (MQ=255)
                                                                                               ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  1:3015510/1‑100 (MQ=255)
                                                                                                 ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  1:441288/100‑1 (MQ=255)
                                                                                                  tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  <  1:1890908/100‑1 (MQ=255)
                                                                                                  |                                                                                                              
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084

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