Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,618,315 G→A 6.4% Q113* (CAA→TAA)  wzy ← O‑antigen polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,618,3150GA 6.4% 185.9 / 2.0 62Q113* (CAA→TAA) wzyO‑antigen polymerase
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base A (2/2);  total (32/30)
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.66e-01

CCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGCCAACAGCAGCAGCA  >  NC_002947/5618167‑5618446
                                                                                                                                                    |                                                                                                                                   
cccAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTg                                                                                                                                     <  3:352074/149‑1 (MQ=255)
cccAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTg                                                                                                                                     <  4:326098/149‑1 (MQ=255)
 ccAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTgg                                                                                                                                    >  3:318586/1‑149 (MQ=255)
         gggTTGGATATCTCTCCGATACCCGCCAGCCTTAAGAGCATAGGCGCACGCTGCAACGCATAGAACTTGATGATGGAGACCAACGCACCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATATCTGTTGATTGCGGGCTTGGCGAAGCTg                                                                                                                            <  4:105554/149‑1 (MQ=255)
           gTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTggg                                                                                                                          >  6:78700/1‑149 (MQ=255)
            tGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGc                                                                                                                         <  7:276767/149‑1 (MQ=255)
                atatCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCCAGCTTGGCCaag                                                                                                                     <  3:105061/149‑2 (MQ=255)
                  atCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAg                                                                                                                   >  7:360735/1‑149 (MQ=255)
                     tcGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCg                                                                                                                 >  8:43232/1‑148 (MQ=255)
                         cGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAAc                                                                                                            >  3:298851/1‑149 (MQ=255)
                           aTACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGc                                                                                                          <  7:43232/149‑1 (MQ=255)
                               ccGCCAGCCTGACCAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCcagc                                                                                                      <  3:257396/149‑1 (MQ=255)
                                  ccAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGCAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCTGCAATAACTGACGAATGCGGGCTTGGCCAAACTGGGCCAACAGCGCAAACGCccgcag                                                                                                    >  6:122950/1‑148 (MQ=255)
                                       cTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTggg                                                                                                                          <  4:135460/121‑1 (MQ=255)
                                       cTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTggg                                                                                                                          >  3:135460/1‑121 (MQ=255)
                                        tGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAacac                                                                                              >  2:105767/1‑148 (MQ=255)
                                        tGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTAGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAacac                                                                                              >  3:168246/1‑148 (MQ=255)
                                        tGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTAGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAacac                                                                                              >  3:162618/1‑148 (MQ=255)
                                         gAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAacac                                                                                              >  2:354171/1‑147 (MQ=255)
                                           aCAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAgg                                                                                          <  5:206490/149‑1 (MQ=255)
                                           aCAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAgg                                                                                          <  8:360735/149‑1 (MQ=255)
                                              gCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATg                                                                                       >  3:177202/1‑149 (MQ=255)
                                              gCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATg                                                                                       >  3:183046/1‑149 (MQ=255)
                                               cAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATg                                                                                       >  6:332672/1‑148 (MQ=255)
                                               cAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATg                                                                                       >  5:345930/1‑148 (MQ=255)
                                                        aCGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAgc                                                                             <  4:298851/149‑1 (MQ=255)
                                                          gCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAgcgc                                                                           >  2:58634/1‑149 (MQ=255)
                                                              cACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGa                                                                       <  6:208557/149‑1 (MQ=255)
                                                                   cATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTcg                                                                  >  5:150162/1‑149 (MQ=255)
                                                                   cATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTcg                                                                  >  7:320872/1‑149 (MQ=255)
                                                                         aCTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATg                                                            >  3:323804/1‑149 (MQ=255)
                                                                         aCTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATg                                                            <  2:285840/149‑1 (MQ=255)
                                                                                  tgGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTcg                                                   <  1:354171/149‑1 (MQ=255)
                                                                                     agaCCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCtt                                                                         >  4:125322/1‑124 (MQ=255)
                                                                                     agaCCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCtt                                                                         <  3:125322/124‑1 (MQ=255)
                                                                                      gaacAACCCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGACTGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTcggcgg                                               <  4:270537/146‑1 (MQ=255)
                                                                                          aaCGCAGCGATGGCCAGCCGCGCGCTGCCAAGCAGCAATAACAGACGAATGCGGGCTAGGCCAAGCTGGGCCAACAGCGGACACGCCAGCAGGCCCACCAGGATGTAGAGCAACCCCTAGAGCTCGCGCATCGGCTCCTCGGCGGGCGa                                           >  5:286833/1‑149 (MQ=255)
                                                                                           aCGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGAc                                          >  7:19863/1‑149 (MQ=255)
                                                                                             gCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTAGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCa                                        <  4:168246/149‑1 (MQ=255)
                                                                                             gCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTAGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCa                                        <  4:162618/149‑1 (MQ=255)
                                                                                                  gATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCa                                   <  4:323804/149‑1 (MQ=255)
                                                                                                       ccAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCt                               >  7:361603/1‑148 (MQ=255)
                                                                                                        cagcagCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCAt                                       >  6:39172/1‑139 (MQ=255)
                                                                                                        cagcagCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCAt                                       <  5:39172/139‑1 (MQ=255)
                                                                                                        cagcagCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCa                             >  4:241769/1‑149 (MQ=255)
                                                                                                         agcagcGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAg                            >  5:52467/1‑149 (MQ=255)
                                                                                                          gcagcGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGc                           <  3:241769/149‑1 (MQ=255)
                                                                                                                 gcTGCCAAGCAGCAATAACTGAGGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGTAACACCAGGATGTAGAGCAAGCGGTTGAGCTCGCGCATGGGctcct                                                     <  7:293887/116‑1 (MQ=255)
                                                                                                                 gcTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCAAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAACCGCTTGAGCTCGCGCATGCGctcct                                                     >  8:293887/1‑116 (MQ=255)
                                                                                                                     ccAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAg                  <  5:78700/147‑1 (MQ=255)
                                                                                                                     ccAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGc                <  6:150162/149‑1 (MQ=255)
                                                                                                                        agcagcAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGCCaa             <  6:52467/149‑1 (MQ=255)
                                                                                                                          cagcaATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAGGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGCCAAca           <  4:183046/149‑1 (MQ=255)
                                                                                                                          cagcaATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAGGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGCCAAca           <  4:177202/149‑1 (MQ=255)
                                                                                                                            gcaATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGCCAAcagc         <  6:345930/149‑1 (MQ=255)
                                                                                                                                 aaCTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGCCAAcagcagcag    >  4:376797/1‑149 (MQ=255)
                                                                                                                                  aCTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGCCAAcagcagcag    >  5:107063/1‑148 (MQ=255)
                                                                                                                                  aCTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGCCAAcagcagcag    <  5:332672/148‑1 (MQ=255)
                                                                                                                                   cTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGc                                                                     >  7:77310/1‑82 (MQ=255)
                                                                                                                                   cTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGc                                                                     <  8:77310/82‑1 (MQ=255)
                                                                                                                                   cTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTcggcgg                                               <  1:172543/104‑1 (MQ=255)
                                                                                                                                   cTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTcggcgg                                               >  2:172543/1‑104 (MQ=255)
                                                                                                                                    tGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGCCAACAgcagcagca  >  8:142307/1‑148 (MQ=255)
                                                                                                                                        gAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAg                            <  2:83052/118‑1 (MQ=255)
                                                                                                                                        gAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAg                            >  1:83052/1‑118 (MQ=255)
                                                                                                                                                    |                                                                                                                                   
CCCAGGATGGGGTGGGATATCTCGCCGATACCCGCCAGCCTGAACAGCAAAGGCGCACGCTGCACCCCATAGAACTTGATGATGGAGACCAACGCAGCGATGGCCAGCAGCGCGCTGCCAAGCAGCAATAACTGACGAATGCGGGCTTGGCCAAGCTGGGCCAACAGCGGAAACGCCAGCAGGAACACCAGGATGTAGAGCAAGCGCTTGAGCTCGCGCATGGGCTCCTCGGCGGGCGACCATGCCAGGCTCAGCCCGCTCCAGGCCAACAGCAGCAGCA  >  NC_002947/5618167‑5618446

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