Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,470,045 C→G 5.2% C33S (TGT→TCT)  algK ← alginate biosynthesis protein AlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,470,0450CG 5.2% 248.0 / 3.0 76C33S (TGT→TCT) algKalginate biosynthesis protein AlgK
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base G (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTGTCGATCAGCCAGTTCAGCGATGGGCCGCGGTC  >  NC_002947/1469902‑1470188
                                                                                                                                               |                                                                                                                                               
gcgTGGCCACCTTGATGTCCGCAAAGCCGAATTGCGCTTCGCTGTAGCACACGACGCCCAGCGACTTGTAGTTGCGCTCGGCCAATGACGTGTCACCGCGCTTCAGCGCCACGATCGCCACGGGCTAGTAGGGCAGGCCGGCACAGCCg                                                                                                                                            >  2:338110/3‑149 (MQ=255)
      ccACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAgg                                                                                                                                      <  3:359845/149‑1 (MQ=255)
       cACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGTGCAGGCCGGCACAGCCGGCCAggg                                                                                                                                     >  8:311972/1‑149 (MQ=255)
       cACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCTGTGCCGGcct                                                                                                                                        <  5:371556/146‑2 (MQ=255)
       cACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCTGTGCCGGcct                                                                                                                                        >  6:371556/1‑145 (MQ=255)
       cACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAggg                                                                                                                                     <  8:98631/149‑1 (MQ=255)
        aCCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGt                                                                                                                                    <  6:84818/149‑1 (MQ=255)
         ccTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGt                                                                                                                                    >  3:355534/1‑148 (MQ=255)
          cTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGa                                                                                                                                  >  2:158223/1‑149 (MQ=255)
          cTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGa                                                                                                                                  >  2:124583/1‑149 (MQ=255)
              aTGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATggc                                                                                                                              >  7:60109/1‑149 (MQ=255)
              aTGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGGTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATggc                                                                                                                              <  4:101818/149‑1 (MQ=255)
                   cGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATggcggcgg                                                                                                                         >  8:32456/1‑149 (MQ=255)
                    gCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATggcggcggc                                                                                                                        >  5:153528/1‑149 (MQ=255)
                      cAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATggcggcggc                                                                                                                        >  2:342249/1‑147 (MQ=255)
                        ggCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGt                                                                                                                    >  7:259867/1‑149 (MQ=255)
                          ccGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGc                                                                                                                  <  7:38702/149‑1 (MQ=255)
                             aCTTGCGCTTCGCTGTAGCCCAGGTCGGCCTGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGcaca                                                                                                               <  1:229724/149‑1 (MQ=255)
                             aCTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGcaca                                                                                                               >  7:251885/1‑149 (MQ=255)
                                  cgcTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAAc                                                                                                           <  3:180487/148‑1 (MQ=255)
                                   gcTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAgg                                                                                                                                      >  1:29626/1‑120 (MQ=255)
                                   gcTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAgg                                                                                                                                      <  2:29626/120‑1 (MQ=255)
                                       cGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCa                                                                                                       <  1:298066/147‑1 (MQ=255)
                                        gCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAgg                                                                                                     <  1:46271/148‑1 (MQ=255)
                                         cTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCACTCGACCTATACCGTGTCACCGCGCTTCAGAGCCAAGTTGGCCAGGCGCTGGTAGGGCAGGCCGGCACAGCCGGACAGGGTGAAGGCGCCGGCCAGTGCAAACACCCCCAGGta                                                                                                   >  8:233015/1‑148 (MQ=255)
                                           gTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCa                                                                                                 >  6:249919/1‑149 (MQ=255)
                                              gCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCAcgc                                                                                              <  7:32456/149‑1 (MQ=255)
                                               cccAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCAcgcg                                                                                             <  1:342249/149‑1 (MQ=255)
                                                 cAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGa                                                                                           <  1:124583/149‑1 (MQ=255)
                                                        gCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTg                                                                                      >  8:69174/1‑147 (MQ=255)
                                                        gCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTa                                                                                    <  8:251885/149‑1 (MQ=255)
                                                         cgaGCGGCTTGTTGTTGCGGTCGGGCATTGCCGTGTCTCCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTTGTCGGGCAGGCCGGCACAGCCGGCCAGGTTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAg                                                                                   <  7:311972/147‑1 (MQ=255)
                                                               ccTTGTTGTTGCGCTGGGCCTATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCa                                                                             <  7:233015/149‑1 (MQ=255)
                                                               ccTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTAt                                                                               >  7:345264/1‑147 (MQ=255)
                                                                cTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTAt                                                                               >  3:78038/1‑146 (MQ=255)
                                                                         gctcTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCAtt                                                                    <  8:80220/145‑1 (MQ=255)
                                                                         gcgcTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCAtt                                                                    <  5:249919/148‑1 (MQ=255)
                                                                         gcgcTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCAtt                                                                    <  8:60109/148‑1 (MQ=255)
                                                                         gcgcTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTg                                                                   >  1:241071/1‑149 (MQ=255)
                                                                          cgcTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTgt                                                                  <  8:345264/149‑1 (MQ=255)
                                                                           gcTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAAcc                                                                                                          <  4:249449/108‑1 (MQ=255)
                                                                           gcTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAAcc                                                                                                          >  3:249449/1‑108 (MQ=255)
                                                                           gcTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTAGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTgtc                                                                 >  4:347986/1‑149 (MQ=255)
                                                                               ggCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGa                                                                                           <  2:242648/119‑1 (MQ=255)
                                                                               ggCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGa                                                                                           >  1:242648/1‑119 (MQ=255)
                                                                               ggCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGAcc                                                                                         >  7:189043/1‑121 (MQ=255)
                                                                               ggCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGAcc                                                                                         <  8:189043/121‑1 (MQ=255)
                                                                               ggCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTgtccgtc                                                             <  4:78038/149‑1 (MQ=255)
                                                                                      gCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCAcgcg                                                                                             <  4:251516/110‑1 (MQ=255)
                                                                                      gCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCAcgcg                                                                                             >  3:251516/1‑110 (MQ=255)
                                                                                      gCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTca                                                        >  3:309413/1‑147 (MQ=255)
                                                                                       ccGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTcaca                                                      <  7:69174/148‑1 (MQ=255)
                                                                                       ccGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAg                                                     <  2:241071/149‑1 (MQ=255)
                                                                                       ccGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGACCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAg                                                     <  4:338159/149‑1 (MQ=255)
                                                                                              aCCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAggccacgg                                              >  1:265444/1‑149 (MQ=255)
                                                                                              aCCGCGCTTCAGGGCATCGTTGGCCAGGCTCTGTTGGGGCAGGCCTGCACAGCCGGCCAGTGTGATGGCGGCTGCCAGTGCACCCAACCGCAGTTTCACGCGGACATTGTATGTATCAGAGATCATTGTCCGTCCCTCACAggcaacgg                                              <  1:338110/149‑1 (MQ=255)
                                                                                               ccGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAggccacggg                                             >  2:161402/1‑149 (MQ=255)
                                                                                                    cTTCAGGGCCTCGTAGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAggccacgggccacg                                        <  3:347986/149‑1 (MQ=255)
                                                                                                     ttCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGgccacgggccacgg                                       >  8:169676/1‑149 (MQ=255)
                                                                                                       cAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACCCACCCCCAGGTGCACGCGGCCGTGGTAGGTAGCAGAGAGCATGTGCCGGCCGGCACAGGCCCCGGGCCCCGGcg                                     >  8:114682/1‑148 (MQ=255)
                                                                                                              tCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTGTCGa                              >  7:330105/1‑149 (MQ=255)
                                                                                                                 ttGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTGTCGATCa                           <  1:158223/149‑1 (MQ=255)
                                                                                                                 ttGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTGTCGATCa                           <  6:153528/149‑1 (MQ=255)
                                                                                                                    gCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCCCGGGCCACGGCTGTGGCGAtgcgc                         >  8:343864/1‑144 (MQ=255)
                                                                                                                    gCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTGTCGATCAGc                         <  2:265444/148‑1 (MQ=255)
                                                                                                                     ccAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTGTCGATCAGc                         >  6:27617/1‑147 (MQ=255)
                                                                                                                       aGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAAcccc                                                                                                        >  7:33454/1‑66 (MQ=255)
                                                                                                                       aGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAAcccc                                                                                                        <  8:33454/66‑1 (MQ=255)
                                                                                                                        ggCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTGTCGATCAGCCAGtt                    <  1:161402/149‑1 (MQ=255)
                                                                                                                         gcgcTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTGACAGGCCACGGGCCACGGCTTTGTCGATCAGCCAGtt                    >  7:91640/1‑148 (MQ=255)
                                                                                                                                    gCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTcac                                                       <  2:37237/102‑1 (MQ=255)
                                                                                                                                    gCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTcac                                                       >  1:37237/1‑102 (MQ=255)
                                                                                                                                      tggCCGGCTGTGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTg                                  >  4:4/2‑121 (MQ=255)
                                                                                                                                      tggCCGGCTGTGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTg                                  <  3:4/120‑1 (MQ=255)
                                                                                                                                        gCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTGTCGATCAGCCAGTTCAGCGATGGGCCGCgg    >  8:269972/1‑149 (MQ=255)
                                                                                                                                          cGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTGTCGATCAGCCAGTTCAGCGATGGGCCGCGGTc  >  2:251522/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGCCAGCGCCTTGTAGTTGCGCTCGGCCAATGCCGTGTCACCGCGCTTCAGGGCCTCGTTGGCCAGGCGCTGGTCGGGCAGGCCGGCACAGCCGGCCAGGGTGATGGCGGCGGCCAGTGCACACAACCCCAGGTTCACGCGGACCTTGTAGGTATCAGAGATCATTGTCCGTCCCTCACAGGCCACGGGCCACGGCTTTGTCGATCAGCCAGTTCAGCGATGGGCCGCGGTC  >  NC_002947/1469902‑1470188

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: