Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,469,147 T→G 7.4% Q476P (CAG→CCG)  algX ← alginate biosynthesis protein AlgX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,469,1470TG 7.4% 654.3 / 6.8 276Q476P (CAG→CCG) algXalginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (133/108);  new base G (4/25);  total (143/133)
Fisher's exact test for biased strand distribution p-value = 1.89e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

TCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGG  >  NC_002947/1468998‑1469296
                                                                                                                                                     |                                                                                                                                                     
  aCCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCctgcctgc                                                                                                                                                     >  1:141054/1‑150 (MQ=255)
  aCCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGActgcctgc                                                                                                                                                     >  1:5777933/1‑150 (MQ=255)
    cGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTgcgc                                                                                                                                                   <  1:5168870/150‑1 (MQ=255)
    cGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTgcgc                                                                                                                                                   <  2:4054796/150‑1 (MQ=255)
       tCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGt                                                                                                                                                <  2:1438969/150‑1 (MQ=255)
       tCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGt                                                                                                                                                >  1:4345714/1‑150 (MQ=255)
        cTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTc                                                                                                                                               <  2:2612121/150‑1 (MQ=255)
            aGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcg                                                                                                                                           >  2:8668383/1‑150 (MQ=255)
              cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCACACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                         <  2:7318683/150‑1 (MQ=255)
              cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTGTgcgct                                                                                                                                         >  1:7601318/1‑149 (MQ=255)
              cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcgcc                                                                                                                                         >  1:6281858/1‑149 (MQ=255)
              cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                         <  1:7368441/150‑1 (MQ=255)
               ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                        <  2:3369195/150‑1 (MQ=255)
               ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                        >  2:7449091/1‑150 (MQ=255)
               ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                        >  2:5819080/1‑150 (MQ=255)
               ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                        >  2:1032729/1‑150 (MQ=255)
                tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGGCTGCGCAgg                                                                                                                                       >  1:2442543/1‑150 (MQ=255)
                   cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCAc                                                                                                                                    >  2:3799968/1‑150 (MQ=255)
                   cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCAc                                                                                                                                    <  1:4373459/150‑1 (MQ=255)
                    gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                   >  2:466417/1‑150 (MQ=255)
                    gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                   >  2:5323723/1‑150 (MQ=255)
                     aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCCCgg                                                                                                                                  >  1:8494650/1‑150 (MQ=255)
                         ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGTAAATTCGTTTTAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                              <  2:6508832/150‑1 (MQ=255)
                         ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGac                                                                                                                              >  2:1492680/1‑149 (MQ=255)
                         ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                              >  1:8023579/1‑150 (MQ=255)
                         ccTGGTAGTAACCCGGCGGCGGTACCAGCGCGGCGGGCGCGGCACCACCAACCACGGCCAACGCGAGCCGGGGGGGGCTGTAATTTCGTTTGAAACGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGGGGTCTGCGCCGGCACGGCGaa                                                                                                                              <  2:5700865/150‑1 (MQ=255)
                            ggTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGCAgtt                                                                                                                           >  1:1124571/1‑150 (MQ=255)
                             gtagtaCCCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACACGGCCACCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAACGGGTCATGGGTCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttg                                                                                                                          <  2:7953804/150‑1 (MQ=255)
                               tgtcaCCCTGCGGCGGTACCCGCGCGGCGTGCGCGGCACCACCAAACACGGCCAAGGCGAGCCGGGGGGGGCTGGAAATTCGTTTGAAATGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAgttgtt                                                                                                                        <  2:8068660/146‑1 (MQ=255)
                                 taACCCTGGGGCGGTACCCGCGCGGCGTGCGCGGCACCACCAACCACGGCCACCGCGAGCCGGGGGGGGCTGGAAATTCGTTTGAACAGGGTCCTGGGTCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGc                                                                                                                      <  2:4191168/150‑1 (MQ=255)
                                   aCCCTGCGGCGTTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCAAACGCGAGCAGGGCGGGGCTGGAAATTCTTTTGAAAAGGGTCTTGGGTCGCCTCACTGGCCTGCCTGCGGGGTCTGCGCCGGCACGGCGAAGTTGTTGCGt                                                                                                                    <  2:3562588/150‑1 (MQ=255)
                                     ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                  >  2:2528377/1‑150 (MQ=255)
                                     ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                  >  2:1801486/1‑150 (MQ=255)
                                     ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                  >  1:5711293/1‑150 (MQ=255)
                                           gcggTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCagag                                                                                                            >  1:5561540/1‑150 (MQ=255)
                                            gggcaCCCGCGCGGCGTGCGCGCCCCCACCAAACCCGGCCAACGCGCGCCGGGGGGGGCTGCAATTTGTTTTGAAATGGTTCTTGGGGCGCCTCACTGGCGGGCCTGCGGGGTCTGCGCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGc                                                                                                           <  2:2508974/146‑1 (MQ=255)
                                            gggcaCCCGCGCGGCGTCCGCGGCCCCACCAAACACGCCCACCCGGCGCCGGGGGGGGCTGCAATTCCGTTTGAAACGGGTCTTGGGCCGCCCCACGGGCCGGCGGGGGCGGCCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGc                                                                                                           <  2:434427/146‑1 (MQ=255)
                                               tACCCGCGCGGCGTGCGCGCCACCACAAAACACGGCCACCGCGAGCCGGGCGGGGCTTTAATTTCTTTTGAAAAGGGTCTTGGGTCGCCTCACGGGCCTGCCGGCGGGGCCGGCCCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTg                                                                                                        <  2:5953008/150‑1 (MQ=255)
                                               tACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTg                                                                                                        >  1:6518019/1‑150 (MQ=255)
                                                  cAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGACATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCt                                                                                                     >  2:7840337/1‑150 (MQ=255)
                                                   tgcgcGGCGTGCGGGGCACCACCAAACACGGCCAACGCGAGCTGGGCGGGGCTGTAATTTGGTTTTAAAAGGGTCTGGGGTCGCCTCCCGGGCCGGCGTGCGCGGTCGGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCtt                                                                                                    <  2:4911184/149‑1 (MQ=255)
                                                   cgcggGGGGTGCGCGCCACCACAAAACACGGCCACCGCGAGCTGGGGGGGGCTGTAAATTGTTTTTAATTGGTTATTGGGTCGCCACACGGGCCTGCCTGGGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCtt                                                                                                    <  2:315492/149‑1 (MQ=255)
                                                    gcgcGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCAGCGCGGTCTGCGCAGGCACGGCGCAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                   >  1:2381488/1‑150 (MQ=255)
                                                            ggcgcGGAACCACCAACAAAGGCCAACGCGAGCTGGGCGGGGCTGGAAATTCGTTTGAAACGGGTCAGGGGTCGCCTCACTGGCCTGCCTGCGCGGCCTGCGAAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCt                                                                                           <  2:1993672/149‑1 (MQ=255)
                                                             gcgcgcCACCACCAAACACGCCCACCGCGAGCCGGGGGGGGCTGCAAATTCGTTTGAAATGGGTCATGGGCCGCCACACGGGCCGGCCTGCGGGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                          <  2:6223916/150‑1 (MQ=255)
                                                             gcgcGGCACCACCAAACACGCCAAACGCGAGCTGGGGGGGGCTTTAAATTCTTTTTAAATGGTTCTGGGGTCGCCTCACGGGCGTGCCTGCGCGTTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                          <  1:7288531/150‑1 (MQ=255)
                                                              cgcgCCACCACAAAACACGCCCACCGCGAGCCGGGCGGGGGTGTAAATTCTTTTGAAAAGGGTCTTGGGTCGCCTCACTGCCCGGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                         <  2:1586025/150‑1 (MQ=255)
                                                                cgccaccacCAAACACGGCCACGGCGCGCCGGGGGGGGCTTTAATTTCTTTTTAAACGGTTCTTGGGGCGCCCCACGGGCCTGCCGGCGCGGTCGGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGt                                                                                       <  1:1974970/147‑1 (MQ=255)
                                                                  gCACCACCAAACACGGCCAACGCGAGCTGGGCGGGGCTGCAATTTCTTTTGAAACGGGTCTTGGGTCGCCTCAGTGGCCGGCCGGCGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGc                                                                                     <  1:7224063/150‑1 (MQ=255)
                                                                  gCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCAGGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGACGTGGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGc                                                                                     >  2:5005294/1‑150 (MQ=255)
                                                                       accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                                >  2:237054/1‑150 (MQ=255)
                                                                       accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGAGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                                >  1:6068888/1‑150 (MQ=255)
                                                                        ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCttt                                                                               >  2:8435203/1‑150 (MQ=255)
                                                                         caAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGGCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCATGGCTTAGAACTTCTGGGTGCCGCGTTc                                                                              >  1:387488/1‑150 (MQ=255)
                                                                          aaccacggCCACCGCGAGCAGGGCGGGGCTGGAAATGCGTTTGAAACGGGTCATGGGGCGCCCCACGGGCCGGCCTGCGGGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTcc                                                                             <  2:1892769/144‑1 (MQ=255)
                                                                           aaCACGGCAACCGCGCGCTGGGGGGGGCTGTAAATTCTTTTCAACTGGGCCTGGGGCCCCCCCCGGGGCCGGCCGGGGGGTTCTGCGCGGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCg                                                                            <  1:330582/150‑1 (MQ=255)
                                                                             cacggCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCGGCCTGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGc                                                                          <  2:5096201/147‑1 (MQ=255)
                                                                              ccggCCAACGCGCGCTGGGGGGGGGTGTAATTTCTTTTGAAACGGGCCCTGGGTCGCCTCCTGGGCCGGCCGGCGGGTTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGcc                                                                         <  2:2183811/148‑1 (MQ=255)
                                                                              aaGGCAAACGCGTGCTGGGGGGGGCTTTAATTTCTTTTTAAATGGTTCTTGGTTCGCATAACTGGCCGGCGTGCGGGTTTTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGcc                                                                         <  1:2703845/150‑1 (MQ=255)
                                                                               cggCCACCGCGCGCGGGGGGGGGGTGTAATTTCTTTTGAAATGGTTCTGGGGTCGCCCCACTGGCGGGCGGGGGGGGTGGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                        <  1:3239671/149‑1 (MQ=255)
                                                                                  ccAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTg                                                                     >  1:1352694/1‑150 (MQ=255)
                                                                                   cAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGt                                                                    >  2:5249530/1‑150 (MQ=255)
                                                                                       gcgAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCt                                                                >  2:6794178/1‑150 (MQ=255)
                                                                                         gAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCg                                                              >  2:6952200/1‑150 (MQ=255)
                                                                                           gctgggCGGGGCTGGAAATTCGTTTGAAATGGGTCATGGGTCGCCTCACGGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGaa                                                            <  2:1361445/147‑1 (MQ=255)
                                                                                           gCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTACGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGaa                                                            <  1:3682858/150‑1 (MQ=255)
                                                                                              ggggggggGTGTAATTTCTTTTTAAATGGGTCTTGGGTCGCCTCAGTGGCCTGCCTGCGCGGTCGGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAcgc                                                         <  1:4864910/150‑1 (MQ=255)
                                                                                              gggCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAcgc                                                         >  2:4130141/1‑150 (MQ=255)
                                                                                               ggCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAcgcg                                                        >  1:6961078/1‑150 (MQ=255)
                                                                                                      cTGTAATTTCTTTTCAAAGGGGCGAGGGGCCGCCTCCGGGGCGTGCCGGCGGGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGgc                                                 <  1:4763814/150‑1 (MQ=255)
                                                                                                         gAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTTGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCACGAACGCGAGCAGGCgct                                              >  1:1831496/1‑150 (MQ=255)
                                                                                                          aatttgttttCAAACGGGTCTGGGGTCGCCTCCGGGGCCTGCCTGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctg                                             <  2:2442528/143‑1 (MQ=255)
                                                                                                          aaaTTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctg                                             >  2:6329851/1‑150 (MQ=255)
                                                                                                                gTTTGAAACGGGTCATGGGTCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCtggctggc                                       <  2:2135657/150‑1 (MQ=255)
                                                                                                                 tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                      >  2:7161810/1‑150 (MQ=255)
                                                                                                                 tttCAAAAGGGTCATGGGTCGCATCACTGCCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                      <  2:7601318/150‑1 (MQ=255)
                                                                                                                  ttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGccc                                     <  2:8470947/150‑1 (MQ=255)
                                                                                                                     aaaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACCGCGAAGTTGTTGCGTTTGCAGAGCTTTGCTTCGACCGTCTGGGTGCCGCTTTACGGCCCCTGTACCTCGAAAGCGAGCAGGCGCAGGCTGGCCCAGt                                  >  1:2488447/1‑150 (MQ=255)
                                                                                                                         ggtgcttGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                              <  2:4799756/144‑1 (MQ=255)
                                                                                                                         gggcctggggCCGCCCCCCGGGCCTGCCTGGGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                              <  1:4029531/143‑1 (MQ=255)
                                                                                                                         gggTCTGGGGTCGCCTCCGGGGCCGGCCGGCGCGGTGGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                              <  1:1211390/150‑1 (MQ=255)
                                                                                                                            tCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTc                           >  2:8583340/1‑150 (MQ=255)
                                                                                                                                 ggTCGCCTCACGGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                      <  2:8223468/150‑1 (MQ=255)
                                                                                                                                  gTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgcg                     >  1:7208004/1‑150 (MQ=255)
                                                                                                                                  gTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgcg                     >  1:8524067/1‑150 (MQ=255)
                                                                                                                                     gCCTCCCGGGCCGGCCGGCGGGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCAt                  <  2:1124571/150‑1 (MQ=255)
                                                                                                                                       cTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCt                >  2:4366436/1‑150 (MQ=255)
                                                                                                                                          tcgggCCTGCGGGCGGGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGa             <  1:6399965/147‑1 (MQ=255)
                                                                                                                                          ccTGGCGTGCCTGGGGGGTCTGCTCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGa             <  1:7477243/149‑1 (MQ=255)
                                                                                                                                          aCTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCACATCTGGa             <  1:8727619/150‑1 (MQ=255)
                                                                                                                                          aCTGGCCTGCCTGCGCGGGCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCATCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCACGTCTTCGCGCAAATGGa             >  1:1775651/1‑150 (MQ=255)
                                                                                                                                           cTGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa            <  2:2805979/150‑1 (MQ=255)
                                                                                                                                              ccctgcctgCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGAc         >  1:7004185/2‑150 (MQ=255)
                                                                                                                                                gtgcctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa       >  2:3735592/2‑150 (MQ=255)
                                                                                                                                                gtgcctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa       >  2:7233911/2‑150 (MQ=255)
                                                                                                                                                ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa       >  1:3851494/1‑150 (MQ=255)
                                                                                                                                                ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa       >  1:4671435/1‑150 (MQ=255)
                                                                                                                                                ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa       >  1:5309763/1‑150 (MQ=255)
                                                                                                                                                   cctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGc    >  2:3450082/1‑150 (MQ=255)
                                                                                                                                                    gtgcGCGTTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCg   <  1:2845702/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGG  >  NC_002947/1468998‑1469296

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: