Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,121,887 (GC)5→4 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.2% 577.3 / 15.1 132coding (409/531 nt)hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base . (72/55);  total (74/58)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
*NC_0009134,121,8800C.95.5% 568.8 / 17.0 132coding (408/531 nt)hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base . (71/55);  total (74/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

CGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGAT  >  NC_000913/4121780‑4121969
                                                                                                   ||                                                                                         
cGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAg                                                                                            <  1:360394/100‑1 (MQ=255)
cGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAg                                                                                            <  1:3913035/100‑1 (MQ=255)
  aTGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcg                                                                                          >  1:2153918/1‑100 (MQ=255)
   tGGTGTGGAAATGGTTGGTATAGATGCAAATGGCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgc                                                                                         <  1:2546362/100‑1 (MQ=255)
      tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                      >  1:2705688/1‑100 (MQ=255)
      tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                      >  2:3934111/1‑100 (MQ=255)
      tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                      >  2:1927807/1‑100 (MQ=255)
      tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                      >  2:1761104/1‑100 (MQ=255)
      tgtgGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcg                                                                                      >  1:477670/1‑100 (MQ=255)
         ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCATGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                   >  1:2373148/1‑98 (MQ=255)
         ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                   >  2:3045306/1‑98 (MQ=255)
         ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                   >  2:2632107/1‑98 (MQ=255)
         ggAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccg                                                                                   >  1:960528/1‑98 (MQ=255)
          gAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgc                                                                                  >  1:1594498/1‑97 (MQ=255)
          gAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAgcgcgcgccgc                                                                                  <  1:2672711/100‑4 (MQ=255)
                gttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTg                                                                          >  1:2237680/1‑100 (MQ=255)
                gttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTg                                                                          <  2:2132888/100‑1 (MQ=255)
                gttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTg                                                                          <  1:34548/100‑1 (MQ=255)
                 ttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTgg                                                                         <  1:2146818/100‑1 (MQ=255)
                 ttggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTgg                                                                         <  2:1992856/100‑1 (MQ=255)
                  tggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTggg                                                                        <  1:639316/100‑1 (MQ=255)
                  tggtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTggg                                                                        >  2:2871327/1‑100 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTgggtg                                                                      <  1:3410630/100‑3 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                      >  1:1774705/1‑100 (MQ=255)
                    gtATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCg                                                                      >  2:362131/1‑100 (MQ=255)
                     tataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTgggtgt                                                                     >  1:2361431/1‑97 (MQ=255)
                     tataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGt                                                                     >  2:861184/1‑100 (MQ=255)
                     tataGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGt                                                                     >  2:755911/1‑100 (MQ=255)
                         gATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAgg                                                                 <  1:2400579/100‑1 (MQ=255)
                          aTGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAggg                                                                >  2:1418665/1‑100 (MQ=255)
                          aTGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAggg                                                                <  2:2935661/100‑1 (MQ=255)
                           tGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgc                                                               >  1:535825/1‑100 (MQ=255)
                           tGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgc                                                               >  2:4008609/1‑100 (MQ=255)
                            gCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgcc                                                              <  1:3229589/100‑1 (MQ=255)
                               aaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccgcn                                                           <  1:4066624/100‑2 (MQ=255)
                               aaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGgccgcc                                                           >  1:864656/1‑100 (MQ=255)
                                aTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGccgccg                                                          <  1:1378251/100‑1 (MQ=255)
                                    cGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGc                                                      >  1:2310238/1‑100 (MQ=255)
                                    cGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTCCAGTGTTTTCTCACA‑‑GCACGCGCCGCAGCCTGGGCGGAAGGGCCGCCGGAGc                                                      >  2:36899/1‑100 (MQ=255)
                                     gCCTGCAATAGCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCATGGGCGTAAGGGCCGCCgg                                                         <  2:1810454/96‑1 (MQ=255)
                                       cTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGa                                                   >  1:108702/1‑100 (MQ=255)
                                       cTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGa                                                   >  2:4106583/1‑100 (MQ=255)
                                        tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGAt                                                  <  1:4021239/100‑1 (MQ=255)
                                         gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                 >  2:3366195/1‑100 (MQ=255)
                                         gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                 >  1:3122290/1‑100 (MQ=255)
                                         gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGCTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATa                                                 >  2:2277907/1‑100 (MQ=255)
                                          cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAg                                                <  2:535821/100‑1 (MQ=255)
                                          cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAACGGCCGCCGGAGCCGATAg                                                >  2:3313002/1‑100 (MQ=255)
                                           aaTATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGc                                               <  2:3617888/100‑1 (MQ=255)
                                           aaTATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGc                                               <  2:1069285/100‑1 (MQ=255)
                                           aaTATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGc                                               >  1:2886306/1‑100 (MQ=255)
                                           aaTATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGc                                               <  1:659631/100‑1 (MQ=255)
                                            atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCa                                              >  1:3980253/1‑100 (MQ=255)
                                            atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCa                                              >  1:128497/1‑100 (MQ=255)
                                             tatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCaa                                             >  1:3533189/1‑100 (MQ=255)
                                              atCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAAt                                            >  2:3482674/1‑100 (MQ=255)
                                              atCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAAt                                            <  2:1209517/100‑1 (MQ=255)
                                                ccAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATaa                                          <  1:1144354/100‑1 (MQ=255)
                                                ccAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATaa                                          >  2:1778590/1‑100 (MQ=255)
                                                  aaCGCCTTTTCAGCAATTTCACGGGCGCTATGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGa                                        >  2:651933/1‑100 (MQ=255)
                                                  aaCGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGa                                        >  2:387390/1‑100 (MQ=255)
                                                   aCGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGCTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGAt                                       <  2:2280854/100‑1 (MQ=255)
                                                      ccTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGt                                    <  2:1593552/100‑1 (MQ=255)
                                                       cTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGTCCGCCGGAGCCGATAGCAATAAGATCGtt                                   >  2:1018631/1‑100 (MQ=255)
                                                       cTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGtt                                   >  1:2611900/1‑100 (MQ=255)
                                                       cTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGtt                                   >  2:4102889/1‑100 (MQ=255)
                                                       cTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGtt                                   >  1:1352497/1‑100 (MQ=255)
                                                            cAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTg                              >  1:1632920/1‑100 (MQ=255)
                                                             aGCAATTTCACGGGCGCTAAGTTCAGCGTTTTCTAACA‑‑GCGCACGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTgg                             >  2:1878925/1‑100 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                            >  1:1528945/1‑100 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                            <  1:3934154/100‑1 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGAGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                            >  1:2826818/1‑100 (MQ=255)
                                                              gCAATTTCACGGGCGCTAAGTTCAGTGTTTTATAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGc                            <  2:2825787/100‑1 (MQ=255)
                                                               cAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCGGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCt                           >  1:3545157/1‑100 (MQ=255)
                                                                aaTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTg                          <  2:3222189/100‑1 (MQ=255)
                                                                aaTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTg                          <  1:362135/100‑1 (MQ=255)
                                                                aaTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTg                          >  1:1422762/1‑100 (MQ=255)
                                                                 aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCt                           >  1:892122/1‑100 (MQ=255)
                                                                     cACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacca                     <  1:1927821/100‑1 (MQ=255)
                                                                     cACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacca                     >  1:1093662/1‑100 (MQ=255)
                                                                     cACGGGCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGcacca                     <  2:108702/100‑1 (MQ=255)
                                                                         ggCGCTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTc                 <  1:500973/100‑1 (MQ=255)
                                                                          gcgcTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCa                <  1:332066/100‑1 (MQ=255)
                                                                          gcgcTAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCa                >  1:1291559/1‑100 (MQ=255)
                                                                              tAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGt            >  1:907602/1‑100 (MQ=255)
                                                                              tAAGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGt            <  1:4102932/100‑1 (MQ=255)
                                                                               aaGTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAGGGGCCGCCGGAGCCGATCGCAGTAAGATCGTTTTCTGGCTGCACCACGTCACCGtt           >  2:2044743/1‑100 (MQ=255)
                                                                                 gTTCAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAc         >  1:544633/1‑100 (MQ=255)
                                                                                    cAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCgg      <  2:1215366/100‑1 (MQ=255)
                                                                                    cAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCgg      >  2:2366290/1‑100 (MQ=255)
                                                                                    cAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCgg      >  2:2375952/1‑100 (MQ=255)
                                                                                    cAGTGTTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCg       >  2:1045201/1‑99 (MQ=255)
                                                                                      gtgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg    >  2:2124273/1‑100 (MQ=255)
                                                                                      gtgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg    <  1:3077676/100‑1 (MQ=255)
                                                                                      gtgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg    >  1:1328389/1‑100 (MQ=255)
                                                                                      gtgtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg    <  1:1089426/100‑1 (MQ=255)
                                                                                        gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat  <  1:2306722/100‑1 (MQ=255)
                                                                                        gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat  <  1:1448344/100‑1 (MQ=255)
                                                                                        gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat  <  2:1746095/100‑1 (MQ=255)
                                                                                        gtTTTCTAACA‑‑GCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat  <  2:265431/100‑1 (MQ=255)
                                                                                                   ||                                                                                         
CGATGGTGTGGAAATGGTTGGTATAGATGCAAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGAT  >  NC_000913/4121780‑4121969

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