Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 33.7% 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,9850TC33.7% 427.9 / 285.3 291intergenic (‑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 (82/111);  new base C (92/6);  total (174/117)
Fisher's exact test for biased strand distribution p-value = 2.63e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-01

Only 100 of 296 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084
                                                                                                  |                                                                                                        
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                         <  1:532975/100‑1 (MQ=255)
 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                        <  2:1280650/100‑1 (MQ=255)
    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                     >  2:634613/1‑100 (MQ=255)
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                      >  2:276027/1‑98 (MQ=255)
      ggATGGCGCTGCACGCAGACGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTGGTTTATGCCGTTTGCGGCTTTAACGCCtttttcg                                                                                                     >  2:1585974/1‑93 (MQ=25)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                  <  1:1761329/100‑1 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                  <  2:3391859/99‑1 (MQ=255)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑aa                                                                                           >  1:1970849/1‑100 (MQ=255)
               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAAt                                                                                        <  1:3629351/100‑1 (MQ=255)
                 gaCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACACAGTACGCTATGTTTATGCCGAATGCGTCTTTAACGCCTTATCCGTCCT‑Acgca                                                                                           >  2:1485222/2‑93 (MQ=25)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                                      <  1:3017425/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                                      <  2:3814777/100‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCg                                                                                      <  1:2869382/100‑1 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGt                                                                                     <  2:834424/100‑1 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGt                                                                                     >  2:1197579/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGt                                                                                     >  2:1686935/1‑99 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTg                                                                                    <  2:814803/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGc                                                                                   >  1:3551049/1‑100 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAAtt                                                                              <  1:3346184/100‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCa                                                                            >  1:2314827/1‑100 (MQ=255)
                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCaa                                                                           <  1:1145316/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACTAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGTCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAAt                                                                          <  2:2364817/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATa                                                                         <  1:377823/100‑1 (MQ=255)
                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTg                                                                    <  2:1773465/100‑1 (MQ=255)
                                       gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑Ga                                                              >  1:819259/1‑100 (MQ=255)
                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑‑AAATCGTGCAAATTCAATACATTG‑CA‑GAAt‑t                                                          >  1:1377222/1‑100 (MQ=255)
                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAt‑tt                                                         <  2:433215/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAt‑ttt                                                        >  1:1104764/1‑100 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAt‑ttt                                                        >  2:1259000/1‑100 (MQ=255)
                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTc                                                       <  2:296087/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCg                                                      <  2:3715525/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCg                                                      <  2:2290613/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCg                                                      <  2:1370813/100‑1 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAgg                                                  <  1:2593722/100‑1 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAgg                                                  <  1:698438/100‑1 (MQ=255)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGc                                                 <  2:2698205/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGcc                                                <  2:2211694/100‑1 (MQ=255)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCt                                               <  1:1896689/100‑1 (MQ=255)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCt                                               <  1:2713677/100‑1 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTg                                              >  1:2901318/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGa                                             >  2:2556248/1‑100 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATa                                           <  1:88539/100‑1 (MQ=255)
                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAg                                         >  1:302455/1‑100 (MQ=255)
                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAg                                         <  1:1216033/100‑1 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGc                                        <  1:723624/100‑1 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCg                                       <  1:2122398/100‑1 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCgg                                      <  1:239360/100‑1 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCgg                                      <  2:1329633/100‑1 (MQ=255)
                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATTGTGTAAATTCAATAAATTG‑AT‑CGAT‑CCACGTAGGCCTGATAAGCGTAGCGCAt                               >  2:2006015/8‑100 (MQ=21)
                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAC‑AAAAACCTGCAAATTCAATATATTG‑CA‑GGAT‑ACGTGTAGGCCTGATAAGCGTAGCGCa                                >  1:3159563/8‑100 (MQ=12)
                                                                     cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAGTTCGTGCAAATTCAATAAATTG‑CA‑ATA‑‑TGACGTAGGCCTGATAAGCGTAGCGCATc                              >  1:909696/8‑100 (MQ=18)
                                                                      ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCa                             >  2:3139033/7‑100 (MQ=18)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CCAT‑CCTCGTAGGCCTGATAAGCGTAGCGCATc                              >  1:3864568/7‑100 (MQ=21)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                              <  2:3528445/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  1:3694331/100‑1 (MQ=255)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                            <  1:2980829/100‑1 (MQ=255)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAATTGTGTAAATTCAATAAATTG‑AT‑CGAT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAgg                           <  1:2006025/97‑1 (MQ=21)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:3274975/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAGCGTGCAAATTCAATACATTG‑CATGGGC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:2734088/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1742172/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAGTTCGTGCAAATTCAATAAATTG‑CA‑‑ATA‑TGACGTAGGCCTGATAAGCGTAGCGCATCAGGcg                         >  2:1141131/3‑99 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCAT‑A‑‑AACAACGCGCAAATTCAATATATTG‑CG‑GAGA‑AAATGTAGGCCTGATAAGCGTAGCGCATCAGGct                         >  2:1960640/3‑99 (MQ=9)
                                                                          aggCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                          >  1:1418089/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGAGATCACGCCGCATCCGGCCT‑AC‑AATACTGTGCAAATTCAATATATTG‑CA‑GGAA‑ACACGTAGGCCTGATAAGCCAAGCGCATCAGGc                          >  2:2434693/3‑100 (MQ=17)
                                                                           ggtcggataaggcgTTCACGCCGCATCCGACCC‑AC‑AAAGTCTTGCAAAATCAATAAATTG‑CA‑GGAA‑CGATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:1897403/11‑100 (MQ=21)
                                                                           ggCCGGATAAGGGGCTTACGCCGCATCCGGCAT‑AC‑AAGAACGTGCAAATTCAACAAATTGTGATGCAT‑CCT‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:746507/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GAGT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:840017/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GAGT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:1631156/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAGATCGTGCAAATTCAATATATTG‑CA‑TGAG‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:20689/2‑100 (MQ=17)
                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GACT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  2:1059684/1‑100 (MQ=255)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:1535773/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:2328186/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAGTCGTATAAATTCAATGTATTG‑CA‑GTTT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                        >  1:1590760/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAGCACGCAAACTCAATATATTG‑CA‑AAGA‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:2144583/1‑99 (MQ=14)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑AC‑AAAAACCTGCAAATTCAATATATTG‑CA‑GGAT‑ACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  1:3472357/1‑99 (MQ=17)
                                                                             ccggataaggcgTTTACGCCTTATCCGGCCT‑AC‑AAAACAGTGCAAATTCAATATATTG‑CA‑ATCT‑TCTGGTAGGCCTGATAAGCATCGCGCATCAGGCaag                       >  1:3911787/9‑98 (MQ=18)
                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCCT‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAACGCATCAGGc                          >  2:1849155/1‑97 (MQ=17)
                                                                              cggataaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAACAGTGCAAATTCAATATATTG‑CA‑ATCT‑TCTGGTAGGCCTGATAAGCATCGCGCATCAGGCaagc                      >  1:1637391/8‑97 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCTGCCT‑AC‑AAAAGCACGCAAATTCAATATATTG‑CA‑GAGA‑TGGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt                     >  1:2191830/7‑100 (MQ=14)
                                                                               ggataaggcgTTCACGCCGCATCCTGCCT‑AC‑AAAAGCACGCAAATTCAATATATTG‑CA‑GAGA‑TGGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt                     >  1:2198275/7‑100 (MQ=14)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:1516479/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑AC‑‑AAATCATGCAAATCCAAAATATTG‑CA‑GGAG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                    >  1:307670/7‑98 (MQ=21)
                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                     <  2:3883383/100‑1 (MQ=255)
                                                                                gataaggcgTTAACGCCGCATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                    >  2:2615176/6‑100 (MQ=18)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:1883779/100‑1 (MQ=255)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:302455/100‑1 (MQ=255)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:3804108/100‑1 (MQ=255)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:3235920/100‑1 (MQ=255)
                                                                                  taaggcgTTTACGCCGCATCCGGCCT‑AC‑AAAAGCGTGCAAATTCAATACATTG‑CA‑TGGG‑CCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattttac                  >  2:1429123/4‑98 (MQ=18)
                                                                                   gcggcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                 <  1:985615/100‑1 (MQ=255)
                                                                                      gcgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt              >  2:2926462/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  2:3595180/100‑1 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  1:1820633/100‑1 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  1:3227671/100‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          <  1:1116322/100‑1 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt         >  1:2086203/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  1:564073/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  2:1502175/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:248586/100‑1 (MQ=255)
                                                                                                  tATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  1:128132/1‑100 (MQ=255)
                                                                                                  |                                                                                                        
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑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: ‑