Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,608,294 T→A 11.7% intergenic (‑156/+32) cyaB ← / ← waaA cyclolysin secretion/processing ATP‑binding protein CyaB/3‑deoxy‑D‑manno‑octulosonic acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,608,2940TA11.7% 151.4 / 11.4 60intergenic (‑156/+32)cyaB/waaAcyclolysin secretion/processing ATP‑binding protein CyaB/3‑deoxy‑D‑manno‑octulosonic acid transferase
Reads supporting (aligned to +/- strand):  ref base T (25/28);  new base A (4/3);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

ACGTCCCACAAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAG  >  NC_002947/5608147‑5608442
                                                                                                                                                    |                                                                                                                                                     
aCGTCCCACAAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCActct                                                                                                                                                       >  8:191783/1‑148 (MQ=255)
 agtgtcacaCGCTTGCGTCGGAAGCTTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTa                                                                                                                                                     <  7:309809/143‑1 (MQ=255)
      cacaAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCgg                                                                                                                                               <  8:74059/149‑1 (MQ=255)
       acaAGCTTCCGACGCAAGCGTCTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACAAGAATCAGACATAACCACTCTAGACCgg                                                                                                                                               >  4:80585/1‑148 (MQ=255)
          aGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCt                                                                                                                                           <  2:38580/149‑1 (MQ=255)
            cTTCCGACTCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCGAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATATCCAGAATCAGAGATATCCACTCTAGACCGGGTCttt                                                                                                                                         <  8:48313/149‑1 (MQ=255)
            cTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGGCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCttt                                                                                                                                         <  6:306047/149‑1 (MQ=255)
             ttCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCttt                                                                                                                                         >  1:166727/1‑148 (MQ=255)
              tCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCttt                                                                                                                                         >  2:190997/1‑147 (MQ=255)
               ccGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTg                                                                                                                                        <  3:177561/147‑1 (MQ=255)
                cGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAg                                                                                                                                      <  1:190997/148‑1 (MQ=255)
                cGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGt                                                                                                                                     <  2:166727/149‑1 (MQ=255)
                cGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGt                                                                                                                                     >  6:311488/1‑149 (MQ=255)
                cGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGACGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGt                                                                                                                                     >  7:156262/1‑149 (MQ=255)
                      aaGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCAAAGACCCGGTCTAGAGTGGatat                                                                                                                               <  7:84444/149‑1 (MQ=39)
                      aaGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGt                                                                                                                                             >  7:201263/1‑135 (MQ=255)
                      aaGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGt                                                                                                                                             <  8:201263/135‑1 (MQ=255)
                      aaGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGt                                                                                                                                     <  3:23454/143‑1 (MQ=255)
                      aaGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGt                                                                                                                                     >  4:23454/1‑143 (MQ=255)
                          gtttgACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTTAGGGCGAGAAATGCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCg                                                                                                                           <  6:278385/146‑1 (MQ=255)
                                           tttAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAg                                                                                                          >  1:398618/1‑149 (MQ=255)
                                             tAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAgg                                                                                                         <  7:191783/148‑1 (MQ=255)
                                               aaGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCAAAGACCCGGTCTAGAGTGGATATGTctga                                                                                                                         <  4:89807/130‑4 (MQ=37)
                                               aaGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCAAAGACCCGGTCTAGAGTGGATATGTctga                                                                                                                         >  3:89807/1‑127 (MQ=37)
                                                         gTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTc                                                                                            <  8:156262/149‑1 (MQ=255)
                                                          tAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTc                                                                                            >  2:58672/1‑148 (MQ=255)
                                                            ggAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGa                                                                                          <  1:58672/148‑1 (MQ=255)
                                                            ggAAGTTTCCTAACGGAGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAAGCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGTATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGa                                                                                          <  3:80585/148‑1 (MQ=255)
                                                                       aaCGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGc                                                                               >  1:349284/1‑148 (MQ=255)
                                                                               gaAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGccccc                                                                      <  2:398618/149‑1 (MQ=255)
                                                                                   tCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGAt                                                                  <  5:311488/149‑1 (MQ=255)
                                                                                   tCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGAt                                                                  >  1:72456/1‑149 (MQ=255)
                                                                                      aGTGTAGGAAGATACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTgg                                                               <  5:150872/149‑1 (MQ=255)
                                                                                          tAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCAACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCACCAGACCCTCGAGCAATCGCTGCAAAGCCACATGATTGGCCCg                                                           >  2:197972/1‑149 (MQ=255)
                                                                                             gAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGAAGATAGCCCCCAGACCATAGAGCAATCGCTGCAAAGCCCCCGGATTGGCCCGCAt                                                        >  1:338106/1‑149 (MQ=255)
                                                                                              aaGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATc                                                       <  7:287777/149‑1 (MQ=255)
                                                                                                 ttACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGTGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGa                                                                                                 <  8:240297/104‑1 (MQ=255)
                                                                                                 ttACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGa                                                                                                 >  7:240297/1‑104 (MQ=255)
                                                                                                 ttACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGCTCTTTGAGTGGATATGTCGGTTATCAGCGCATCAACGCCCCCAGa                                                                                                 >  7:246689/1‑104 (MQ=255)
                                                                                                       cTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGCCACTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTg                                              >  5:41458/1‑149 (MQ=255)
                                                                                                        ttACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGcc                                                 <  6:43947/145‑1 (MQ=255)
                                                                                                        ttACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGcc                                                 >  5:43947/1‑145 (MQ=255)
                                                                                                          aCGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGcc                                            <  1:197972/148‑1 (MQ=255)
                                                                                                          aCGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGc                                             <  2:349284/147‑1 (MQ=255)
                                                                                                            gCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGGGGATCAAGGCGCCCAGGCCCTCGAGCAATCGCTGGCAAGCCCCCTGAGAGGCCCGCATCCCGgcccggcccga                                         >  6:224777/1‑148 (MQ=255)
                                                                                                                   gAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCAc                                    >  5:70151/1‑147 (MQ=255)
                                                                                                                   gAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCAc                                    <  6:70151/147‑1 (MQ=255)
                                                                                                                      taaataaCCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATg                               >  4:66368/1‑149 (MQ=255)
                                                                                                                         ataaCCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGa                                                                                          <  6:414767/87‑1 (MQ=255)
                                                                                                                         ataaCCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGa                                                                                          >  5:414767/1‑87 (MQ=255)
                                                                                                                            aCCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTTGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTgg                         <  5:224777/149‑1 (MQ=255)
                                                                                                                               aGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGc                                                                                        <  3:336599/83‑1 (MQ=255)
                                                                                                                               aGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGc                                                                                        >  4:336599/1‑83 (MQ=255)
                                                                                                                                     cgACATATCCACTCAAAGACCCGGTCTAGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGc                >  4:153161/2‑149 (MQ=38)
                                                                                                                                      gACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCg               >  2:412985/1‑149 (MQ=255)
                                                                                                                                      gACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCg               >  4:320993/1‑149 (MQ=255)
                                                                                                                                       aCATATCCACTCAAAGACCCGGTCTAGAGTGGAATTGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCATAGCCCACTGATAGGCCCCCATCACCGCCCTGCCCGCCACACCCAAGCGCTGCGCGTCCAGCgg              >  1:262760/1‑149 (MQ=255)
                                                                                                                                            tCCACTCAAAGACCCGGTCTAGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCg               >  5:98311/1‑143 (MQ=39)
                                                                                                                                            tCCACTCAAAGACCCGGTCTAGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCg               <  6:98311/143‑1 (MQ=39)
                                                                                                                                             ccACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGt          >  4:401038/1‑147 (MQ=255)
                                                                                                                                             ccACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGt          <  3:401038/147‑1 (MQ=255)
                                                                                                                                                 tctAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCgg              <  1:339721/139‑1 (MQ=255)
                                                                                                                                                 tctAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCgg              >  2:339721/1‑139 (MQ=255)
                                                                                                                                                    tAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCCGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAg  >  7:160060/1‑149 (MQ=255)
                                                                                                                                                    tAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTgg                         >  5:348181/1‑126 (MQ=255)
                                                                                                                                                    tAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTgg                         <  6:348181/126‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACGTCCCACAAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAG  >  NC_002947/5608147‑5608442

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