Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 138,949 G→A 11.2% R519W (CGG→TGG)  kinB ← alginate biosynthesis sensor protein KinB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947138,9490GA11.2% 260.9 / 18.2 98R519W (CGG→TGG) kinBalginate biosynthesis sensor protein KinB
Reads supporting (aligned to +/- strand):  ref base G (46/41);  new base A (6/5);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

ACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCGATCTGC  >  NC_002947/138801‑139094
                                                                                                                                                    |                                                                                                                                                  
acaATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGCTGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGCCGCTGATGATCACCCGTTCTGCATGCCg                                                                                                                                                    >  7:339437/1‑149 (MQ=255)
 caATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGt                                                                                                                                                   >  4:105093/1‑149 (MQ=255)
        cTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCAt                                                                                                                                           <  8:235172/149‑1 (MQ=255)
         ttGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATg                                                                                                                                          >  4:206048/1‑149 (MQ=255)
            cacaGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCAGACGGCATGCa                                                                                                                                        <  5:245984/148‑1 (MQ=255)
             acaGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCATCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAgg                                                                                                                                      >  6:361516/1‑149 (MQ=255)
              caGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGc                                                                                                                                     >  1:246826/1‑149 (MQ=255)
               aGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCg                                                                                                                                    >  4:126722/1‑149 (MQ=255)
               aGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCg                                                                                                                                    >  1:77964/1‑149 (MQ=255)
               aGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTCCTGCATGCCAGACGGCATGCAgaac                                                                                                                                    <  6:375495/149‑4 (MQ=255)
                gcgcGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTCCTGCATGCCAGACGGCATGCAgaacg                                                                                                                                   <  5:189731/149‑5 (MQ=255)
                gcgcGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTCCTGCATGCCAGACGGCATGCAgaacg                                                                                                                                   >  6:189731/1‑145 (MQ=255)
                  gcgAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGAt                                                                                                                                 <  3:134173/149‑1 (MQ=255)
                  gcgAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGAt                                                                                                                                 <  5:289228/149‑1 (MQ=255)
                    gAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATgc                                                                                                                               <  3:105093/149‑1 (MQ=255)
                        cccagcccagCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCAGACGGCATGCa                                                                                                                                        <  6:412453/136‑1 (MQ=255)
                        cccagcccagCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCAGACGGCATGCa                                                                                                                                        >  5:412453/1‑136 (MQ=255)
                          cagcccagcGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCAGACGGCATGCa                                                                                                                                        >  8:374782/1‑134 (MQ=255)
                          cagcccagcGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCAGACGGCATGCa                                                                                                                                        <  7:374782/134‑1 (MQ=255)
                              ccagcGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGAATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGtt                                                                                                                     >  8:369639/1‑149 (MQ=255)
                               cagcGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATACAGGCGGATGCGCCCGCCGTTg                                                                                                                    <  1:187940/149‑1 (MQ=255)
                                agcGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGGTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTgg                                                                                                                   <  6:38658/149‑1 (MQ=255)
                                  cgccgcCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCg                                                                                                                 <  8:243426/149‑1 (MQ=255)
                                  cgccgcCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCg                                                                                                                 <  3:206048/149‑1 (MQ=255)
                                       ccTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCAAGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATg                                                                                                            <  6:225153/149‑1 (MQ=255)
                                             ttGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTTGAGATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGa                                                                                                       >  1:135172/1‑149 (MQ=255)
                                             ttGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTATCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGAt                                                                                                      >  7:323878/1‑149 (MQ=255)
                                               gCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGa                                                                                                       <  1:400488/146‑1 (MQ=255)
                                               gCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATgg                                                                                                    >  4:371622/1‑149 (MQ=255)
                                                gactctacctGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGc                                                                                                   <  1:144851/145‑1 (MQ=255)
                                                     tacctGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTtgt                                                                                              >  6:180392/1‑149 (MQ=255)
                                                       cctGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTtgtg                                                                                             <  1:172750/148‑1 (MQ=255)
                                                       cctGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTgc                                                                                            <  3:126722/149‑1 (MQ=255)
                                                        ctGCACTAAAGGCTCGAAGATGAGCGCTTGTTGCCCATAGGAAATGCCCTCTCAGTTGTCCTCGACGATGATGATCACCCGTTCTGCATCACGTCTGGAATGCAGCAGGATGCGCCCGCCCATGGCGGTATGGCGGATGGCGTtgtg                                                                                             >  3:289870/1‑147 (MQ=255)
                                                         tGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTg                                                                                                                                               <  2:86962/96‑1 (MQ=255)
                                                         tGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTg                                                                                                                                               >  1:86962/1‑96 (MQ=255)
                                                             cGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAgg                                                                                      <  4:394633/149‑1 (MQ=255)
                                                              gAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGt                                                                                     >  6:87190/1‑149 (MQ=255)
                                                                aaGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTtgtg                                                                                             >  1:202054/1‑139 (MQ=255)
                                                                aaGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGAATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCCGGTTg                                                                                   <  7:369639/149‑1 (MQ=255)
                                                                   gCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGTCGTTGTGCAGCAGGTTGTc                                                                                 >  8:385498/1‑148 (MQ=255)
                                                                   gCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTcc                                                                                <  1:261774/149‑1 (MQ=255)
                                                                    cTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCg                                                                                                  <  8:292130/130‑1 (MQ=255)
                                                                    cTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCg                                                                                                  >  7:292130/1‑130 (MQ=255)
                                                                    cTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCa                                                                               <  2:254016/149‑1 (MQ=255)
                                                                     tCGAAGATGCGCGCTTGTTGCCCGTAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAg                                                                              >  1:234790/1‑149 (MQ=255)
                                                                     tCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAg                                                                              >  3:18217/1‑149 (MQ=255)
                                                                     tCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAg                                                                              >  4:199712/1‑149 (MQ=255)
                                                                               gcgcTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTc                                                                     >  3:179063/1‑148 (MQ=255)
                                                                                   ttgttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACt                                                                >  2:119235/1‑149 (MQ=255)
                                                                                      ttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTggg                                                             >  1:257428/1‑149 (MQ=255)
                                                                                      ttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTgg                                                              <  4:179063/148‑1 (MQ=255)
                                                                                      ttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTgg                                                              <  5:87190/148‑1 (MQ=255)
                                                                                          ccATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGAc                                                         >  3:350985/1‑149 (MQ=255)
                                                                                            aTAGGAAATGCCCTGTCCGTGGGCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTTTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGt                                                       <  4:289870/149‑1 (MQ=255)
                                                                                               ggAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCg                                                     <  4:18217/148‑1 (MQ=255)
                                                                                                gAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGc                                                   <  5:180392/149‑1 (MQ=255)
                                                                                                gAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGc                                                   <  2:246826/149‑1 (MQ=255)
                                                                                                 aaaTGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGcc                                                  >  7:352669/1‑149 (MQ=255)
                                                                                                  aaTGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCt                                                 <  2:234790/149‑1 (MQ=255)
                                                                                                      cccTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGa                                              >  8:133574/1‑148 (MQ=255)
                                                                                                        ctctCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATAc                                           >  6:265294/1‑149 (MQ=255)
                                                                                                        ctctCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGAACTCGGTCCCACTGGGCGCCGTCGGCCTGGATAc                                           >  8:119561/1‑149 (MQ=255)
                                                                                                         tctcCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATAcg                                          >  8:2533/1‑149 (MQ=255)
                                                                                                          ctcCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATAcgc                                         >  5:171222/1‑149 (MQ=255)
                                                                                                           tcCGTTGTCCTCGACGCTGATGATCACCCATTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATAcgcg                                        >  7:346402/1‑149 (MQ=255)
                                                                                                             cGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCg                                                                                                                 <  2:303024/74‑1 (MQ=255)
                                                                                                             cGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCg                                                                                                                 >  1:303024/1‑74 (MQ=255)
                                                                                                             cGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAAc                                                                 >  1:186286/1‑122 (MQ=255)
                                                                                                             cGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAAc                                                                 <  2:186286/122‑1 (MQ=255)
                                                                                                             cGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATAcgcg                                        <  2:135172/147‑1 (MQ=255)
                                                                                                                    cTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAgcggcg                               >  1:57931/1‑149 (MQ=255)
                                                                                                                    cTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAgcggcg                               >  4:8775/1‑149 (MQ=255)
                                                                                                                    cTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAgcggcg                               >  3:395357/1‑149 (MQ=255)
                                                                                                                       gACGCTGATGATCCCCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCATGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGt                            <  8:323878/149‑1 (MQ=255)
                                                                                                                       gACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCggg                             <  2:25192/148‑1 (MQ=255)
                                                                                                                            tgatgaTCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGt                       >  2:172955/1‑149 (MQ=255)
                                                                                                                            tgatgaTCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGt                       >  1:141050/1‑149 (MQ=255)
                                                                                                                                gatCACCCTTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTcc                    <  8:346402/148‑1 (MQ=255)
                                                                                                                                gatCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCCGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGAGGCGGGTCCACTTACt                   >  8:179343/1‑149 (MQ=255)
                                                                                                                                 atCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCtt                  <  1:172955/149‑1 (MQ=255)
                                                                                                                                 atCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCt                   <  5:361516/148‑1 (MQ=255)
                                                                                                                                  tCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCtt                  >  6:28967/1‑148 (MQ=255)
                                                                                                                                  tCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCtt                  <  5:28967/148‑1 (MQ=255)
                                                                                                                                  tCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTg                 >  1:215496/1‑149 (MQ=255)
                                                                                                                                  tCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTg                 <  5:265294/149‑1 (MQ=255)
                                                                                                                                  tCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTg                 >  1:45930/1‑149 (MQ=255)
                                                                                                                                    aCCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCaa                                                                  >  1:132694/1‑98 (MQ=255)
                                                                                                                                    aCCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCaa                                                                  <  2:132694/98‑1 (MQ=255)
                                                                                                                                        gTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTcc                    >  4:132114/1‑140 (MQ=255)
                                                                                                                                        gTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTcc                    <  3:132114/140‑1 (MQ=255)
                                                                                                                                           cTGCATGCCAGACGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTc         >  8:230358/1‑148 (MQ=255)
                                                                                                                                           cTGCATGCCAGACGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTc         >  6:307775/1‑148 (MQ=255)
                                                                                                                                           cTGCATGCCAGACGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGAAGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCg        >  5:270773/1‑149 (MQ=255)
                                                                                                                                              cATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCGAt      <  2:257428/148‑1 (MQ=255)
                                                                                                                                              cATGCCGTCTGGCATGCAGGCGGATGCGCACGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATCACCGGTAGCGGCGGGTGCACTTCCTAGATCAACTCgct      >  8:141887/1‑146 (MQ=255)
                                                                                                                                                 gCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCGATCTGc  >  6:73640/1‑149 (MQ=255)
                                                                                                                                                  ccGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCGATCTGc  <  3:8775/148‑1 (MQ=255)
                                                                                                                                                  ccGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCGATCTGc  <  4:395357/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACAATTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCGATCTGC  >  NC_002947/138801‑139094

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: