Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,887 (GC)5→4 100% coding (400‑401/531 nt) hslV ← peptidase component of the HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,121,8790G.97.7% 594.7 / 5.4 131coding (409/531 nt)hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base . (55/73);  total (55/76)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01
*NC_0009134,121,8800C.96.2% 576.6 / 13.7 131coding (408/531 nt)hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base . (54/72);  total (55/76)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 161 total aligned reads displayed.

CGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGAT  >  NC_000913/4121780‑4121972
                                                                                                   ||                                                                                            
cGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAg                                                                                               <  2:3674787/100‑1 (MQ=255)
    ggTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcg                                                                                           <  1:3087839/100‑1 (MQ=255)
     gtgtGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgc                                                                                          >  2:2245097/1‑100 (MQ=255)
     gtgtGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGAGCTAAGTTCAGTGTTTTCTAACAgagcg                                                                                           >  1:3584031/1‑99 (MQ=255)
      tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAACAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                         <  2:2553472/100‑1 (MQ=255)
       gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgc                                                                                        >  1:674993/1‑100 (MQ=255)
       gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgc                                                                                        >  2:2749746/1‑100 (MQ=255)
       gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGCTTTCTAACAgcgcgcgc                                                                                        >  1:1268974/1‑100 (MQ=255)
        tgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgcc                                                                                       >  1:348098/1‑99 (MQ=255)
        tgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgcc                                                                                       >  1:596408/1‑99 (MQ=255)
        tgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgcc                                                                                       <  1:1957009/100‑2 (MQ=255)
         ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                      >  2:89913/1‑98 (MQ=255)
         ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                      >  2:342537/1‑98 (MQ=255)
            aaTGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTATAACAgcgcgcgccgca                                                                                    >  1:3243533/1‑95 (MQ=255)
              tggtttgtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTGTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgcagcc                                                                                 <  1:2216306/100‑8 (MQ=255)
              tggttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgcagcc                                                                                 >  2:2019087/1‑93 (MQ=255)
              tggttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgcagcc                                                                                 >  1:1470237/1‑93 (MQ=255)
               ggttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCt                                                                              >  2:1109474/1‑100 (MQ=255)
               ggttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCt                                                                              >  1:910697/1‑100 (MQ=255)
                attggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTg                                                                             >  1:1583107/2‑100 (MQ=255)
                 ttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTgg                                                                            <  1:3705808/100‑1 (MQ=255)
                 ttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTgg                                                                            <  1:536366/100‑1 (MQ=255)
                 ttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTgg                                                                            <  1:536692/100‑1 (MQ=255)
                  tggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTggg                                                                           >  2:3396984/1‑100 (MQ=255)
                   ggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGc                                                                          >  2:537457/1‑100 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                         <  2:145193/100‑1 (MQ=255)
                      ataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTa                                                                       >  2:188631/1‑100 (MQ=255)
                      ataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTGCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGTGCGCCGCAGCCTGGGCGTa                                                                       <  1:76287/100‑1 (MQ=255)
                       taGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTACAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTaa                                                                      >  2:2218380/1‑100 (MQ=255)
                         gATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAgg                                                                    >  2:1445637/1‑100 (MQ=255)
                         gATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAgg                                                                    >  2:2300926/1‑100 (MQ=255)
                           tGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgc                                                                  >  1:2337693/1‑100 (MQ=255)
                            gCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAggg                                                                   <  2:2472836/100‑1 (MQ=255)
                            gCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgcc                                                                 >  2:3539699/1‑100 (MQ=255)
                             cAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccg                                                                >  1:1086653/1‑100 (MQ=255)
                              aaaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAAC‑‑CGCGCGCGCCGCAGCCTGGGCGTAAGGgccgc                                                               >  2:3787663/1‑100 (MQ=255)
                              aaaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGCTCCGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccgc                                                               >  1:1829560/1‑100 (MQ=255)
                               aaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccgcc                                                              >  2:2744408/1‑100 (MQ=255)
                                  gTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGa                                                           <  1:1778247/100‑1 (MQ=255)
                                  gTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGa                                                           >  2:1180845/1‑100 (MQ=255)
                                    cGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGc                                                         <  2:2172629/100‑1 (MQ=255)
                                     gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        <  1:1320353/100‑1 (MQ=255)
                                     gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        <  1:1592625/100‑1 (MQ=255)
                                      ccTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCg                                                       >  2:3241922/1‑100 (MQ=255)
                                       cTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGa                                                      >  2:398792/1‑100 (MQ=255)
                                        tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     <  1:3249157/100‑1 (MQ=255)
                                         gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                    >  1:1906215/1‑100 (MQ=255)
                                         gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                    <  1:3104133/100‑1 (MQ=255)
                                          cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAg                                                   <  1:469602/100‑1 (MQ=255)
                                            atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCa                                                 >  2:446333/1‑100 (MQ=255)
                                              atCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAAt                                               >  1:2903602/1‑100 (MQ=255)
                                               tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATa                                              <  2:887756/100‑1 (MQ=255)
                                               tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATa                                              <  2:73040/100‑1 (MQ=255)
                                                ccAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATaa                                             >  2:2110004/1‑100 (MQ=255)
                                                 cAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATaa                                             <  1:3633095/99‑1 (MQ=255)
                                                 cAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAg                                            >  2:1800884/1‑100 (MQ=255)
                                                  aaCGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGa                                           <  1:617264/100‑1 (MQ=255)
                                                   aCGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGAt                                          <  1:2243169/100‑1 (MQ=255)
                                                   aCGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGAt                                          <  1:2726768/100‑1 (MQ=255)
                                                    cGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATc                                         <  1:1171891/100‑1 (MQ=255)
                                                     gCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCg                                        <  1:137811/100‑1 (MQ=255)
                                                       cTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGtt                                      >  1:2440777/1‑100 (MQ=255)
                                                         tttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGtttt                                    <  2:492597/100‑1 (MQ=255)
                                                          ttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTc                                   <  1:1921858/100‑1 (MQ=255)
                                                          ttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTc                                   <  1:2379301/100‑1 (MQ=255)
                                                          ttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTc                                   >  2:1750083/1‑100 (MQ=255)
                                                            cAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCt                                  >  1:1897318/1‑99 (MQ=255)
                                                            cAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                                 >  2:1586341/1‑100 (MQ=255)
                                                            cAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                                 <  2:503310/100‑1 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               >  2:1714322/1‑100 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               <  1:42916/100‑1 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               <  2:2403486/100‑1 (MQ=255)
                                                               cAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCt                              <  2:1118518/100‑1 (MQ=255)
                                                                 aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGc                            >  2:1841349/1‑100 (MQ=255)
                                                                   ttCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcac                          >  2:3644218/1‑100 (MQ=255)
                                                                     cACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacca                        >  1:271425/1‑100 (MQ=255)
                                                                     cACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGATTTCTGGCTGcacca                        >  1:987941/1‑100 (MQ=255)
                                                                         ggCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGTAGCCGATAGCAATAAGATCGTTTTCGGGCTGCACCACGTc                    <  2:3216998/100‑1 (MQ=255)
                                                                         ggCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTc                    >  1:2045941/1‑100 (MQ=255)
                                                                         ggCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTc                    <  1:2629992/100‑1 (MQ=255)
                                                                         ggCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTc                    <  2:3350972/100‑1 (MQ=255)
                                                                          gctcTAAGTGCAGCTTTTTCTACCACCGCGTGCCGCAGCCTGGGCGTAAGGGTCGCACGAGCCGATAGCAATAAGATCGGTTTCTGGCTGcaccta                       <  1:720130/93‑3 (MQ=255)
                                                                             cTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGGTTTCTGGCTGCACCACATCACCg                >  2:3025377/1‑100 (MQ=255)
                                                                                 gTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAc            >  1:1532470/1‑100 (MQ=255)
                                                                                  ttCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAcc           <  2:1980610/100‑1 (MQ=255)
                                                                                   tCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCg          <  2:3097475/100‑1 (MQ=255)
                                                                                   tCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCg          <  1:1002048/100‑1 (MQ=255)
                                                                                   tCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCg          >  2:3004936/1‑100 (MQ=255)
                                                                                    cAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCgg         <  1:1109469/100‑1 (MQ=255)
                                                                                     aGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt        <  1:3200497/100‑1 (MQ=255)
                                                                                     aGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt        <  1:883635/100‑1 (MQ=255)
                                                                                     aGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt        <  1:2515467/100‑1 (MQ=255)
                                                                                      gtgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg       >  2:560481/1‑100 (MQ=255)
                                                                                       tgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtga      <  2:674997/100‑1 (MQ=255)
                                                                                       tgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtga      >  1:286323/1‑100 (MQ=255)
                                                                                         ttttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatg    <  2:1906222/100‑1 (MQ=255)
                                                                                          tttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga   <  2:1796930/100‑1 (MQ=255)
                                                                                          tttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacgt   <  1:1714314/100‑4 (MQ=255)
                                                                                           ttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat     <  1:2153294/97‑1 (MQ=255)
                                                                                             ttatcaccgcgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCGGGCTGCACCACGTCTCCGTTACCGGTgatgat  <  1:1619055/93‑1 (MQ=255)
                                                                                                   ||                                                                                            
CGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGAT  >  NC_000913/4121780‑4121972

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