Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,608,300 G→C 6.2% intergenic (‑162/+26) 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,3000GC 6.2% 262.4 / 2.1 83intergenic (‑162/+26)cyaB/waaAcyclolysin secretion/processing ATP‑binding protein CyaB/3‑deoxy‑D‑manno‑octulosonic acid transferase
Reads supporting (aligned to +/- strand):  ref base G (34/43);  new base C (4/2);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACAAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAGGCGGCG  >  NC_002947/5608153‑5608448
                                                                                                                                                    |                                                                                                                                                    
cacaAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCgg                                                                                                                                                     >  1:334659/1‑149 (MQ=255)
cacaAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCgg                                                                                                                                                     >  2:195349/1‑149 (MQ=255)
 acaAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTATATTCTTAGGATGTTTCCTAAAGGCGAGAAATACAGTGTAGGAAGTTACAACTTAAGCCGGAGGAATAAATAACGAGAATAATACATTTACACTCTAGAGCggg                                                                                                                                                    >  6:385694/1‑149 (MQ=255)
      cTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCttt                                                                                                                                               >  5:199538/1‑149 (MQ=255)
      cTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCttt                                                                                                                                               >  6:159942/1‑149 (MQ=255)
         ctgACGCAAGCGTGTGACAATTCTGACTTTTAAAGTATATTCGTACGAAGTTTCATAACGGCGAGAAATCCAGTGTAGGAAGTTACAAATTTCGGCGTAGGTTTAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGa                                                                                                                                             <  3:72821/146‑1 (MQ=255)
         ccGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCttt                                                                                                                                               <  5:86474/146‑1 (MQ=255)
         ccGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGa                                                                                                                                             <  2:263643/148‑1 (MQ=255)
         ccGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTAGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCAAAGACCCGGTCta                                                                                                                                                <  5:380288/146‑2 (MQ=255)
           gACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTg                                                                                                                                          <  5:159942/149‑1 (MQ=255)
           gACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTg                                                                                                                                          <  2:342760/149‑1 (MQ=255)
            aCGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTgg                                                                                                                                         <  8:26436/149‑1 (MQ=255)
             cGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGa                                                                                                                                        <  5:290597/149‑1 (MQ=255)
             cGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGa                                                                                                                                        <  3:148953/149‑1 (MQ=255)
              gCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGat                                                                                                                                       <  3:94085/149‑1 (MQ=255)
               cAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGata                                                                                                                                      <  2:52556/149‑1 (MQ=255)
                aaGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGatat                                                                                                                                     >  4:138612/1‑149 (MQ=255)
                 aGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATg                                                                                                                                    <  2:376033/149‑1 (MQ=255)
                  gCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTg                                                                                                                                          >  3:203712/1‑142 (MQ=255)
                  gCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTg                                                                                                                                          <  4:203712/142‑1 (MQ=255)
                     tgtgACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCgg                                                                                                                                <  6:35942/149‑1 (MQ=255)
                       tgACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACAAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGt                                                                                                                               >  5:312112/1‑148 (MQ=255)
                        gACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTa                                                                                                                             >  6:209462/1‑149 (MQ=255)
                         acacTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTAt                                                                                                                            >  2:228917/1‑149 (MQ=255)
                             ttCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGt                                                                                                                                           <  2:281421/130‑1 (MQ=255)
                             ttCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGt                                                                                                                                           >  1:281421/1‑130 (MQ=255)
                             ttCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGc                                                                                                                        >  3:39856/1‑149 (MQ=255)
                               cTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCg                                                                                                                       >  3:287768/1‑148 (MQ=255)
                                tGACTTTTAAAGCCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGa                                                                                                                     >  5:125746/1‑149 (MQ=255)
                                    ttttAAATTCTATTCGTACGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCaa                                                                                                                 <  6:247153/149‑1 (MQ=255)
                                      ttAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGat                                                                                                                                       >  2:226994/1‑125 (MQ=255)
                                      ttAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAg                                                                                                                <  5:359999/148‑1 (MQ=255)
                                      ttAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAg                                                                                                                <  6:199538/148‑1 (MQ=255)
                                            tCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGccccc                                                                                                          >  7:214828/1‑148 (MQ=255)
                                            tCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCCGGTATTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCGCCa                                                                                                         >  7:80929/1‑149 (MQ=255)
                                             cTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCa                                                                                                         <  2:79569/148‑1 (MQ=255)
                                                    tAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCg                                                                                                 >  1:254941/1‑149 (MQ=255)
                                                      ggAAGTTTCCTAACGTCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGa                                                                                                <  3:807/148‑1 (MQ=255)
                                                      ggAAGTTTCCTAACGGCGATAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGa                                                                                                <  3:138612/148‑1 (MQ=255)
                                                      ggAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGa                                                                                                <  6:125746/148‑1 (MQ=255)
                                                      ggAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGa                                                                                                <  4:287768/148‑1 (MQ=255)
                                                      ggAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAg                                                                                               >  7:275104/1‑149 (MQ=255)
                                                      ggAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAg                                                                                               <  8:214828/149‑1 (MQ=255)
                                                       gAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCACGGCCCCCAGACCCTCGAg                                                                                               >  2:113820/1‑148 (MQ=255)
                                                       gAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGc                                                                                              <  1:195349/149‑1 (MQ=255)
                                                               cTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTg                                                                                      >  5:35062/1‑149 (MQ=255)
                                                                     gCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCaa                                                                                   >  8:46148/1‑146 (MQ=255)
                                                                     gCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCaa                                                                                   <  7:46148/146‑1 (MQ=255)
                                                                        agaAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGcccc                                                                                                           <  2:31891/119‑1 (MQ=255)
                                                                        agaAATCCAGTGTAGGAAGTTACAACTTACCCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGcccc                                                                                                           >  1:31891/1‑119 (MQ=255)
                                                                          aaaTCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCt                                                                           >  6:210877/1‑149 (MQ=255)
                                                                           aaTCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCt                                                                           <  6:312112/148‑1 (MQ=255)
                                                                             tCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGAt                                                                        <  5:209462/149‑1 (MQ=255)
                                                                                 gtgtAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGc                                                                    <  1:228917/149‑1 (MQ=255)
                                                                                   gtAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGccc                                                                  >  8:262255/1‑149 (MQ=255)
                                                                                      ggAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATc                                                                                                                           >  7:347460/1‑89 (MQ=255)
                                                                                      ggAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATc                                                                                                                           <  8:347460/89‑1 (MQ=255)
                                                                                           ttACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCAcc                                                          >  1:153050/1‑149 (MQ=255)
                                                                                           ttACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCAcc                                                          <  2:153050/149‑1 (MQ=255)
                                                                                           ttACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCAcc                                                          >  3:128537/1‑149 (MQ=255)
                                                                                               aaCTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTg                                                                                      >  5:354584/1‑117 (MQ=255)
                                                                                               aaCTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTg                                                                                      <  6:354584/117‑1 (MQ=255)
                                                                                                 cTTAGGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACACGGTCTTTGAGTGGATATGTCGGTTATCGGCGGATCCAGGCCCCCGGACCCGCGCGCCAAAGCGGCCAAGCCCCCGGAGTGGCCGGCATCACCGCACt                                                     >  2:293167/1‑148 (MQ=255)
                                                                                                      gCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTATACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGc                                               <  8:80929/149‑1 (MQ=255)
                                                                                                      gCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTGGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCGCCGACCCGCGAGGCATCGCTGCACAGCCCCGGGAGTGGCCCGCCTCACCGCCCTGCCCGc                                               >  8:174388/1‑149 (MQ=255)
                                                                                                       ccGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGcc                                              <  5:210877/149‑1 (MQ=255)
                                                                                                        gttAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCt                                             <  6:26108/147‑1 (MQ=255)
                                                                                                             gAATAAATAACCAGATTCAGTCATATCGACTGTAGACCGGGGCTTTGTGTGTATATGTCGGTTATCAGCGGATCATGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGcc                                                                               <  8:221112/110‑1 (MQ=255)
                                                                                                             gAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCATCGGATCAAGGCCCCCAGACCATCGAGCAATCGATGCAAAGcc                                                                               >  7:221112/1‑110 (MQ=255)
                                                                                                                  aataaCCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGCCCCTCGGGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCACTgcg                                   >  3:219190/1‑149 (MQ=255)
                                                                                                                  aataaCCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATgcg                                   <  7:262255/149‑1 (MQ=255)
                                                                                                                   ataaCCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCGAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCCCCCATgcgc                                  >  4:259732/1‑149 (MQ=255)
                                                                                                                   ataaCCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTAACCCATgcgc                                  >  1:241663/1‑149 (MQ=255)
                                                                                                                      aacAGAATCAGAAATATCCACTCTAGACCGGTTCTTTGAGTTGATATTTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTg                                <  5:385694/146‑1 (MQ=255)
                                                                                                                       ccAGAATCAGACATATCCACTCTTGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGGCGCCTGATTGGCCCGCCTCAGCGCCCTGCCCGCCTCACCCATGCGCTggg                              <  1:293167/149‑1 (MQ=255)
                                                                                                                        cAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCGCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGc                             >  5:208523/1‑149 (MQ=255)
                                                                                                                        cAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGc                             >  6:72528/1‑149 (MQ=255)
                                                                                                                               cgACATATCCACTCAAAGACCCGGTCTAGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCAGCCTCACCCATGCGCTGGGCGTCCTGc                      >  2:30329/2‑149 (MQ=38)
                                                                                                                                gACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGAccc                                                                                                    <  3:36315/70‑1 (MQ=255)
                                                                                                                                gACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGAccc                                                                                                    >  4:36315/1‑70 (MQ=255)
                                                                                                                                 agatatCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTTGCCCGCATCACCGCCCTGCCCGCCTCa                                           <  5:228223/124‑1 (MQ=255)
                                                                                                                                 aCATATCCACTCAAAGACCCGGTCTAGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCgg                    >  8:213037/1‑149 (MQ=38)
                                                                                                                                         aCTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCCCCCCTGCGCTGGGCGTCCTGCGGCAGTTccc            >  8:105719/1‑148 (MQ=255)
                                                                                                                                               aGACCCGGTCTAGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAggc      <  7:213037/149‑1 (MQ=255)
                                                                                                                                                gACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAg        <  4:128537/147‑1 (MQ=255)
                                                                                                                                                  ccGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAggcgg    <  7:304724/149‑1 (MQ=255)
                                                                                                                                                    gggTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAGgcggcg  <  4:39856/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACAAGCTTCCGACGCAAGCGTGTGACACTTCTGACTTTTAAAGTCTATTCGTAGGAAGTTTCCTAACGGCGAGAAATCCAGTGTAGGAAGTTACAACTTACGCCGTAGGAATAAATAACCAGAATCAGACATATCCACTCTAGACCGGGTCTTTGAGTGGATATGTCGGTTATCAGCGGATCAAGGCCCCCAGACCCTCGAGCAATCGCTGCAAAGCCCCCTGATTGGCCCGCATCACCGCCCTGCCCGCCTCACCCATGCGCTGGGCGTCCTGCGGCAGTTCCACCAGGCGGCG  >  NC_002947/5608153‑5608448

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: