Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,117,903 Δ1 bp 5.9% coding (502/1281 nt) wcaK ← colanic acid biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,117,9030G. 5.9% 223.0 / 7.7 68coding (502/1281 nt)wcaKcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base . (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

CTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCG  >  NC_000913/2117756‑2118048
                                                                                                                                                   |                                                                                                                                                 
cttctGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTg                                                                                                                                                   >  4:272775/1‑148 (MQ=255)
  tctGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACTTAGTTCGCCAGTTGGt                                                                                                                                                 >  1:225949/1‑148 (MQ=255)
  tctGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGt                                                                                                                                                 >  6:221896/1‑148 (MQ=255)
   ctGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGt                                                                                                                                                 >  8:233188/1‑147 (MQ=255)
      tgtgGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTaa                                                                                                                                              <  3:117019/147‑1 (MQ=255)
      tgtgGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTaa                                                                                                                                              >  4:117019/1‑147 (MQ=255)
      tgtgGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACAGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTaa                                                                                                                                              <  6:162463/147‑1 (MQ=255)
           tgATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCt                                                                                                                                        <  2:225949/148‑1 (MQ=255)
              tCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCAt                                                                                                                                     <  8:32007/148‑1 (MQ=255)
                gACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATcc                                                                                                                                   <  2:354538/148‑1 (MQ=255)
                 aCCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCt                                                                                                                                  <  6:192019/148‑1 (MQ=255)
                  ccAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTg                                                                                                                                          <  3:94130/139‑1 (MQ=255)
                  ccAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCACCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTGGGTTAAATGGCTCATCCTg                                                                                                                                 >  4:390122/1‑148 (MQ=255)
                  ccAGCAACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTg                                                                                                                                          >  4:94130/1‑139 (MQ=255)
                          gcgGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGc                                                                                                                         >  6:3284/1‑148 (MQ=255)
                                 cGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACt                                                                                                                  >  2:153145/1‑148 (MQ=255)
                                 cGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACt                                                                                                                  >  5:410712/1‑148 (MQ=255)
                                         tGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCaa                                                                                                          >  6:329757/1‑148 (MQ=255)
                                          gTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGTATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATcc                                                                                                                                   <  5:70259/122‑1 (MQ=255)
                                          gTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGTATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATcc                                                                                                                                   >  6:70259/1‑122 (MQ=255)
                                            tCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAac                                                                                                                     >  8:209933/1‑134 (MQ=255)
                                            tCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAac                                                                                                                     <  7:209933/134‑1 (MQ=255)
                                            tCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCa                                                                                                       >  2:218238/1‑148 (MQ=255)
                                            tCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATAGCTCATCCTGGAACGGGCCAACACTGAGCCCAATCa                                                                                                       >  2:251557/1‑148 (MQ=255)
                                               aCTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATaa                                                                                                    <  7:233188/148‑1 (MQ=255)
                                                            tgATATTGCTGCGTTTCCTCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCttttt                                                                                       <  2:74827/148‑1 (MQ=255)
                                                            tgATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCttttt                                                                                       <  3:86669/148‑1 (MQ=255)
                                                                 ttGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGcc                                                                                   <  5:3284/147‑1 (MQ=255)
                                                                 ttGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCa                                                                                  <  6:371965/148‑1 (MQ=255)
                                                                    ctgcGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCATCACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATaa                                                                               <  1:251557/148‑1 (MQ=255)
                                                                      gcGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATcc                                                                                                                                   <  5:9677/94‑1 (MQ=255)
                                                                      gcGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATcc                                                                                                                                   >  6:9677/1‑94 (MQ=255)
                                                                       cGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACg                                                                            <  6:410712/148‑1 (MQ=255)
                                                                             atcaGAGCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCaa                                                                       <  5:221896/147‑1 (MQ=255)
                                                                              tcaGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACg                                                                            <  6:35671/141‑1 (MQ=255)
                                                                              tcaGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACg                                                                            >  5:35671/1‑141 (MQ=255)
                                                                                aGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTg                                                                   <  3:272775/148‑1 (MQ=255)
                                                                                aGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTg                                                                   >  1:213274/1‑148 (MQ=255)
                                                                                 gATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCGTTTTCGCCATAAACGTGCAAAGTGc                                                                  <  2:393257/148‑1 (MQ=255)
                                                                                  aTCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGa                                                                                                                               <  8:102044/86‑1 (MQ=255)
                                                                                  aTCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGa                                                                                                                               >  7:102044/1‑86 (MQ=255)
                                                                                      aaGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGtt                                                             <  1:153145/148‑1 (MQ=255)
                                                                                            aCCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACt                                                       >  1:422188/1‑148 (MQ=255)
                                                                                            aCCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACt                                                       >  3:384973/1‑148 (MQ=255)
                                                                                              cGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGc                                                     <  3:153482/148‑1 (MQ=255)
                                                                                                          gATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCa                                            >  6:99608/1‑145 (MQ=255)
                                                                                                           aTCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATaga                                        <  5:329757/148‑1 (MQ=255)
                                                                                                                         gTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAAc                          >  1:412091/1‑148 (MQ=255)
                                                                                                                           gACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTg                                                                          <  7:190511/98‑1 (MQ=255)
                                                                                                                           gACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTg                                                                          >  8:190511/1‑98 (MQ=255)
                                                                                                                           gACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACg                         >  4:181804/1‑147 (MQ=255)
                                                                                                                            agcAAAAACGTAGTTCGCCAGTTTGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGAt                       <  3:259426/146‑1 (MQ=255)
                                                                                                                            aCCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGAt                       <  5:271494/148‑1 (MQ=255)
                                                                                                                             ccAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATc                      >  4:392641/1‑148 (MQ=255)
                                                                                                                             ccAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATc                      >  5:369761/1‑148 (MQ=255)
                                                                                                                                aaaaCGTAGTTCGCCAGTTGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACaaaaa                            <  2:360109/138‑1 (MQ=255)
                                                                                                                                aaaaCGTAGTTCGCCAGTTGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACaaaaa                            <  2:354106/138‑1 (MQ=255)
                                                                                                                                aaaaCGTAGTTCGCCAGTTGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACaaaaa                            >  1:360109/1‑138 (MQ=255)
                                                                                                                                aaaaCGTAGTTCGCCAGTTGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACaaaaa                            >  1:354106/1‑138 (MQ=255)
                                                                                                                                 aaaCGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATccgcc                  >  1:15858/1‑148 (MQ=255)
                                                                                                                                  aaCGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATccgccg                 <  4:384973/148‑1 (MQ=255)
                                                                                                                                      tAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCt             <  2:213274/148‑1 (MQ=255)
                                                                                                                                          tCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGa          >  4:102305/1‑147 (MQ=255)
                                                                                                                                          tCGCCACTTCGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGa          >  1:343382/1‑147 (MQ=255)
                                                                                                                                            gCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATa        <  5:99608/147‑1 (MQ=255)
                                                                                                                                             ccAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGAt                                     <  4:277918/117‑1 (MQ=255)
                                                                                                                                             ccAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGAt                                     >  3:277918/1‑117 (MQ=255)
                                                                                                                                             ccAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAAt      <  3:390122/148‑1 (MQ=255)
                                                                                                                                                 ttGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCAAGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGcg  >  2:326915/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAGGATCAGCGCGTCGCAGTGACCAAAAACGTAGTTCGCCAGTTGGTTAAATTGCTCATCCTGGAACGGGCCAACACTGTGACCAATCATAAACAGCGGCTTTTTCGCCATAAACGTGCAAAGTGCATGTTCAAACTGCGGCACGCCATAGAGATCGACAAAAAACGATCCGCCGACCTGGATAATGGCG  >  NC_000913/2117756‑2118048

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