Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,446 A→T 15.9% intergenic (‑35/+24) dtpC ← / ← cadA dipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,4460AT15.9% 582.7 / 91.4 234intergenic (‑35/+24)dtpC/cadAdipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (112/83);  new base T (17/20);  total (131/103)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-02

Only 100 of 239 total aligned reads displayed.

AGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCG  >  NC_000913/4356348‑4356545
                                                                                                  |                                                                                                   
  tACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                  <  2:2928093/100‑1 (MQ=255)
  tACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTAAAAAAAAAAAAAAAAGGGAAGTGGCAAGcc                                                                                                  <  1:1580005/100‑1 (MQ=255)
    cTCCCAGATTTGGATCGCCACGATATAGTATATCTCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                                >  2:1967741/1‑100 (MQ=255)
    cTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                  <  2:2530447/98‑1 (MQ=255)
    cTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                  <  1:2725514/98‑1 (MQ=255)
    cTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                                >  1:1474716/1‑100 (MQ=255)
       ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTc                                                                                             <  2:3577109/100‑1 (MQ=255)
        cAGATTTGGATCGCCACGATATAGTCTATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGTGCAAAAAAGGGACGTGGCAAGCCCCTTcc                                                                                            >  1:2481498/1‑100 (MQ=255)
         aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTTCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTccc                                                                                           <  1:2354458/100‑1 (MQ=255)
         aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                           <  2:2807134/100‑1 (MQ=255)
         aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                           <  2:3965818/100‑1 (MQ=255)
         aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTACTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                           <  2:4066702/100‑1 (MQ=255)
         aGATTGGGATCGCCACGATATAGAATATCGCGCGCGGATGTGAGGGTGTTTTCATGTGTTCTCATAATAAAAAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                           <  1:2960374/100‑1 (MQ=255)
          gATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCt                                                                                          <  2:1069143/100‑1 (MQ=255)
          gATTGGGAGCGTCGCGATATATTATATCGCGCTCGGCTGTGGGTTTGTTTTCATGTTTTCTCCTTATGATCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                           <  1:3640011/99‑1 (MQ=255)
           aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtt                                                                                         <  2:211154/100‑1 (MQ=255)
            tttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTg                                                                                        <  1:4288216/100‑1 (MQ=255)
            tttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTAAGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtt                                                                                         <  1:3471370/99‑1 (MQ=39)
             ttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttt                                                                                       <  2:1689106/100‑3 (MQ=255)
             ttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGt                                                                                       <  2:459425/100‑1 (MQ=255)
              tGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                                                                      >  1:3972193/1‑100 (MQ=255)
              tGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                                                                      <  2:450046/100‑1 (MQ=255)
               ggATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttt                                                                                       <  1:3906974/98‑3 (MQ=255)
                 aTCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttttgc                                                                                   <  1:2022852/100‑7 (MQ=38)
                  cgcgcgcGATATACTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAg                                                                                  >  2:1616368/7‑100 (MQ=255)
                    gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCt                                                                                >  2:3285866/1‑100 (MQ=255)
                    gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCt                                                                                >  1:577755/1‑100 (MQ=255)
                     ccACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTa                                                                               >  2:724117/1‑100 (MQ=255)
                     ccACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTa                                                                               >  2:1161935/1‑100 (MQ=255)
                         gATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatta                                                                           <  1:631310/100‑1 (MQ=255)
                              aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttttt                                                                      <  1:46267/100‑1 (MQ=255)
                              aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttttt                                                                      >  1:1870801/1‑100 (MQ=255)
                              aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAttt                                                                        >  1:1693525/1‑98 (MQ=255)
                                 atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGc                                                                   <  2:2394897/100‑1 (MQ=255)
                                 atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGc                                                                   <  2:492910/100‑1 (MQ=255)
                                 atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGc                                                                   <  1:1807397/100‑1 (MQ=255)
                                  tatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCt                                                                  >  1:3304452/1‑100 (MQ=255)
                                  tatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTACt                                                                  >  2:2053070/1‑100 (MQ=255)
                                   atCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCtt                                                                 >  1:1207541/1‑100 (MQ=255)
                                     cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCttt                                                                <  1:2399051/99‑1 (MQ=255)
                                     cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttc                                                               <  1:1761037/100‑1 (MQ=255)
                                        gcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttcttc                                                            <  1:1526570/100‑1 (MQ=255)
                                           cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                         >  2:1718427/1‑100 (MQ=255)
                                           cgGCTGTGAGGGTGGTTTCATGTGGTCTCCTTATGAGCAAAAAAGGGAAGTGGTAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                         <  1:1370560/100‑1 (MQ=255)
                                            ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                        <  2:2760107/100‑1 (MQ=255)
                                             gCTGTGAGATTGTTTTCATGTGCTCTCCTTATGAGCAAAAAATTGAAGTGTCAAGCCACTTCCCTTCTACGAGCTAATTATTCTTTGCTTTCTTCTTTCa                                                       <  1:3737769/100‑1 (MQ=255)
                                              cTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                      <  1:2876203/100‑1 (MQ=255)
                                                gtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATa                                                    >  2:1760303/1‑100 (MQ=255)
                                                 tgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCACGCCACTTCCCTTGTCCGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAc                                                   >  1:4202349/1‑100 (MQ=255)
                                                  gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                  <  2:3149376/100‑1 (MQ=255)
                                                  gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                  >  1:3943521/1‑100 (MQ=255)
                                                  gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                  <  2:2978041/100‑1 (MQ=255)
                                                    gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCtt                                                <  2:1927303/100‑1 (MQ=255)
                                                    gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCtt                                                <  2:754244/100‑1 (MQ=255)
                                                       tgtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                             >  2:4369704/1‑100 (MQ=255)
                                                        gtTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                             >  1:2595833/1‑99 (MQ=255)
                                                        gtTTGCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTTGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCGTTCAATACCTTAACg                                            <  1:1002663/100‑1 (MQ=255)
                                                         ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACgg                                           <  2:1319836/100‑1 (MQ=255)
                                                         ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACgg                                           <  2:630603/100‑1 (MQ=255)
                                                           ttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGta                                         <  1:3243887/100‑1 (MQ=255)
                                                                 tgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCgg                                   >  2:3537227/1‑100 (MQ=255)
                                                                 tgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCgg                                   >  1:3251078/1‑100 (MQ=255)
                                                                    tctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCa                                <  2:3504694/100‑1 (MQ=255)
                                                                     ctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCAt                               >  2:3423139/1‑100 (MQ=255)
                                                                       ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                             <  1:3063240/100‑1 (MQ=255)
                                                                       ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                             >  1:3443919/1‑100 (MQ=255)
                                                                        cTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAg                            >  2:3119555/1‑100 (MQ=255)
                                                                        cTTATGAGCAAAAAAAAAAAAAAGAAAACCACTTCACAAGTACAAGCTAATTATTTTTTACTTACTTATTTCAATACATTAACGGtataac                                     >  1:1318491/1‑89 (MQ=255)
                                                                         ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                           >  1:2045906/1‑100 (MQ=255)
                                                                          tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                          <  2:2702153/100‑1 (MQ=255)
                                                                          tATGAGCAAAAAAAGGAAGTGGCAAGCCACTTCCCTTGTACGAGCGAATTAGTTTTTGCTTTCTTCTTGCAATACCTTAACGCTATAGCGGCCATCAGcc                          <  2:3833389/100‑1 (MQ=255)
                                                                          tATGAAAAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                          <  1:696580/100‑1 (MQ=255)
                                                                            tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCt                         >  1:3871314/1‑99 (MQ=255)
                                                                             gAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                        >  2:2298842/1‑99 (MQ=255)
                                                                             gAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGa                       >  2:1417876/1‑100 (MQ=255)
                                                                             gAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGATAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCAATCAGCCTGa                       >  2:1239857/1‑100 (MQ=255)
                                                                             gAGCAAAAAAGGGAAGAGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGGTTTCTTCTTTCAATACCTTAACGGTATAGTGGCCATCAGCGTg                        >  1:4301500/1‑99 (MQ=255)
                                                                               tcgtgcaaGGGAAGTGGCTTTCAACTTCCCTTGGAGGGGCTAATTATTTTTTGCTTTATACTTTCAATACATTACCGGTATAGCGGCCATCAGCCTGACg                     <  1:2515086/94‑1 (MQ=255)
                                                                               gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                     >  1:3178168/1‑100 (MQ=255)
                                                                               gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                     >  1:3258697/1‑100 (MQ=255)
                                                                                cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                     >  2:4233406/1‑99 (MQ=255)
                                                                                 aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                   <  2:333574/100‑1 (MQ=255)
                                                                                 aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                   <  1:893658/100‑1 (MQ=255)
                                                                                  aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTCTTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTa                  >  2:3564139/1‑100 (MQ=255)
                                                                                  aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTTTTTTTTTTTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTa                  >  1:1056883/1‑100 (MQ=255)
                                                                                    caaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg                >  2:1490665/2‑100 (MQ=255)
                                                                                       gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc             >  2:4025352/1‑100 (MQ=255)
                                                                                       gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc             >  2:157314/1‑100 (MQ=255)
                                                                                       gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc             >  2:2579401/1‑100 (MQ=255)
                                                                                       gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTATTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc             >  2:4113780/1‑100 (MQ=255)
                                                                                        ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTTTTTTTTTCTTTCAATACACTAACGGTATAGAATCCATTAGCCTGACGGTAGGcccc            >  1:2183505/1‑100 (MQ=255)
                                                                                        ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGTTTTTTTTTTTCAATACATTAACGGTATAGCGGCCATCAGCTTGACGGTATGCAcc            >  1:4397362/1‑100 (MQ=255)
                                                                                        ggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc            <  2:1285459/100‑1 (MQ=255)
                                                                                        ggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc            <  1:1711710/100‑1 (MQ=255)
                                                                                          aaGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGt          <  1:714109/100‑1 (MQ=255)
                                                                                              ggCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAat      >  2:2609801/1‑100 (MQ=255)
                                                                                              ggCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAat      >  2:3482390/1‑100 (MQ=255)
                                                                                                cAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAatat    >  2:2183969/1‑100 (MQ=255)
                                                                                                  aGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCg  >  2:3896985/1‑100 (MQ=255)
                                                                                                  aGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCg  <  2:431823/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
AGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCG  >  NC_000913/4356348‑4356545

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