Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,464,207 T→G 10.0% Q473P (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,464,2070TG10.0% 535.5 / 6.4 232Q473P (CAG→CCG) algXalginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (119/75);  new base G (2/27);  total (123/109)
Fisher's exact test for biased strand distribution p-value = 1.62e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 234 total aligned reads displayed.

TTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGGTG  >  NC_002947/1464058‑1464356
                                                                                                                                                      |                                                                                                                                                     
ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                                       >  2:1666052/1‑150 (MQ=255)
 tGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTg                                                                                                                                                      <  1:8367709/150‑1 (MQ=255)
 tGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTg                                                                                                                                                      >  2:8218640/1‑150 (MQ=255)
  gAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                                     >  2:2948192/1‑150 (MQ=255)
    gCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcgc                                                                                                                                                   >  1:3203254/1‑150 (MQ=255)
    gCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGGCTgcgc                                                                                                                                                   >  1:5589380/1‑150 (MQ=255)
     cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                                  <  1:473056/150‑1 (MQ=255)
      ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                                 >  1:6322516/1‑150 (MQ=255)
      ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                                 >  2:5888444/1‑150 (MQ=255)
      ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                                 >  2:5513867/1‑150 (MQ=255)
      ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                                 >  2:7990009/1‑150 (MQ=255)
      ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                                 <  2:3722783/150‑1 (MQ=255)
       tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAgg                                                                                                                                                <  1:100779/150‑1 (MQ=255)
        ttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCGTGCCTGCGCAGACCGCGCAgg                                                                                                                                                <  1:4035470/150‑1 (MQ=255)
         tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCa                                                                                                                                              <  1:2680314/150‑1 (MQ=255)
          cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCAc                                                                                                                                             <  1:6003289/150‑1 (MQ=255)
          cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCAc                                                                                                                                             <  1:2110555/150‑1 (MQ=255)
           gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                            >  1:66538/1‑150 (MQ=255)
            aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCGTGCCTGCGCAGACCGCGCAGGCaggc                                                                                                                                           <  1:7645052/150‑4 (MQ=255)
             tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGc                                                                                                                                          >  1:5076148/1‑150 (MQ=255)
             tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGc                                                                                                                                          >  2:1089659/1‑150 (MQ=255)
               cccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTGTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCCGGCACGGCGa                                                                                                                                        >  2:2187072/1‑150 (MQ=255)
                ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAATGGGTCATGGGTCGCCTCACTGGCCTGCCAGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                       >  1:1227689/1‑150 (MQ=255)
                 cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAg                                                                                                                                      >  2:2661286/1‑150 (MQ=255)
                 cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAg                                                                                                                                      >  2:6117589/1‑150 (MQ=255)
                  tGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgt                                                                                                                                     >  2:8070002/1‑150 (MQ=255)
                   tttataaaCCCTCCGGCGTTACCAGCGCGGCGTGCGCGCAACAACAAAAAAAGCAAAACGCAAGCAGGGCGGGGATGAAAATTCTTTTTAAAAGGGTCATGGGTCCCCTAACTGGCCTCCCTGCGCGCTCTGCCCCGGCACGGCGAAgtt                                                                                                                                    <  2:156702/144‑1 (MQ=255)
                    gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttg                                                                                                                                   >  2:8046706/1‑150 (MQ=255)
                     tagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGACgttgt                                                                                                                                  >  2:696875/1‑150 (MQ=255)
                     tagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgt                                                                                                                                  >  1:4527047/1‑150 (MQ=255)
                        taACCCTGCGGCGGTACCCGCGCGGCGTGCGGGGCACCACCACCCACGGCCAACGCGAGCCGGGGGGGGCTGGAATTTCGTTGGAAAAGGGTCATGGGTCGCCCCACGGGCCTGCCTGGGGGGTGGGCGCAGGCACGGCGAAGTTGTTGc                                                                                                                               <  2:5106065/150‑1 (MQ=255)
                          aCCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGt                                                                                                                             >  1:1838842/1‑150 (MQ=255)
                            cgggcggcgGAACCCGCGCGGCGTGCGCGGCACCACCAACCACGGCCAACGCGAGCTGGGGGGGGCTGTAATTTCTTTTGAAAAGGGTCTGGGGTCGCCTCAGGGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                           <  2:3503460/147‑1 (MQ=255)
                            ccggcggcgGTACCCGCGCGGGGGGCGCGCCCCCACAAAACCCGGCCACCGCGCGCGGGGGGGGGCTGCAATTGCTTTTCAAAGGGGTCTTGGGTCGCCGCCCGGGCCGGCGTGGGCGGTGGGCGCTGGCACGGCGAAGTTGTTGCGttt                                                                                                                           <  2:3298636/147‑1 (MQ=255)
                            ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                           >  1:2793169/1‑150 (MQ=255)
                             cTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCAAGGCGAAGTTGTTGCGTTTg                                                                                                                          >  1:1309522/1‑150 (MQ=255)
                              tGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGc                                                                                                                         >  2:7082716/1‑150 (MQ=255)
                                        ccAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGc                                                                                                               >  1:3779981/1‑150 (MQ=255)
                                        ccAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGc                                                                                                               >  2:8937369/1‑150 (MQ=255)
                                         cAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCATAGCTTGGCt                                                                                                              >  1:2184702/1‑150 (MQ=255)
                                             gcggggGGCGCGGCACCCCCAACCCCGGCCACCGGGCGCGGGGGGGGGGTGGAATTTGGTTTGAACCGGGTCCTGGGGCGCCTCCGGGGCCGGGGGGGGGGGTGGGCGCGGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGa                                                                                                          <  2:2563756/150‑1 (MQ=255)
                                                  gTGCGCGGCACCACCAACCAAGGCCACCGCGAGCTGGGCGGGGCTGGAAATTCTTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTc                                                                                                     <  2:9066148/150‑1 (MQ=255)
                                                   ggcgcgcCCCCACCAAACACGGCAACGGGGAGGTGGGGGGGGCTGTAATTTGTTTTGAATGGGTTCTGGGGTCGCCACCCGGGCCGGCGGGGGCGGTGTGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCt                                                                                                    <  1:6261845/149‑1 (MQ=255)
                                                    gcgcGGCCCCACCAAACCCGGCCACCGCGCGCCGGGGGGGGCTGGAATTTGTTTTCAAACGGGTCCTGGGTCGCCTCCGGGGCCTGCCTGCGGGGCCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                                   <  2:3164546/150‑1 (MQ=255)
                                                         gCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCTTTTGAAAAGGGTCTTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGc                                                                                              <  2:8168505/150‑1 (MQ=255)
                                                           accaccAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTATGGGTGCCg                                                                                            >  1:7786989/1‑150 (MQ=255)
                                                              accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                                         >  1:6230127/1‑150 (MQ=255)
                                                                caAACAAGGCCAACGCGAGCAGGGCGGGGCTGTAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTTCAGAGCTTGGCTTAGACCTTCTGGGTGCCGCTTTc                                                                                       >  2:3907130/1‑150 (MQ=255)
                                                                  accacgcccACCGCGCGCGGGGGGGGGCTGTAATTTCGTTTTAAATGGGTCTTGGGCCGCCTCACGGGCGTGCGTGGGGGGTCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCg                                                                                     <  1:7888882/143‑1 (MQ=255)
                                                                  aaCAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCg                                                                                     >  1:3760465/1‑150 (MQ=255)
                                                                      cggCCAAGGCGAGCAGGGCGGGGCTGGAAATTGTTTTGAAACGGGTCTTGGGGCGCCTCCCTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                                 <  2:5570456/149‑1 (MQ=255)
                                                                      cgcccAACGCGAGCTGGGCGGGGCTGGAAATTCGTTTTAAATGGGTCTTGGGTCGCCCCACTGGCCGGCCTGCGGGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                                 <  2:2397509/147‑1 (MQ=255)
                                                                          caccgcgAGCGGGGGGGGGGTTTAATTTCGTTTCAAACGGGTATGGGGGCGCCTCACTGGCGGGCCTGCGGGGTCGGCGCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGt                                                                             <  1:4757742/147‑1 (MQ=255)
                                                                          cAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACCGCGAAGTTGTTTCGTTTGCAGAGCTTGGCTTAGACCATATGGGTGCCGCATTACGGCCCCAGt                                                                             >  1:8518892/1‑150 (MQ=255)
                                                                           accgcgcgCCGGGCGGGGGTGGAATTGGGTTGTAAACGGGTCTGGGGTCGCCTCCCTGGCCTGCGTGCGGGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccccagaa                                                                            <  2:2419470/148‑5 (MQ=255)
                                                                           accgcgAGCTGGGGGGGGCTGTAAATTCGTTTTAAAAGGGTCTTGGGTCGCCACACTGGCCTGCCGGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTa                                                                            <  1:7408653/148‑1 (MQ=255)
                                                                            cggcgcgcCGGGGGGGGGTGTAAATTCGTTTGAAACGGGTCTGGGGTCGCCTCACTGGCCGGCCGGGGGGGTGTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTAc                                                                           <  2:7814017/148‑1 (MQ=255)
                                                                            aCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTAc                                                                           >  2:4255706/1‑150 (MQ=255)
                                                                             cgcgcgCTGGGCGGGGCGGGAAATTCGTTTCAAATGGGGCGGGGGGCGCCTCCGGGGCGGGCGGGCGCGGGCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTAcc                                                                          <  1:5740407/150‑1 (MQ=255)
                                                                               cgAGCTGGGCGGGGCTGTAAATTCGTTTGAAATGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                                        <  2:5339986/150‑1 (MQ=255)
                                                                               cgAGCTGGGCGGGGCTGGAAATTCGTTTGAAATGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                                        <  2:8362191/150‑1 (MQ=255)
                                                                                gAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCg                                                                       >  2:44431/1‑150 (MQ=255)
                                                                                gAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCAGCGCGGTCTGCGCAGTCACGGCGAAGTTGTTTCGTTTTCAGAGCTTGGCTTCGACCATCTGGGTGCCGCTTTCCGGCCCCGGTACCTCg                                                                       >  2:44709/1‑150 (MQ=255)
                                                                                     gggCGGGGCTGGAAATTCGTTTGAAAAGGGTCAGGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAcgc                                                                  <  2:2981888/150‑1 (MQ=255)
                                                                                            gCTGTAAATTCGTTTTAATTGGTTCTGGGGTCGCCTCACTGGCCTGCGTGCGCGGTCGGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAgg                                                           <  2:5005634/150‑1 (MQ=255)
                                                                                              tGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGgcg                                                         >  2:1301758/1‑150 (MQ=255)
                                                                                                gAATTTCGTTTTAAATGGGTCTTGGGTCGCCTCACTGGCCGGCCTGCGCGGGCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgct                                                       <  2:2056436/150‑1 (MQ=255)
                                                                                                     tatttttAAATGGGTCATGGGTCCCAACCCTGGCCTGCCTGCGCGTTAGGCCCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctggctg                                                  <  1:4834696/147‑1 (MQ=255)
                                                                                                     tCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctggctg                                                  >  1:2915820/1‑150 (MQ=255)
                                                                                                        tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                               >  1:2960449/1‑150 (MQ=255)
                                                                                                        tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                               >  2:5488389/1‑150 (MQ=255)
                                                                                                        tttAAAATGGTTCCTGGGCCGCCCCACTGGCGTGCGTGCGGGGTCTGCGCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                               <  1:4900627/150‑1 (MQ=255)
                                                                                                           gAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTTCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCACGAACGCGAGCAGGCGCTGGCTGGCCCAg                                            >  1:531886/1‑150 (MQ=255)
                                                                                                           gAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCAGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAg                                            >  1:6640203/1‑150 (MQ=255)
                                                                                                            aaaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTTCAGAGCTTTGCTTCGACCTTATGGGTGCCGCTTTACGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCAGGCCCact                                           >  1:8803906/1‑148 (MQ=255)
                                                                                                            aaaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGt                                           >  2:8693770/1‑150 (MQ=255)
                                                                                                             aacgggTCTGGGGTCGCCTCACGGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTc                                          <  2:6907377/147‑1 (MQ=255)
                                                                                                             aaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTAGACATTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTc                                          >  2:923310/1‑150 (MQ=255)
                                                                                                              aaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTcc                                         >  2:7572896/1‑150 (MQ=255)
                                                                                                              aaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGGTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTcc                                         >  2:7455187/1‑150 (MQ=255)
                                                                                                                gggTCTTGGGTCGCCGCCGGGGCCGGCGTGGGGGGTCTGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                       <  1:2549248/150‑1 (MQ=255)
                                                                                                                gggTCATGGGTCGCCTCCCGGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                       <  2:6757159/150‑1 (MQ=255)
                                                                                                                   tcttGGGTCGCCTCCCGGGCGTGCCTGCGCGTTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTc                                    <  1:2414426/147‑1 (MQ=255)
                                                                                                                   tCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTTTTGCGTTTGCAGAGCTTTGCTTAGACCGTCTGGGTGCCGCTTTCCGGCCCCAGTACCACGAACGCGAGCAGGCGCAGGCTGGCCCAGTCCGCGTc                                    >  1:6999309/1‑150 (MQ=255)
                                                                                                                   tCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTc                                    >  2:4694848/1‑150 (MQ=255)
                                                                                                                     aTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCtt                                  >  2:6583486/1‑150 (MQ=255)
                                                                                                                        ggTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                               <  2:5320468/150‑1 (MQ=255)
                                                                                                                          tCGCCTCAGTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCgcgc                             <  1:2948192/150‑1 (MQ=255)
                                                                                                                                  ggggCCTGCGGGCGGGGCCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa                     <  1:2661286/148‑1 (MQ=255)
                                                                                                                                   tGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAg                    >  2:6727026/1‑150 (MQ=255)
                                                                                                                                    ggCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGa                   >  2:8621923/1‑150 (MQ=255)
                                                                                                                                       ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa                >  1:8103158/1‑150 (MQ=255)
                                                                                                                                            ctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCg            >  1:5347522/1‑150 (MQ=255)
                                                                                                                                             tgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCgg           >  1:1766561/1‑150 (MQ=255)
                                                                                                                                              gcgcGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGc          <  2:3663455/150‑1 (MQ=255)
                                                                                                                                              gcgcGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGc          <  2:2702235/150‑1 (MQ=255)
                                                                                                                                                gcgGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCg        <  2:5791274/150‑1 (MQ=255)
                                                                                                                                                gcgGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCg        >  1:8654107/1‑150 (MQ=255)
                                                                                                                                                   gTCTGCGCAGGCGCGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCg     <  2:5352593/150‑1 (MQ=255)
                                                                                                                                                   gTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCg     <  1:7082716/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGGTG  >  NC_002947/1464058‑1464356

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: