Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,464,207 T→G 11.9% 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,2070TG11.9% 358.9 / 9.3 166Q473P (CAG→CCG) algXalginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (74/55);  new base G (2/21);  total (82/84)
Fisher's exact test for biased strand distribution p-value = 1.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

TTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGGTG  >  NC_002947/1464058‑1464356
                                                                                                                                                      |                                                                                                                                                     
ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                                       >  2:1903009/1‑150 (MQ=255)
  gAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                                     >  1:1360342/1‑150 (MQ=255)
     cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGTGCGGTCTGCGCa                                                                                                                                                  <  1:3129185/150‑1 (MQ=255)
     cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                                  >  1:3704092/1‑150 (MQ=255)
      ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                                 <  1:1194633/150‑1 (MQ=255)
      ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                                 >  1:28857/1‑150 (MQ=255)
      ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                                 <  1:2440636/150‑1 (MQ=255)
       tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                                >  2:1632410/1‑150 (MQ=255)
        ttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGc                                                                                                                                               >  1:3494000/1‑150 (MQ=255)
         tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCa                                                                                                                                              <  1:3615222/150‑1 (MQ=255)
          cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                             <  2:3958398/150‑2 (MQ=255)
          cGAGCCCCTGGGAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                             <  1:3253246/150‑2 (MQ=255)
           gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                            >  2:4512365/1‑150 (MQ=255)
            aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACACGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAATGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                           <  2:2203089/150‑1 (MQ=255)
            aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                           >  2:3198123/1‑150 (MQ=255)
            aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                           >  2:3387337/1‑150 (MQ=255)
             tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACCAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCaggcc                                                                                                                                          <  2:4273746/150‑5 (MQ=255)
             tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGc                                                                                                                                          >  1:2282346/1‑150 (MQ=255)
             tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCaggcc                                                                                                                                          <  2:377768/150‑5 (MQ=255)
             tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGc                                                                                                                                          >  1:1608796/1‑150 (MQ=255)
               cccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGa                                                                                                                                        >  2:256236/1‑150 (MQ=255)
               cccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATGCGTTGGAAATGGGTGATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCGGCGCAGGCACGGCgc                                                                                                                                        >  2:2822098/1‑149 (MQ=255)
                ccTGTTAGTAAACCTGCGGCGGTACCATCTCGGCGTGCGCGTCAACACCAAAAAATTCAAACGCGATCAGGGAGGTGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGACTCACTGGCCTGCATGAGCGGTCTGCTCAGTCACGGCGaa                                                                                                                                       >  2:772830/1‑150 (MQ=255)
                ccTGGTAGTACCCCTGGGGCGGTACCAGCGCGGGGTGCGCGGCCCCACCAAACACGGCCACCGCGAGCAGGGGGGGGCTGTAATTTGGTTTGAAACGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGGGGTCTGCGCAGGCACGGCGaa                                                                                                                                       <  2:2787562/150‑1 (MQ=255)
                ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                       >  2:1871423/1‑150 (MQ=255)
                   gttattaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGCCACCACCAAACACGGCCACCGCGAGCGGGGGGGGGCTGGAAATTCTTTTTAAACGGGTCTTGGGTCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAgtt                                                                                                                                    <  2:1774963/145‑1 (MQ=255)
                      agtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAACCACGGCCAACGCGAGCTGGGCGGGGCTGTAAATTCGTTTGAAAAGGGTCTTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgtt                                                                                                                                 <  2:3343073/150‑1 (MQ=255)
                      agtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACACGGCCAACGCGAGCCGGGCGGGGCTGGAAATTCTTTTTAAACGGGTAATGGGTCGCCTCACTGGCCTGCCTGCGCGTTCTGCGCAGGCACGGCGAAgttgtt                                                                                                                                 <  2:4431315/150‑1 (MQ=255)
                          aCCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGt                                                                                                                             >  2:2797205/1‑150 (MQ=255)
                             cTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTGGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                          >  1:2332733/1‑150 (MQ=255)
                               gcggcgGTACCAGCGGGGCGTGCGGGGCACCACCAAACACGGCCACCGCGAGCAGGGTGGGGGTGTAAATTGTTTTGAAATGGTTCATGGGTCGCCTCACTGGCCTGCCTGGGGGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCa                                                                                                                        <  2:2651832/150‑1 (MQ=255)
                               gcggcgGAACCAGCGCGGCGTGCCCGCCACAACAAAACACCCCAAACCCCAGCAGGGCGGGGATGAAAATTCTTTTAAAAAGGCTCATGGGACCCAAAACTGCCCCGCCTGCGCGCTCTCCCAAGCCACGGCGAAGTTGTTGCGTTTGCa                                                                                                                        <  2:905644/150‑1 (MQ=255)
                                     gTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCtt                                                                                                                  >  2:2322716/1‑150 (MQ=255)
                                         cAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCt                                                                                                              >  1:788308/1‑150 (MQ=255)
                                            cgcgGCGTGCGCGGCACCACCAAACACGCCCAACGCGAGCTGGGGGGGGGTGGAAATTCGTTTGAACTGGGTATGGGGGCGCCTCACGGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCg                                                                                                           <  1:3731448/150‑1 (MQ=255)
                                             gcggcgTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGa                                                                                                          >  2:663048/1‑150 (MQ=255)
                                                 cgTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCtt                                                                                                      >  2:374360/1‑150 (MQ=255)
                                                  gTGCGCGCCACCACCACACAAGGCCACCGCGAGCCGGGGGGGGGTGGAATTTCGTTTTAAATGGTTCTTGGGTCGCCCCACGGGCCTGCCTGCGCGGCCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTc                                                                                                     <  2:3492221/150‑1 (MQ=255)
                                                   tGCGCGGCCCCACCAAACACGGCCACCGCGAGCCGGGGGGGGCTGTAATTTCTTTTGAAACGGGTCTTGGGTCGCCTCACGGGCCGGCCTGCGCGGTGTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCt                                                                                                    <  2:1501144/150‑1 (MQ=255)
                                                   tGCGCGGCACCACCAAACACGGCAAACGCGAGCAGGGCGGGGCTGTAAATTCTTTTGAACAGGGTCATGGGTCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCt                                                                                                    <  2:3695201/150‑1 (MQ=255)
                                                   ggcgcgcCCCCACAAAACACGCCCACCGCGAGCCGGGGGGGGGTGTAAATTCGTTTCAAACGGGTCTGGGGTCGCCCCACGGGCCGGCGGGCGGGGACGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCt                                                                                                    <  2:1501075/149‑1 (MQ=255)
                                                    gcgcGGCACCACCAAACAAGGCAAACGCGTGCCGGGGGGGGCTGTAATTTCTTTTTAAAAGGGTCTTGGGTCGCCTCACTGGCCGGCCTGGGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                                   <  2:497539/150‑1 (MQ=255)
                                                      gcgGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGGCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTggg                                                                                                 >  1:50820/1‑150 (MQ=255)
                                                          caccacCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGGGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGcc                                                                                             >  1:531921/1‑150 (MQ=255)
                                                             caccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCt                                                                                          >  2:1765870/1‑150 (MQ=255)
                                                              accaAACCCGCCCACCGCGAGCCGGGGGGGGCTGGAATTTGGTTTTAAACGGGTCATGGGGCGCCCCACGGGCCGGCCGGCGCGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                                         <  2:2442659/150‑1 (MQ=255)
                                                               ccaAACACGCCCACGGCGCGCGGGGGGGGGGTGCAATTTCTTTTTAAATGGGAGTTGGGTCGCCACCCGGCCGTGCCTGGGGGGTCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCttt                                                                                        <  1:4190996/150‑1 (MQ=255)
                                                               ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCGTCGACCGTCTGGGTGCCGCttt                                                                                        >  1:746183/1‑150 (MQ=255)
                                                                 aaaCAAGCCAAACGCGAGCCGGGCGGGGCTGTAAATTCTTTTTAAAAGGGTAATGGGTCGCCTCACTGGCCTGCCTGCGCTTTCTGCTCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTcc                                                                                      <  1:674170/150‑1 (MQ=255)
                                                                      cggCCACGGCGAGCTGGGGGGGGGTGTAATTTGTTTTTAAATGGGTCTTGGGTCGCCTCAGGGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                                 <  1:2553448/149‑1 (MQ=255)
                                                                      cggCCACCGCGCGCCGGGGGGGGGTGTAATTTCTTTTGAAACGGGCCCTGGGTCGCCTCACGGGCGGGCCTGCGGGGGGGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                                 <  1:2162560/149‑1 (MQ=255)
                                                                      cggCCACCGCGAGCTGGGGGGGGCTGTAATTTGTTTTTAACTGGGTATTGGGTCGCCTCAGTGGCCTGCCTGGGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                                 <  2:4287782/149‑1 (MQ=255)
                                                                       ggCCACCGCGCGCCGGGGGGGGCTGTATATTCGTTTTAACTGGGTCCTGGGGCGCCTCCCGGGCCGGCGTGCGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGcccc                                                                                <  2:1224450/150‑1 (MQ=255)
                                                                       cgCCAACGGCGCGCCGGGGGGGGCTGTATATTTTTTTTAATTGGTTCTGGGGTCGCCAAAGTGGCCTGGCGGGGGGGTGTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                                 <  1:2050853/149‑1 (MQ=255)
                                                                             cgcgAGCTGGGGGGGGCTGGAAATTCGTTTGAAATGGGTCTTGGGTCGCCCCACTGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccccagaagg                                                                          <  2:3836391/150‑7 (MQ=255)
                                                                                cgcgcgggggggggTTGTAATTTGTTTTCAAATGGGTCTTGGGTCGCCTCAGTGGCCTGCCTGGGGGGTCTGCGCTGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                                        <  1:3241886/144‑1 (MQ=255)
                                                                                cgcgccgggcggggcgGTAAATTCGTTTGAAAAGGGTCTGGGGTCGCCTCACTGGCCGGCCGGCGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                                        <  2:2389016/144‑1 (MQ=255)
                                                                                cgAGCAGGGCGGGGCTGGAAATTCGTTTGAAACGGGTCTTGGGTCGCCTCACTGGCCGGCCTGCGCGTTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                                        <  2:2195428/150‑1 (MQ=255)
                                                                                   gccggggggggCTGCAATTTCTTTTGAAATGGGTCTTGGGTCGCCTCACTGGCCGGCGGGGGGGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGaa                                                                     <  2:2413359/147‑1 (MQ=255)
                                                                                      ggggggggTTGTAATTTCGTTTGAAAGGGGTCTTGGGCCGCCCCCCGGCCGGGCCTGGGGGGTCGGCGCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAcgc                                                                  <  1:264750/150‑1 (MQ=255)
                                                                                        ggggggTTGGAAATTCTTTTGAAATGGTTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                                <  2:1845893/148‑1 (MQ=255)
                                                                                        ggggggGTGTAATTTCGTTTGAAAGGGGTCTGGGGGGGCCGCTGGGGCGGGCGGGGGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                                <  1:4307726/148‑1 (MQ=255)
                                                                                        ggggggGTGGATATTCGTTTTAAAAGGGTCTTGGGTCGCCTCCCTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                                <  2:1977590/148‑1 (MQ=255)
                                                                                        ggggggCTGTAATTTGGTTTGAAACGGGCCATGGGCCGCCCCCCGGGCCGGCGTGCGCGGTCGCCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                                <  2:116532/148‑1 (MQ=255)
                                                                                          ggggCTGCAATTTTTTTTTAAATGGGCCTGGGGCCGCCCCCCGGGCCGGCGGGGGCGGGCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGc                                                              <  1:1235621/150‑1 (MQ=255)
                                                                                              cTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAACTTGTTGCGTTTGCAGAGCTTTGCTTAGACCATCATGGTGCCGCATTACGGCCCCAGTACCTACAACCCGAGCAGgc                                                          >  1:1919273/1‑150 (MQ=255)
                                                                                                ttaatttCGTTTGAAACGGGACATGGGTCGCCTACCTGGCCTGCCTGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGgcgc                                                        <  2:2279695/145‑1 (MQ=255)
                                                                                                 gAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgct                                                       >  2:3135051/1‑150 (MQ=255)
                                                                                                   aaTTCTTTTGAAACGGGTCATGGGCCGCCTCACGGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctgg                                                     <  2:3085600/150‑1 (MQ=255)
                                                                                                     ttCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctggct                                                   >  2:3288324/1‑150 (MQ=255)
                                                                                                       gttttCAAACGGGTCTGGGGCCGCCTCCGGGGCCGGCGGGCGGGGTGGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGctggctgg                                                 <  1:38743/148‑1 (MQ=255)
                                                                                                       cGTTTGAAACGGGTCATGGGGCGCCTCCCGGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGctggctgg                                                 <  2:662068/150‑1 (MQ=255)
                                                                                                        gTTTGAAACGGGACTGGGGTCGCCTCCCGGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTGGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAATGCGAGCAGGCGCtggctggc                                                <  2:2509690/150‑1 (MQ=255)
                                                                                                         tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                               >  2:4055926/1‑150 (MQ=255)
                                                                                                         tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACCGCGAAGTTGTTTCGTTTGCAGAGCTTGGCTTCGACCATCTGGGTGCCGCTTTACGGCCCCAGTAACACGAACGCGAGCAGGCGCAGGCTGGcc                                               >  1:1919542/1‑150 (MQ=255)
                                                                                                           tGAAACGGGTCTTGGGTCGCCCCCCGGGCCTGCCTGCGGGGGCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa                                             <  1:2713457/150‑1 (MQ=255)
                                                                                                           tGAAAAGGGTCTGGGGTCGCCCCACTGGCCGGCGTGCGCGGTCTGCGCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa                                             <  2:2282346/150‑1 (MQ=255)
                                                                                                             aaaaGGGTCATGGGTCGCCTCACTGGCCTGCCTTCGCGGTCAGCGCACGCACCGCGAAGTTGTTTCGTTTGCAGAGATTGGCTTAGACCTTCTGGGTGCCGCTTTACGGCCCCAGTACCACGAAAGCGAGCACGCGCAGGCTGGCCCAGt                                           >  1:2367961/1‑150 (MQ=255)
                                                                                                                 ggggctggggTCGCCTCACGGGCCTGGGGGGGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                       <  1:3462516/143‑1 (MQ=255)
                                                                                                                 gggTCTTGGGGCGCCTCACTGGCCTGCCTGCGCGGCCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                       <  2:1195529/150‑1 (MQ=255)
                                                                                                                   gTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                                     >  1:1017259/1‑150 (MQ=255)
                                                                                                                   gTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                                     >  1:1017238/1‑150 (MQ=255)
                                                                                                                       ggggTCGCCTCCCGGGCCTGCCTGCGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCGTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTc                                 <  2:2101685/149‑1 (MQ=255)
                                                                                                                         ggTCGCCTCACTGGCCTGCCGGCGCGGACTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                               <  2:4237929/150‑1 (MQ=255)
                                                                                                                           tCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCgcgc                             >  2:2752921/1‑150 (MQ=255)
                                                                                                                                  agtGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGa                      <  2:4358627/148‑1 (MQ=255)
                                                                                                                                    tGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAg                    >  2:4021672/1‑150 (MQ=255)
                                                                                                                                     gccgtgcctgcGCAGACCGCGCAGGCAGGCGAAGTTGTTGCGTTTGCAGAGCGTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCAGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCACCGCTTCGCGCATATGGAAGAc                  >  2:3318436/5‑150 (MQ=255)
                                                                                                                                        gtgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa                >  2:1022122/2‑150 (MQ=255)
                                                                                                                                           cctgcGCAGACCGCGCAGGCAGGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGc             >  2:3024454/1‑150 (MQ=255)
                                                                                                                                           cctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAg                    <  1:3083619/143‑1 (MQ=255)
                                                                                                                                            ctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCg            <  2:3617861/150‑1 (MQ=255)
                                                                                                                                            ctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCg            >  2:1358522/1‑150 (MQ=255)
                                                                                                                                                  ggTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTc      <  2:1663431/150‑1 (MQ=255)
                                                                                                                                                  ggTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTc      <  1:2781197/150‑1 (MQ=255)
                                                                                                                                                   cccTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCg     >  2:3543281/3‑150 (MQ=255)
                                                                                                                                                    tCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCgg    >  1:2452242/1‑150 (MQ=255)
                                                                                                                                                     cTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGgt   >  1:4281108/1‑150 (MQ=255)
                                                                                                                                                     cTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGgt   >  1:4281317/1‑150 (MQ=255)
                                                                                                                                                      tGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGgtg  >  2:2016288/1‑150 (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: