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.92.7% 498.4 / 31.5 123coding (409/531 nt)hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base G (3/6);  new base . (54/60);  total (57/66)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01
*NC_0009134,121,8800C.92.7% 502.7 / 31.4 123coding (408/531 nt)hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base C (3/6);  new base . (54/60);  total (57/66)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 137 total aligned reads displayed.

CGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGAT  >  NC_000913/4121780‑4121972
                                                                                                   ||                                                                                            
cGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAg                                                                                               >  1:1525955/1‑100 (MQ=255)
cGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAg                                                                                               >  2:2839215/1‑100 (MQ=255)
     gtgtGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgc                                                                                          >  1:2408049/1‑100 (MQ=255)
      tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                         >  2:1214088/1‑100 (MQ=255)
       gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTCACAgcac                                                                                            >  1:1064688/1‑94 (MQ=255)
       gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgc                                                                                        >  2:1625372/1‑100 (MQ=255)
       gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgc                                                                                        >  2:1009006/1‑100 (MQ=255)
        tgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgcc                                                                                       <  2:1792354/100‑2 (MQ=255)
         ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                      <  2:2642268/100‑3 (MQ=255)
          gAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGTCGCTAAGGTCAGTGTTTTCTAACAgcgcgcgccgc                                                                                     >  1:701942/1‑97 (MQ=255)
           aaaTGGTTGGTAAAGATGCAAATGTCGCCTGCAATAGCCAACGCCTTTTCAGCAATGTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgca                                                                                    <  1:2706333/100‑5 (MQ=255)
              tggttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgcagcc                                                                                 <  2:1265069/100‑8 (MQ=255)
              tggttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgcagcc                                                                                 <  1:2693125/100‑8 (MQ=255)
                gttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTg                                                                             >  1:609535/1‑100 (MQ=255)
                 ttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTgg                                                                            <  2:1001652/100‑1 (MQ=255)
                   ggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAGTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTgg                                                                            >  1:2738269/1‑100 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                         <  2:2397475/100‑1 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                         <  2:2314079/100‑1 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                         <  1:583369/100‑1 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                         <  2:2756203/100‑1 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                         <  1:1237030/100‑1 (MQ=255)
                      ataGATGCAAATGTCGCCTGCATTATCCAACGCCTTTGCAGCAATTTTACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCg                                                                         <  1:1963865/100‑1 (MQ=255)
                      ataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTa                                                                       >  1:45256/1‑100 (MQ=255)
                      ataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTa                                                                       <  2:45257/100‑1 (MQ=255)
                      ataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTa                                                                       >  1:903294/1‑100 (MQ=255)
                       taGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTaa                                                                      >  1:2110888/1‑100 (MQ=255)
                        aGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAg                                                                     >  2:2541718/1‑100 (MQ=255)
                        aGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAg                                                                     <  2:534412/100‑1 (MQ=255)
                         gATGCAAGTGTCGCCGGCAATATCCAACGCCTTTTCCGCAATTTCACGGGCGCTAAGTTCGGTGTTTTCTAACA‑‑TCGCGCGCCGCAGCGTGGGCTTa                                                                       >  1:719101/1‑97 (MQ=255)
                           tGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgc                                                                  >  1:104510/1‑100 (MQ=255)
                           tGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgc                                                                  >  1:2202972/1‑100 (MQ=255)
                            gCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgcc                                                                 <  2:2366598/100‑1 (MQ=255)
                            gCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgcc                                                                 <  2:2781764/100‑1 (MQ=255)
                             cAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccg                                                                >  2:1312494/1‑100 (MQ=255)
                               aaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCGAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccgcc                                                              <  2:890653/100‑1 (MQ=255)
                                 tGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTTGGCGTAAGGGCCGCCgg                                                            >  1:1736686/1‑100 (MQ=255)
                                 tGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCGAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCgg                                                            <  1:1844275/100‑1 (MQ=255)
                                    cGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGc                                                         >  2:2097911/1‑100 (MQ=255)
                                     gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        <  1:249255/100‑1 (MQ=255)
                                     gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        <  2:1410069/100‑1 (MQ=255)
                                     gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        <  2:1644778/100‑1 (MQ=255)
                                     gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        >  2:2693723/1‑100 (MQ=255)
                                     gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        >  1:74259/1‑100 (MQ=255)
                                      ccTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTACGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGACCCg                                                       >  1:2516706/1‑100 (MQ=255)
                                        tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     <  2:967717/100‑1 (MQ=255)
                                        tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     >  1:1971766/1‑100 (MQ=255)
                                        tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     <  1:1647959/100‑1 (MQ=255)
                                        tGCAATATCCAACACCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGCGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     >  1:2122187/1‑100 (MQ=255)
                                         gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                    >  2:481696/1‑100 (MQ=255)
                                         gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                    >  2:2845737/1‑100 (MQ=255)
                                         gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTGCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                    <  2:1754272/100‑1 (MQ=255)
                                          cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAg                                                   <  1:2195326/100‑1 (MQ=255)
                                           aaTATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGc                                                  <  2:543526/100‑1 (MQ=255)
                                           aaTATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGc                                                  <  1:2162642/100‑1 (MQ=255)
                                            atatCCAACGCCTTTTTTTTTTTTTTTCGGGCGCTAAGATCAGTGATTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGc                                                  >  1:1396043/1‑99 (MQ=255)
                                            atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAg                                                   <  2:2000151/100‑1 (MQ=255)
                                            atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCa                                                 <  2:2188480/100‑1 (MQ=255)
                                            atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCa                                                 >  1:141175/1‑100 (MQ=255)
                                                ccAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCTCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATaa                                             >  1:2031667/1‑100 (MQ=255)
                                                ccAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATaa                                             <  2:1327140/100‑1 (MQ=255)
                                                 cAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGCTAGCAATAAg                                            >  1:851215/1‑100 (MQ=255)
                                                   aCGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGAt                                          <  1:506720/100‑1 (MQ=255)
                                                     gCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCg                                        <  2:2515233/100‑1 (MQ=255)
                                                        ttttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGttt                                     <  1:2079362/100‑1 (MQ=255)
                                                         tttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGtttt                                    <  2:870968/100‑1 (MQ=255)
                                                            cAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                                 >  1:1710982/1‑100 (MQ=255)
                                                            cAGCAAGTTCACGGGCGCTAAGTTCAGTGTTTGCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                                 <  1:1315420/100‑1 (MQ=255)
                                                             aGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCt                                  <  1:1281124/100‑1 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               <  1:1059551/100‑1 (MQ=255)
                                                                aaTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTg                             <  2:1024660/100‑1 (MQ=255)
                                                                 aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGc                            >  2:1312832/1‑100 (MQ=255)
                                                                 aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGc                            >  2:2733639/1‑100 (MQ=255)
                                                                 aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATAGTTTTCTGGCTGc                            <  1:1801067/100‑1 (MQ=255)
                                                                  tttCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGca                           <  2:337437/100‑1 (MQ=255)
                                                                   ttCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcac                          <  1:2527550/100‑1 (MQ=255)
                                                                    tCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacc                         <  2:2192180/100‑1 (MQ=255)
                                                                    tCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacc                         >  2:366804/1‑100 (MQ=255)
                                                                    tCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacc                         <  1:2387308/100‑1 (MQ=255)
                                                                     cACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCTGACTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcac                          >  1:1588983/1‑98 (MQ=255)
                                                                      aCGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcaccac                       <  1:443942/100‑1 (MQ=255)
                                                                       cGGGCGCTAAGTTCAGTGTTTTCTACCA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACg                      <  2:1583951/100‑1 (MQ=255)
                                                                       cGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACg                      <  2:1796445/100‑1 (MQ=255)
                                                                       cGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACg                      <  1:1835907/100‑1 (MQ=255)
                                                                       cGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACg                      <  2:1915934/100‑1 (MQ=255)
                                                                       cGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACg                      <  2:405203/100‑1 (MQ=255)
                                                                         ggCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTc                    <  2:903296/100‑1 (MQ=255)
                                                                           cgcCAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCAc                  >  1:147225/1‑100 (MQ=255)
                                                                            gcTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCAcc                 >  2:2880839/1‑100 (MQ=255)
                                                                                 gTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCCATAAGATCATTTTCTGGCTGCACCACGTCAACGTTa             >  2:865395/1‑99 (MQ=255)
                                                                                  ttCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAcc           <  2:2615926/100‑1 (MQ=255)
                                                                                   tCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCg          <  2:2285266/100‑1 (MQ=255)
                                                                                      gtgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAGGGGCCGCCGGAGCCGATAGCAATAAGATTGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg       <  2:1929139/100‑1 (MQ=255)
                                                                                        gtTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg       >  2:2210343/1‑100 (MQ=255)
                                                                                        gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat     <  1:1944420/100‑1 (MQ=255)
                                                                                        gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat     >  2:2378707/1‑100 (MQ=255)
                                                                                        gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat     >  1:950579/1‑100 (MQ=255)
                                                                                          tttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga   >  1:1937566/1‑100 (MQ=255)
                                                                                          tttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga   >  2:1832222/1‑100 (MQ=255)
                                                                                           ttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatg    <  2:2018711/100‑1 (MQ=255)
                                                                                           ttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTgatgat  >  2:438669/1‑100 (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: