Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 138,949 G→A 6.9% 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,9490GA 6.9% 251.3 / 5.6 88R519W (CGG→TGG) kinBalginate biosynthesis sensor protein KinB
Reads supporting (aligned to +/- strand):  ref base G (41/39);  new base A (3/3);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-02

TTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCG  >  NC_002947/138805‑139088
                                                                                                                                                |                                                                                                                                            
tttCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTg                                                                                                                                          >  8:58069/1‑149 (MQ=255)
   gcTTGCACAGCGCGAGCCCGAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGTAAGGCTCTAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGc                                                                                                                                        <  5:148258/147‑1 (MQ=255)
   ccTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCa                                                                                                                                       >  8:145124/1‑149 (MQ=255)
       gCACAGCGCGAGCCCCAGCCCAGCGCCGACTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCACTCTCCGTTGTCCTCGACGCTGCTGATAACCCGTTCTCCATGCCGTcagccaag                                                                                                                                     >  3:226391/1‑140 (MQ=255)
        cacaGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCAGACGGCATGCa                                                                                                                                   <  2:299818/148‑1 (MQ=255)
         acaGCGCGAGCCCCAGCCCAGTGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCCTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATCCAgg                                                                                                                                 >  4:281790/1‑149 (MQ=255)
         acaGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAgg                                                                                                                                 >  2:245538/1‑149 (MQ=255)
          caGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGc                                                                                                                                >  5:261177/1‑149 (MQ=255)
            gcgcGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCg                                                                                                                               >  4:151991/1‑148 (MQ=255)
                        gcccagcGCCGCCTTTCTTGCGACGTACCTCCACCAAAGGCTCGAACATGCGCGCTTGTTGCCCATAGGAAATCCCCCCTCCGTTGTCCTCGACGCTCATGATCACCGGTTCTCCATGCCGTCTGGCATGCAGGCGCATCcgcctgcat                                                                                                                  <  2:327525/149‑3 (MQ=255)
                          ccagcGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGtt                                                                                                                >  6:208370/1‑149 (MQ=255)
                           cagcGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTg                                                                                                               <  1:245538/149‑1 (MQ=255)
                               gccgccTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCgg                                                                                                          >  8:98210/1‑149 (MQ=255)
                                   ccTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATg                                                                                                      >  7:213099/1‑149 (MQ=255)
                                    cTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATgg                                                                                                     >  3:186471/1‑149 (MQ=255)
                                     tttCTTGCGACCTAGCTGCAGGAAAGGCTCGAAGATGCGCTCTTGTTGCCCATTGGAAATGCCCTGTCCGTTGTCCTCGACGCTGATGATCACCGGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGc                                                                                                    <  4:174176/149‑1 (MQ=255)
                                         ttGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGc                                                                                                    >  1:19906/1‑145 (MQ=255)
                                         ttGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGc                                                                                                    <  2:19906/145‑1 (MQ=255)
                                              acctacctGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGt                                                                                           >  6:84213/1‑149 (MQ=255)
                                              acctacctGCACGAAAGGCTCGAAGATGCCCGCTTGTAGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGt                                                                                           >  6:94303/1‑149 (MQ=255)
                                               cctacctGCACGTAAGGCTCGAAGCTGCGCGCTTGTTGCCCATAGGAAATGCCCTCTGCGTTGTCCTGGACGCTGATGATCACCCGTTGTGCATGCCGTCTGTCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGtt                                                                                          <  3:281790/149‑1 (MQ=255)
                                               cctacctGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGtt                                                                                          <  2:214954/149‑1 (MQ=255)
                                               cctacctGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGtt                                                                                          <  6:261177/149‑1 (MQ=255)
                                                   cctGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTgc                                                                                      <  3:151991/149‑1 (MQ=255)
                                                   cctGCACGAAAGGCTCGAAGATGCGCGCTTGGTGCCCATAGGAAATGCCCTCTCCGTTCTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTtgt                                                                                        <  5:216945/147‑1 (MQ=255)
                                                    ctGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTtgtg                                                                                       >  6:109669/1‑147 (MQ=255)
                                                       cACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTgcagca                                                                                  <  7:58069/149‑1 (MQ=255)
                                                       cACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTgcagca                                                                                  <  5:208370/149‑1 (MQ=255)
                                                       cACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTgcag                                                                                    >  8:172782/1‑147 (MQ=255)
                                                          gAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTgcagca                                                                                  >  8:181243/1‑146 (MQ=255)
                                                          gAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTgcagca                                                                                  <  7:181243/146‑1 (MQ=255)
                                                          gAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGt                                                                               >  2:49851/1‑149 (MQ=255)
                                                          gAAAGGCTCGAAGATGCGCGATTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCAGTTCTGCATGCCGTCTGGCATGCAGGCGGATGAGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGt                                                                               >  8:238203/1‑149 (MQ=255)
                                                            aaGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTg                                                                             <  8:213099/149‑1 (MQ=255)
                                                               gCTCGAAGATGCGCGCTTGTTGCCCATAGTAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTc                                                                           <  4:253522/148‑1 (MQ=255)
                                                               gCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTcc                                                                          <  1:49851/149‑1 (MQ=255)
                                                               gCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTcc                                                                          <  4:186471/149‑1 (MQ=255)
                                                               gCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTc                                                                           >  5:278565/1‑148 (MQ=255)
                                                                cTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTc                                                                           >  3:270232/1‑147 (MQ=255)
                                                                cTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCa                                                                         <  6:13743/149‑1 (MQ=255)
                                                                 tCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTg                                                                             >  4:224851/1‑144 (MQ=255)
                                                                 tCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTg                                                                             <  3:224851/144‑1 (MQ=255)
                                                                 tCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAg                                                                        >  4:221931/1‑149 (MQ=255)
                                                                      gATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTc                                                                   >  8:85963/1‑149 (MQ=255)
                                                                           gcgcTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTc                                                               >  7:16429/1‑148 (MQ=255)
                                                                               ttgttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACt                                                          >  5:154980/1‑149 (MQ=255)
                                                                               ttgttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACt                                                          <  7:145124/149‑1 (MQ=255)
                                                                                 gttgCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTgg                                                        >  2:215622/1‑149 (MQ=255)
                                                                                    gCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGc                                                      <  5:109669/148‑1 (MQ=255)
                                                                                      ccATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGAc                                                   >  4:51247/1‑149 (MQ=255)
                                                                                       cATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACg                                                  >  5:13326/1‑149 (MQ=255)
                                                                                        aTAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGt                                                 >  5:193464/1‑149 (MQ=255)
                                                                                                   ccTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGc                                             >  8:270316/1‑142 (MQ=255)
                                                                                                   ccTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGc                                             <  7:270316/142‑1 (MQ=255)
                                                                                                     tctcCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGc                                             <  2:228273/140‑1 (MQ=255)
                                                                                                     tctcCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGc                                             >  1:228273/1‑140 (MQ=255)
                                                                                                         cGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATAcgcg                                  <  6:278565/147‑1 (MQ=255)
                                                                                                         cGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATAcgc                                   <  6:154980/146‑1 (MQ=255)
                                                                                                             gTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAgc                             <  8:16429/148‑1 (MQ=255)
                                                                                                              tCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAgcg                            <  8:119502/148‑1 (MQ=255)
                                                                                                                cTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGTTCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAgcggc                          >  1:260569/1‑149 (MQ=255)
                                                                                                                cTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAgcggcg                         >  5:195904/1‑149 (MQ=255)
                                                                                                                 tCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGAGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGcggcgg                        >  5:230546/1‑149 (MQ=255)
                                                                                                                       cTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAg                  <  7:98210/149‑1 (MQ=255)
                                                                                                                          atgatCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTc               <  3:261409/149‑1 (MQ=255)
                                                                                                                           tgatCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTcc              <  6:195904/149‑1 (MQ=255)
                                                                                                                              tCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCtt            <  3:51247/148‑1 (MQ=255)
                                                                                                                               cACCCGTTCTGCTTGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTGGTGCAGCAGGTTGTCCTGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTg           <  5:94303/148‑1 (MQ=255)
                                                                                                                               cACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCg                                                     >  8:292793/1‑106 (MQ=255)
                                                                                                                               cACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCg                                                     <  7:292793/106‑1 (MQ=255)
                                                                                                                                aCCCGTTCTGCTTGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGAt         <  2:201203/149‑1 (MQ=255)
                                                                                                                                aCCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGa          <  3:276677/148‑1 (MQ=255)
                                                                                                                                 cccGTTGTGCATGCCTGCGGGCATGCAGGCGGATGCTGGGGCCGTTGGCGGTATGGCGGTTGGCGTTGTGCAGCAGGTTTTACAGCACTCTGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGAta        <  4:226391/149‑2 (MQ=255)
                                                                                                                                 cccGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCg                                               >  2:250868/1‑110 (MQ=255)
                                                                                                                                 cccGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCg                                               <  1:250868/110‑1 (MQ=255)
                                                                                                                                 cccGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTgct         >  3:117526/1‑146 (MQ=255)
                                                                                                                                    gTTCTGCATGCCTTCTGGCATGCAGGCGTATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAAc     <  7:172782/149‑1 (MQ=255)
                                                                                                                                     gccTGCATGCCAGACGGCATGCAGACGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGa                                        >  4:239397/3‑113 (MQ=255)
                                                                                                                                     gccTGCATGCCAGACGGCATGCAGACGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGa                                        <  3:239397/111‑1 (MQ=255)
                                                                                                                                     gccTGCATGCCAGACGGCATGCAGACGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACt    >  7:97891/3‑149 (MQ=255)
                                                                                                                                       cTGCATGCCAGACGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTg           <  4:108916/140‑1 (MQ=255)
                                                                                                                                       cTGCATGCCAGACGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTg           >  3:108916/1‑140 (MQ=255)
                                                                                                                                        tGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCt                                           >  6:83208/1‑107 (MQ=255)
                                                                                                                                        tGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCt                                           <  5:83208/107‑1 (MQ=255)
                                                                                                                                         gCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTc                     >  1:128267/1‑128 (MQ=255)
                                                                                                                                          cATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCg  <  7:85963/146‑1 (MQ=255)
                                                                                                                                             gCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAgcggc                          >  6:96104/1‑119 (MQ=255)
                                                                                                                                             gCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAgcggc                          <  5:96104/119‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
TTTCCTTGCACAGCGCGAGCCCCAGCCCAGCGCCGCCTTTCTTGCGACCTACCTGCACGAAAGGCTCGAAGATGCGCGCTTGTTGCCCATAGGAAATGCCCTCTCCGTTGTCCTCGACGCTGATGATCACCCGTTCTGCATGCCGTCTGGCATGCAGGCGGATGCGCCCGCCGTTGGCGGTATGGCGGATGGCGTTGTGCAGCAGGTTGTCCAGCACTCGGTCCAACTGGGCGACGTCGGCCTGGATACGCGGTAGCGGCGGGTCCAGTTCCTTGATCAACTCG  >  NC_002947/138805‑139088

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