Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 237,457 T→C 6.2% L113S (TTG→TCG)  algR → alginate biosynthesis regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947237,4570TC 6.2% 188.3 / 2.1 64L113S (TTG→TCG) algRalginate biosynthesis regulatory protein
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base C (2/2);  total (30/34)
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.49e-01

CAGGTTGCCGCTCGCCTGTGCGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAGGCATCGAACTGATCCCGTTACCCCAGGTGATCTATTTC  >  NC_002947/237309‑237599
                                                                                                                                                      |                                                                                                                                              
cAGGTTGCCGCTCGCCTGTGCGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCtt                                                                                                                                                >  6:254514/1‑149 (MQ=255)
  ggTTGCCGCTCGCCTGTGCGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTgc                                                                                                                                              <  8:51344/149‑1 (MQ=255)
   gTTGCCGCTCGCCTGTGCGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTgcg                                                                                                                                             <  2:34495/149‑1 (MQ=255)
   gTTGCCGCTCGCCTGTGCGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTgcg                                                                                                                                             >  6:106027/1‑149 (MQ=255)
        cGCTCGCCTGTGCGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGTCTTGCGCAAgg                                                                                                                                        >  4:209001/1‑149 (MQ=255)
        cGCTCGCCTGTGCGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGTCTTGCGCAAgg                                                                                                                                        >  6:187237/1‑149 (MQ=255)
               cTGTGCGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGaaa                                                                                                                                  >  3:286847/1‑148 (MQ=255)
               cTGTGCGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGaaa                                                                                                                                  >  2:24217/1‑148 (MQ=255)
                    cGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAAccc                                                                                                                             >  7:281686/1‑148 (MQ=255)
                    cGAACGCGAGGCGCCGCCCGCGGTGGTGGTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAAccc                                                                                                                             >  1:189242/1‑148 (MQ=255)
                     gAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCaa                                                                                                                           <  4:6431/149‑1 (MQ=255)
                      caccgcgggcggcGCCTCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCGAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCaa                                                                                                                           >  3:187222/8‑149 (MQ=255)
                      caccgcgggcggcGCCTCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCaa                                                                                                                           >  3:91002/8‑149 (MQ=255)
                      aaCGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAAc                                                                                                                          <  1:217787/149‑1 (MQ=255)
                       aCGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAacccaccc                                                                                                                         >  8:263598/1‑149 (MQ=255)
                       aCGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAAcc                                                                                                                         >  3:55595/1‑149 (MQ=255)
                         gcgAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACcgc                                                                                                                       >  2:324757/1‑149 (MQ=255)
                           gAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCgcgc                                                                                                                     >  2:269689/1‑149 (MQ=255)
                             ggCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCGTTGCGCGATGCCTTGCGCTCGGCCGAAAAACCCAACAGCGcc                                                                                                                    >  4:240778/1‑148 (MQ=255)
                                   gcccgcGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTgg                                                                                                             >  7:19301/1‑149 (MQ=255)
                                    cccgcGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCACCCAATTggc                                                                                                            <  5:254514/149‑1 (MQ=255)
                                     ccgcGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTggcg                                                                                                           <  3:174381/149‑1 (MQ=255)
                                          ggtggtGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGgcgc                                                                                                       <  2:251216/149‑1 (MQ=255)
                                            tggtgTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCt                                                                                                      >  3:198269/1‑148 (MQ=255)
                                              gtgtTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTa                                                                                                    <  5:171791/148‑1 (MQ=255)
                                               tgtTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCATGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTAc                                                                                                   <  2:196594/148‑1 (MQ=255)
                                                 ttttGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTAccc                                                                                                 <  2:313542/148‑1 (MQ=255)
                                                 ttttGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACccg                                                                                                >  5:200175/1‑149 (MQ=255)
                                                        ccGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCgcgc                                                                                                                     <  5:5552/121‑1 (MQ=255)
                                                        ccGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCgcgc                                                                                                                     >  6:5552/1‑121 (MQ=255)
                                                         cGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCcggccggc                                                                                          <  3:240778/147‑1 (MQ=255)
                                                         cGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCggccggcc                                                                                         <  5:253727/148‑1 (MQ=255)
                                                          ggggATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCaa                                                                                       <  6:200175/149‑1 (MQ=255)
                                                           gggATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAAt                                                                                      <  1:324757/149‑1 (MQ=255)
                                                             gatgatGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGa                                                                                    <  8:57470/149‑1 (MQ=255)
                                                             gatgatGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCAAGGCATCGCGCAAGGcgtgg                                                                                                                                   >  4:235884/1‑98 (MQ=38)
                                                             gatgatGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCAAGGCATCGCGCAAGGcgtgg                                                                                                                                   <  3:235884/102‑5 (MQ=38)
                                                              atgatgAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGaa                                                                                   <  7:61635/149‑1 (MQ=255)
                                                                gatgaATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCAAGGCATCGCGCaa                                                                                                                                          >  2:37053/1‑92 (MQ=255)
                                                                gatgaATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCAAGGCATCGCGCaa                                                                                                                                          <  1:37053/92‑1 (MQ=255)
                                                                       aCGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAggcggcggc                                                                          >  3:79506/1‑149 (MQ=255)
                                                                                agGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGc                                                                 >  3:41207/1‑149 (MQ=255)
                                                                                agGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGc                                                                 >  3:45453/1‑149 (MQ=255)
                                                                                      tCAAGGACAGTACCCTCAGCCATGTGGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCatat                                                            >  1:169575/1‑149 (MQ=255)
                                                                                        aaGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCa                                                          <  3:229660/148‑1 (MQ=255)
                                                                                         aGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAgc                                                        <  4:79506/149‑1 (MQ=255)
                                                                                         aGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAgc                                                        <  5:61983/149‑1 (MQ=255)
                                                                                          ggACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAgcg                                                       >  1:177335/1‑149 (MQ=255)
                                                                                              aGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCg                                                   <  4:55595/149‑1 (MQ=255)
                                                                                                   ccTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCAcgc                                              >  3:86473/1‑149 (MQ=255)
                                                                                                   ccTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCAcgc                                              >  2:234359/1‑149 (MQ=255)
                                                                                                        gCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCaaa                                         <  1:24217/149‑1 (MQ=255)
                                                                                                         ccATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAg                                        <  4:198269/149‑1 (MQ=255)
                                                                                                         ccATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAg                                        <  2:189242/149‑1 (MQ=255)
                                                                                                           aTGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCatat                                                            >  6:235464/1‑127 (MQ=255)
                                                                                                           aTGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCatat                                                            <  5:235464/127‑1 (MQ=255)
                                                                                                           aTGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAgg                                       >  6:324177/1‑148 (MQ=255)
                                                                                                           aTGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAGGc                                      <  2:169575/149‑1 (MQ=255)
                                                                                                           aTGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGAGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAGGc                                      <  3:9530/149‑1 (MQ=255)
                                                                                                              tgACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAGGCAt                                    <  3:228126/148‑1 (MQ=255)
                                                                                                              tgACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAGGCAt                                    <  8:281686/148‑1 (MQ=255)
                                                                                                                   cAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAGGCATCGAAc                               >  4:109141/1‑149 (MQ=255)
                                                                                                                              aGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCACGCAGCCATATCAGCGCCCGCACGCGCAAAGGCATCGAACTGATCCCGTTAc                   <  4:187222/149‑1 (MQ=255)
                                                                                                                                             gcgATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAGGCATCGAACTGATCCCGTTACCCCAGGTGATCTAt     <  8:19301/149‑1 (MQ=255)
                                                                                                                                                aTGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCATGCGCAAAGGCATCGAACTGATCCCGTTACCCCAGGTGATCTATTTc  <  1:234359/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
CAGGTTGCCGCTCGCCTGTGCGAACGCGAGGCGCCGCCCGCGGTGGTGTTTTGCTCCGGGGATGATGAATACGGTGCAGAGGCATTCAAGGACAGTACCCTCAGCCATGTGACCAAGCCCTTCCAGGCCCACGCCTTGCGCGATGCCTTGCGCAAGGCCGAAAAACCCAACCGCGCCCAATTGGCGGCGCTTACCCGGCCGGCCAATGAAGGCGGCGGCCCACGCAGCCATATCAGCGCCCGCACGCGCAAAGGCATCGAACTGATCCCGTTACCCCAGGTGATCTATTTC  >  NC_002947/237309‑237599

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: