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.96.4% 605.3 / 13.3 137coding (409/531 nt)hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base . (67/65);  total (70/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01
*NC_0009134,121,8800C.96.4% 602.8 / 13.7 137coding (408/531 nt)hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base . (67/65);  total (70/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 167 total aligned reads displayed.

GATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGAT  >  NC_000913/4121781‑4121972
                                                                                                  ||                                                                                            
gATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGAGCGCTAAGTTCAGTGTTTTCTAACAgc                                                                                              >  1:2879913/1‑100 (MQ=255)
    gtgtGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgc                                                                                          <  1:60250/100‑1 (MQ=255)
     tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgct                                                                                         <  2:3139237/100‑2 (MQ=255)
     tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                         <  2:2195010/100‑1 (MQ=255)
     tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                         >  1:3895537/1‑100 (MQ=255)
     tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                         <  2:1486079/100‑1 (MQ=255)
     tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgc                                                                                          >  2:489970/1‑99 (MQ=255)
      gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgc                                                                                        >  1:2684849/1‑100 (MQ=255)
      gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgc                                                                                        >  1:2601415/1‑100 (MQ=255)
      gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgc                                                                                        >  1:2387991/1‑100 (MQ=255)
      gtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGATAAGTTCAGTGTTTTCTAACAgcgcgcgc                                                                                        >  1:3618007/1‑100 (MQ=255)
       tgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgcc                                                                                       <  1:2022390/100‑2 (MQ=255)
       tgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgcc                                                                                       <  1:2258099/100‑2 (MQ=255)
        ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                      >  1:1839593/1‑98 (MQ=255)
        ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                      >  2:3036386/1‑98 (MQ=255)
        ggAAATGGTTGGTATAGATACAAATGTAGCATGCAATATACAACACCTTATCAGCAATTTAACAAACGCTAAATTCAGTGTTTTATAACAtcgatcgcc                                                                                       >  2:1522048/1‑93 (MQ=255)
         gAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgn                                                                                     <  1:1436173/100‑4 (MQ=255)
          aaaTGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgca                                                                                    <  1:1529927/100‑5 (MQ=255)
           aaTGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgcag                                                                                   <  1:2740047/100‑6 (MQ=255)
             tggttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgcagcc                                                                                 <  1:1347755/100‑8 (MQ=255)
                ttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTg                                                                             <  2:650681/99‑1 (MQ=255)
                 tggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTggg                                                                           <  2:68681/100‑1 (MQ=255)
                  ggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGc                                                                          >  2:881100/1‑100 (MQ=255)
                  ggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGc                                                                          <  2:917020/100‑1 (MQ=255)
                  ggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTCGGc                                                                          >  2:3382971/1‑100 (MQ=255)
                   gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                         <  2:255901/100‑1 (MQ=255)
                   gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                         >  2:2828021/1‑100 (MQ=255)
                      taGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTaa                                                                      >  2:411919/1‑100 (MQ=255)
                       aGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAg                                                                     >  1:1070842/1‑100 (MQ=255)
                        gATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAgg                                                                    <  2:3024303/100‑1 (MQ=255)
                          tGCAAATGTCGCCTGCAATATCCAAGGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgc                                                                  >  2:1027106/1‑100 (MQ=255)
                          tGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgc                                                                  >  1:3385377/1‑100 (MQ=255)
                           gCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGATAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgcc                                                                 >  1:29759/1‑100 (MQ=255)
                            cAAATGTCGCCTGCAATATCCAACGCCTGTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccg                                                                <  2:1591358/100‑1 (MQ=255)
                              aaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccgcc                                                              >  1:3688988/1‑100 (MQ=255)
                              aaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccgcc                                                              >  1:1439814/1‑100 (MQ=255)
                               aTGTCGCCTGCAATATCCAACGCCTATTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGccgccg                                                             >  2:1777152/1‑100 (MQ=255)
                                  tCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAg                                                          <  2:1137713/100‑1 (MQ=255)
                                   cGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTCAGGGCCGCCGGAGc                                                         >  2:463949/1‑100 (MQ=255)
                                    gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        >  1:1142499/1‑100 (MQ=255)
                                    gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGcc                                                        <  1:1381801/100‑1 (MQ=255)
                                     ccTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCg                                                       <  1:1595831/100‑1 (MQ=255)
                                       tGCAATATCCAACGCCTTTTCAGCAATTTCACGTGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGa                                                      <  2:204784/99‑1 (MQ=255)
                                       tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     <  1:1061506/100‑1 (MQ=255)
                                       tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     <  1:3673072/100‑1 (MQ=255)
                                       tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                     >  1:2025883/1‑100 (MQ=255)
                                        gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                    <  1:766214/100‑1 (MQ=255)
                                        gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                    >  2:1246341/1‑100 (MQ=255)
                                         cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAg                                                   <  1:3290824/100‑1 (MQ=255)
                                           atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCACGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCa                                                 >  1:1802019/1‑100 (MQ=255)
                                             atCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAAt                                               >  1:1206246/1‑100 (MQ=255)
                                              tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATa                                              >  1:173769/1‑100 (MQ=255)
                                              tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATa                                              <  2:2684825/100‑1 (MQ=255)
                                              tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATa                                              >  1:2314015/1‑100 (MQ=255)
                                                cAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCACGCGCCGCAACCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATaa                                             >  1:214961/1‑99 (MQ=255)
                                                 aaCGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGa                                           >  2:3089296/1‑100 (MQ=255)
                                                  aCGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGAt                                          <  1:1671946/100‑1 (MQ=255)
                                                    gCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCg                                        <  2:3851923/100‑1 (MQ=255)
                                                     ccTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCg                                        >  1:3692660/1‑99 (MQ=255)
                                                      cTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGt                                       <  1:881887/99‑1 (MQ=255)
                                                       ttttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGttt                                     <  2:2109002/100‑1 (MQ=255)
                                                        tttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGtttt                                    >  1:3586545/1‑100 (MQ=255)
                                                         ttCAGCAATTTCACGGGCGAAAAGTTCAGTGTTGTTTAATA‑‑GCGAGCGCCGCAGCCTGGGCGAAATGACAACAGGAGCGGATATCAATTAGagaggttt                                    >  2:3576245/1‑92 (MQ=255)
                                                          tCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCt                                  >  1:2892579/1‑100 (MQ=255)
                                                          ccgatcttttCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCt                                  <  1:381058/93‑1 (MQ=255)
                                                           cAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                                 <  2:1752343/100‑1 (MQ=255)
                                                            aGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGACTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTgg                                >  2:3786816/1‑100 (MQ=255)
                                                             gCAATTTCACGTGCGCTATGTTCATTTTTTTCTAACA‑‑GCGTGTGCCGCAGGCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               <  1:720254/100‑1 (MQ=255)
                                                             gCAATTTCACGGGCGCTACGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTACGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                                 <  2:1066252/100‑1 (MQ=255)
                                                             gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               >  1:1024662/1‑100 (MQ=255)
                                                             gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               <  2:230172/100‑1 (MQ=255)
                                                             gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAAC‑‑CGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                                 >  1:547848/1‑98 (MQ=255)
                                                               aaTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                               >  1:3541752/1‑100 (MQ=255)
                                                                aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGc                            >  1:1519431/1‑100 (MQ=255)
                                                                 tttCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCtagatcgg                           >  2:381052/1‑93 (MQ=255)
                                                                   tCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacc                         >  2:3486042/1‑100 (MQ=255)
                                                                   tCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacc                         <  2:507727/100‑1 (MQ=255)
                                                                     aCGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcaccac                       >  1:1028393/1‑100 (MQ=255)
                                                                     aCGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcaccac                       >  1:225716/1‑100 (MQ=255)
                                                                      cGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcaccac                       >  2:122887/1‑99 (MQ=255)
                                                                      cGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACg                      <  1:3486081/100‑1 (MQ=255)
                                                                      cGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACg                      <  1:393546/100‑1 (MQ=255)
                                                                        ggCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCTGCAGCCTGTGTGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTc                    <  1:131670/100‑1 (MQ=255)
                                                                           gcTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCAcc                 <  1:3911918/100‑1 (MQ=255)
                                                                           gcTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCAcc                 <  2:3064960/100‑1 (MQ=255)
                                                                           gcTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCAcc                 <  2:2871871/100‑1 (MQ=255)
                                                                            cTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCAGCAGCCTGGGCGTAAGGGCCGCCGTAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACg                      >  2:2493049/1‑94 (MQ=255)
                                                                               aGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTa             >  2:1346899/1‑100 (MQ=255)
                                                                               aGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTa             >  1:3316114/1‑100 (MQ=255)
                                                                                gTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAc            >  1:1581251/1‑100 (MQ=255)
                                                                                  tCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCg          <  1:2420072/100‑1 (MQ=255)
                                                                                     gtgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg       >  1:503901/1‑100 (MQ=255)
                                                                                     gtgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg       <  1:2518185/100‑1 (MQ=255)
                                                                                     gtgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGATTTCTGGCTGCACCACGTCACCGTTACCGGtg       >  1:1610387/1‑100 (MQ=255)
                                                                                       gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat     <  2:3586506/100‑1 (MQ=255)
                                                                                       gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat     >  2:260293/1‑100 (MQ=255)
                                                                                       gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat     >  1:3617950/1‑100 (MQ=255)
                                                                                        ttttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatg    >  2:1717541/1‑100 (MQ=255)
                                                                                          ttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacgtg  <  2:2459351/100‑5 (MQ=255)
                                                                                          ttCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacgtg  <  1:1448270/100‑5 (MQ=255)
                                                                                                  ||                                                                                            
GATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGAT  >  NC_000913/4121781‑4121972

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: