Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,464,205 T→A 10.4% T474S (ACC→TCC) ‡ algX ← alginate biosynthesis protein AlgX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,464,2050TA10.4% 356.0 / 28.4 166T474S (ACC→TCC) ‡algXalginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (74/66);  new base A (5/10);  total (81/85)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-02

Only 100 of 168 total aligned reads displayed.

TCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGG  >  NC_002947/1464056‑1464354
                                                                                                                                                     |                                                                                                                                                     
tCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGt                                                                                                                                                       <  1:1304414/150‑1 (MQ=255)
  ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                                     <  2:4283214/150‑1 (MQ=255)
  ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                                     >  2:1903009/1‑150 (MQ=255)
    gAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                                   >  2:3672130/1‑150 (MQ=255)
       cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGTGCGGTCTGCGCa                                                                                                                                                <  1:3129185/150‑1 (MQ=255)
       cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                                >  1:3704092/1‑150 (MQ=255)
        ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                               >  1:28857/1‑150 (MQ=255)
        ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                               <  1:1194633/150‑1 (MQ=255)
        ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                               <  1:2440636/150‑1 (MQ=255)
         tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                              >  2:1632410/1‑150 (MQ=255)
         tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                              <  1:2356818/150‑1 (MQ=255)
         tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                              >  2:1632891/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)
              aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCACCGCGCGCAGGGGGGGGCTGGAAATTCTTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                         <  2:2203083/150‑1 (MQ=255)
              aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                         >  2:3387337/1‑150 (MQ=255)
              aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                         >  2:3198123/1‑150 (MQ=255)
               tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGc                                                                                                                                        >  1:2282346/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)
                   cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAg                                                                                                                                    >  2:2304494/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)
                              ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGACATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                         >  2:3462516/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)
                                                     tGCGCGGCCCCACCAAACACGGCCACCGCGAGCCGGGGGGGGCTGTAATTTCTTTTGAAACGGGTCTTGGGTCGCCTCACGGGCCGGCCTGCGCGGTGTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCt                                                                                                  <  2:1501144/150‑1 (MQ=255)
                                                     tGCGCGGCACCACCAAACACGGCAAACGCGAGCAGGGCGGGGCTGTAAATTCTTTTGAACAGGGTCATGGGTCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCt                                                                                                  <  2:3695201/150‑1 (MQ=255)
                                                      gcgcGGCACCACCAAACAAGGCAAACGCGTGCCGGGGGGGGCTGTAATTTCTTTTTAAAAGGGTCTTGGGTCGCCTCACTGGCCGGCCTGGGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                                 <  2:497539/150‑1 (MQ=255)
                                                        gcgGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGGCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTggg                                                                                               >  1:50820/1‑150 (MQ=255)
                                                           gCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTTAAAAGGGTCAGGGGTCGCCTAAGTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGc                                                                                            <  2:3755701/150‑1 (MQ=255)
                                                            caccacCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGGGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGcc                                                                                           >  1:531921/1‑150 (MQ=255)
                                                               caccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCGTGGCTTCGACCTTCTGGGTGCCGCt                                                                                        >  1:1545808/1‑150 (MQ=255)
                                                                accaAACCCGCCCACCGCGAGCCGGGGGGGGCTGGAATTTGGTTTTAAACGGGTCATGGGGCGCCCCACGGGCCGGCCGGCGCGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                                       <  2:2442659/150‑1 (MQ=255)
                                                                 ccaAACACGCCCACGGCGCGCGGGGGGGGGGTGCAATTTCTTTTTAAATGGGAGTTGGGTCGCCACCCGGCCGTGCCTGGGGGGTCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCttt                                                                                      <  1:4190996/150‑1 (MQ=255)
                                                                 ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCttt                                                                                      >  2:4028031/1‑150 (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)
                                                                          gCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTACGGCCCCt                                                                             >  1:4276181/1‑150 (MQ=255)
                                                                             aaCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTa                                                                          >  2:3782199/1‑150 (MQ=255)
                                                                               cgcgAGCTGGGGGGGGCTGGAAATTCGTTTGAAATGGGTCTTGGGTCGCCCCACTGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccccagaagg                                                                        <  2:3836391/150‑7 (MQ=255)
                                                                                 cgcgcgggggggggTTGTAATTTGTTTTCAAATGGGTCTTGGGTCGCCTCAGTGGCCTGCCTGGGGGGTCTGCGCTGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                                      <  1:3241886/144‑1 (MQ=255)
                                                                                 cgcgccgggcggggcgGTAAATTCGTTTGAAAAGGGTCTGGGGTCGCCTCACTGGCCGGCCGGCGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                                      <  2:2389016/144‑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)
                                                                                                      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)
                                                                                                            tGAAACGGGTCTTGGGTCGCCCCCCGGGCCTGCCTGCGGGGGCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa                                           <  1:2713457/150‑1 (MQ=255)
                                                                                                              aaaaGGGTCATGGGTCGCCTCACTGGCCTGCCTTCGCGGTCAGCGCACGCACCGCGAAGTTGTTTCGTTTGCAGAGATTGGCTTAGACCTTCTGGGTGCCGCTTTACGGCCCCAGTACCACGAAAGCGAGCACGCGCAGGCTGGCCCAGt                                         >  1:2367961/1‑150 (MQ=255)
                                                                                                               aatgggTCTGGGGGCGCCTCCGGGGCCGGCCGGCGGGGCGGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTc                                        <  1:1919742/147‑1 (MQ=255)
                                                                                                                  ggggctggggTCGCCTCACGGGCCTGGGGGGGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                     <  1:3462516/143‑1 (MQ=255)
                                                                                                                  gggTCTGGGGGCGCCTCCGGGGCCGGCCGGGGGGGGCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                     <  1:2489953/150‑1 (MQ=255)
                                                                                                                    gTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                                   >  1:1017238/1‑150 (MQ=255)
                                                                                                                    gTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                                   >  1:1017259/1‑150 (MQ=255)
                                                                                                                        ggggTCGCCTCCCGGGCCTGCCTGCGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCGTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTc                               <  2:2101685/149‑1 (MQ=255)
                                                                                                                          ggTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                             <  1:1830361/150‑1 (MQ=255)
                                                                                                                          ggTCGCCTCACTGGCCTGCCGGCGCGGACTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                             <  2:4237929/150‑1 (MQ=255)
                                                                                                                            tCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCgcgc                           >  2:2752921/1‑150 (MQ=255)
                                                                                                                               acTCACTGGCCTGCCTTAGCGGTCTTCGCAGGCACTGCGAAGTTGAATCGTTTGCAGAGCTTGTCATCTACCTTCTTGGTGCCCCTTTACGCCCCCTGTACCTCGAACGAGAGCAGGCGCATGCTGGCCCAGTCCTCGTCTTCGCGCaac                        >  2:1759557/2‑148 (MQ=255)
                                                                                                                                cTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCACGAAAGCGAGCAGGCGCTGGCTGGCCCAGTCCACGTCTTAGCGCAtag                       >  1:2064454/1‑148 (MQ=255)
                                                                                                                                   agtGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGa                    <  2:4358627/148‑1 (MQ=255)
                                                                                                                                   aCTGGCCGGCCTGCGCGTTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGa                    <  2:1981439/150‑1 (MQ=255)
                                                                                                                                    cTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa                   <  2:1845961/150‑1 (MQ=255)
                                                                                                                                     tGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAg                  >  2:3543571/1‑150 (MQ=255)
                                                                                                                                      gccgtgcctgcGCAGACCGCGCAGGCAGGCGAAGTTGTTGCGTTTGCAGAGCGTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCAGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCACCGCTTCGCGCATATGGAAGAc                >  2:3318436/5‑150 (MQ=255)
                                                                                                                                         gtgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa              >  2:1022122/2‑150 (MQ=255)
                                                                                                                                          tgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGaa             >  2:1190253/1‑150 (MQ=255)
                                                                                                                                          tgcctgcGCAGACCGCGCAGGCAGGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGaa             >  2:1402435/1‑150 (MQ=255)
                                                                                                                                            cctgcGCAGACCGCGCAGGCAGGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGc           >  2:3024454/1‑150 (MQ=255)
                                                                                                                                            cctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAg                  >  2:3083619/1‑143 (MQ=255)
                                                                                                                                            cctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAg                  <  1:3083619/143‑1 (MQ=255)
                                                                                                                                             ctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCg          >  2:1358522/1‑150 (MQ=255)
                                                                                                                                                   ggTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTc    <  1:2781197/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGG  >  NC_002947/1464056‑1464354

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: