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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,464,1980TG13.2% 317.6 / 5.3 137Q476P (CAG→CCG) algXalginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (64/49);  new base G (2/20);  total (68/69)
Fisher's exact test for biased strand distribution p-value = 3.43e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

TCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGC  >  NC_002947/1464049‑1464345
                                                                                                                                                     |                                                                                                                                                   
tCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGgcctgcct                                                                                                                                                     >  1:2751647/1‑150 (MQ=255)
tCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGgactgcct                                                                                                                                                     >  1:311850/1‑150 (MQ=255)
   ccGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTgcg                                                                                                                                                  >  1:4558608/1‑150 (MQ=255)
    cGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTgcgc                                                                                                                                                 <  2:2967634/150‑1 (MQ=255)
    cGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCGgcgc                                                                                                                                                 >  1:1878832/1‑150 (MQ=255)
       tCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGt                                                                                                                                              <  2:2751647/150‑1 (MQ=255)
         ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                            >  2:24822/1‑150 (MQ=255)
         ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                            >  1:3711161/1‑150 (MQ=255)
         ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                            >  2:1738301/1‑150 (MQ=255)
         ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                            >  1:2855604/1‑150 (MQ=255)
          tGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTg                                                                                                                                           <  1:2271038/150‑1 (MQ=255)
           gAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGACATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                          >  1:351903/1‑150 (MQ=255)
           gAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                          >  2:3517921/1‑150 (MQ=255)
             gCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcgc                                                                                                                                        <  2:1249480/150‑1 (MQ=255)
              cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                       <  2:3624383/150‑1 (MQ=255)
              cTTTTCGAGCCCCTGGTAGTAACCCGGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCACCCACGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCCTGGGTCGCCTCCCTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                       <  1:192315/150‑1 (MQ=255)
               ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                      >  2:3165577/1‑150 (MQ=255)
               ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                      >  1:4847863/1‑150 (MQ=255)
                tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                     <  2:2537502/150‑1 (MQ=255)
                    gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                 <  2:4998912/150‑1 (MQ=255)
                    gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                 >  2:3453805/1‑150 (MQ=255)
                     aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCaggc                                                                                                                                <  1:3130305/150‑4 (MQ=255)
                     aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCaggc                                                                                                                                <  2:3559290/150‑4 (MQ=255)
                      tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGc                                                                                                                               >  2:167641/1‑150 (MQ=255)
                       ccccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCg                                                                                                                              >  1:2661261/1‑150 (MQ=255)
                        cccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGa                                                                                                                             >  1:4588572/1‑150 (MQ=255)
                         ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAAGTCGTTTGAAAAGGGTCAGGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                            >  2:1797697/1‑150 (MQ=255)
                          cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAg                                                                                                                           >  1:5212477/1‑150 (MQ=255)
                               attaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGCCCACCGCGAGCTGGGGGGGGCTGGAAATTCGTTTGAAATGGGTCTTGGGTCGCCCCACGGGCCGGCCGGCGCGGTCGGCGCAGGCACGGCGAAgttgtt                                                                                                                      <  2:2227239/148‑1 (MQ=255)
                               ataaaCCCTGCGGCGTTACCAGCGCGGCGTGCGCGCCACAACAAAAAAAGCCAAACCCTAGCAGGGCGGGGATGTAAATTCTTTTCAAAAGGTTAATGGGTCCCCTCAATGGCCTGCCTGCGCGGTCGGCGCAGGCACGGCGAAgttgtt                                                                                                                      <  2:2647759/147‑1 (MQ=255)
                               agtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgtt                                                                                                                      >  1:1201136/1‑150 (MQ=255)
                               agtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgtt                                                                                                                      >  2:5346800/1‑150 (MQ=255)
                                gtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCGTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGtggttg                                                                                                                     >  2:2755147/1‑150 (MQ=255)
                                gtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGt                                                                                                                                              >  1:2671792/1‑125 (MQ=255)
                                gtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGt                                                                                                                                              <  2:2671792/125‑1 (MQ=255)
                                gtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGt                                                                                                                                              <  2:2671786/125‑1 (MQ=255)
                                 taACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGc                                                                                                                    >  2:4107734/1‑150 (MQ=255)
                                  aaCCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCg                                                                                                                   >  1:5378201/1‑150 (MQ=255)
                                     ccggcggggGTACCCGCGGGGGGTGCGCGGCACCACCAAACACGGCCAACGCGAGCCGGGCGGGGCTGGAAATTGTTTTGAAATGGGTCTTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                <  2:330121/147‑1 (MQ=255)
                                     ccTGCGGCGTTACCTGCGCGGCGTGCGCGGCACCACCAAACAAGCCCACCGCTAGCCGGGCGGGGCTGTAATTGCGTTTCAAACGGGTCAGGGGTCGCCACACGGGCCGGCCTGCGGGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                <  2:5390478/150‑1 (MQ=255)
                                      cTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGGAGttgttgagatcg                                                                                                               >  1:498091/1‑144 (MQ=255)
                                                 ccAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGGTTGCAGAGCTTGGc                                                                                                    >  2:3649564/1‑150 (MQ=255)
                                                  cAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGGCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCt                                                                                                   >  1:536946/1‑150 (MQ=255)
                                                    gcgcGGCGTGCGCGGCACCACCACACACGGCCACCGCGAGCAGGGCGGGGCTGGAAATTGGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCGGCGGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                 <  2:4909490/150‑1 (MQ=255)
                                                    gcgcGGCGTGCGCGGCACCACCAACCACGGCCACCGGGCGCCGGGCGGGGCTGTAAATTCTTTTGAAATGGGTCTTGGGGCGCCCCCGTGGCCGGCCTGCGGGGCCGGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                 <  2:1624594/150‑1 (MQ=255)
                                                    gcgcGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                 >  2:4466427/1‑150 (MQ=255)
                                                    gcgcGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                 >  2:4466447/1‑150 (MQ=255)
                                                    gcgcGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGGTTGCAGAGCTTGGCTTc                                                                                                 >  2:1090090/1‑150 (MQ=255)
                                                      gcggcgTGCGCGCCACCCCCAAACCAGCCCACCGCGCGCGGGGCGGGGCTGGAATTTCTTTTTAAACGGGTCTTGGGTCGCCTCAGGGGCCTGCCTGGGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGa                                                                                               <  2:3711161/150‑1 (MQ=255)
                                                             gggcgcCACCACAAACCCCGGCCACCGCGTGCGGGGGGGGGCTTGATTTTGGTTTCAAATGGGTCCGGGGTCGCCTCCGGGGCCGGCCTGCGGGGCCTGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                        <  1:5172176/148‑1 (MQ=255)
                                                              cgcgGCACCACAAAACACGGCCACCGCGCGCCGGGGGGGGCTGTAATTTGTTTTGAAATGGTACTGGGGGCGCCCCACGGGCCTGCGGGCGGGGTCGGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                       <  1:4881273/150‑1 (MQ=255)
                                                                    accacaAAACACGGCCACCGCGAGCAGGGGGGGGCTGTAAATTGGTTTGAAACGGGTCTTGGGTCGCCTCACTGGCCGGCCGGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCg                                                                                 <  2:4756155/150‑1 (MQ=255)
                                                                      caccaAACCCGGCCAACGCGCGCGGGGGGGGGGTGGAAATTCTTTTGAAACGGGTCATGGGTCGCCTCACGGGCCTGCCTGCGGGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCt                                                                               <  2:2122507/150‑1 (MQ=255)
                                                                      caccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCag                                                                                                                                  <  1:908950/99‑2 (MQ=255)
                                                                       accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                              >  2:742649/1‑150 (MQ=255)
                                                                       acaaaaCAAGGCCAACGCGAGCAGGGCGGGGCTGTAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                              <  2:4910094/147‑1 (MQ=255)
                                                                        ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGACGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCttt                                                                             >  2:665273/1‑150 (MQ=255)
                                                                        ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCttt                                                                             >  2:1460086/1‑150 (MQ=255)
                                                                         caAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCAGCCCGGTCTGCGCAGGCACGGCGAAGATGTTGCGTTTGCAGAGCATGGCATCGAACTTCTGGGTGCCGCTATc                                                                            >  2:3006660/1‑150 (MQ=255)
                                                                          aaaCAAGCCCAACGCTAGCAGGGCGGGGCTGTAATTTCTTTTGAAAAGGGTCTTGGGTCGCCTCAGTGGCCTGCCTGCGCGGACTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTcc                                                                           <  2:799411/150‑1 (MQ=255)
                                                                           aaCAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCg                                                                          >  2:1060886/1‑150 (MQ=255)
                                                                           aaCAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCATCTGGGTGCCGCGTTCCg                                                                          >  1:3203725/1‑150 (MQ=255)
                                                                            aCAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCgg                                                                         >  1:3003131/1‑150 (MQ=255)
                                                                               cgcccACCGCCAGCAGGGCGGGGCGGGTAATTGGTTTGAAATGGGTCATGGGTCGCCTCCCGGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                      <  2:3061878/147‑1 (MQ=255)
                                                                               aGGCCAACGCGAGCAGGGCGGGGCTGTAAATTCGTTTTAAATGGGTCATGGGTCGCCTCACGGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                      <  2:4350331/150‑1 (MQ=255)
                                                                                 gCCACCGCGAGCTGGGCGGGGCTGTAAATTCTTTTTAAAAGGGTCATGGGTCGCCTCCCTGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCt                                                                    <  2:393133/150‑1 (MQ=255)
                                                                                  ccAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                            <  2:1544959/77‑1 (MQ=255)
                                                                                  ccAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                            >  1:1544959/1‑77 (MQ=255)
                                                                                    accgcgAGCTGGGGGGGGCTGTAATTTCTTTTCAAATGGGTCTTGGGCCGCCTCAGGGGCCTGCCGGGGGGGTGTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTa                                                                 <  1:4416049/148‑1 (MQ=255)
                                                                                    aaCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTAGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTa                                                                 >  1:4329282/1‑150 (MQ=255)
                                                                                     ccgcgcgCCGGGGGGGGCTGCAAATTCTTTTCAAATGGGTCTTGGGTCGCCACACTGCCCGGCCGGCGCGGTCGCCCCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTAc                                                                <  1:3895449/149‑1 (MQ=255)
                                                                                      cgcgAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTAcc                                                               <  2:4975789/150‑1 (MQ=255)
                                                                                       gcgcgcGGGGGGGGGCTGGAAATTGTTTTTAAAGGGGTCAGGGGTCGCCGCCGGGGCGTGCCGGGGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCt                                                              <  1:5263609/150‑1 (MQ=255)
                                                                                       gcgAGCTGGGGGGGGCTGTAATTTGTTTTGAAATGGGTCATGGGTCGCCTCCGTGGCCTGCCGGCGGGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCt                                                              <  2:1201136/150‑1 (MQ=255)
                                                                                        cgcgccgggCGGGGCTGGAAATTCGTTTGAAACGGGTCTTGGGTCGCCTCACTGGCCGGCGTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                             <  2:4466557/144‑1 (MQ=255)
                                                                                            ctgggCGGGGCTGGAAATTCTTTTGAAATGGGTCTTGGGTCGCCTCAGTGGCCTGCCTGCGGGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAc                                                         <  2:4243477/148‑1 (MQ=255)
                                                                                               ggCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAcgcg                                                      >  2:832322/1‑150 (MQ=255)
                                                                                                    gggttcaaaTTGTTTTGAATCGGGTCTTGGGTCGCCACCGCGGCCGCCCGGCGGGGCCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAg                                                 <  1:2856072/144‑1 (MQ=255)
                                                                                                       tGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCAGTACCTCGAACGCGAGCAGgcg                                              >  2:1587452/1‑150 (MQ=255)
                                                                                                         taaaTTCGTTTGACATGGGTCAGGGGGCGCCTCACTGGCCGGCGGGCGCGGCCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgct                                            <  2:1979962/149‑1 (MQ=255)
                                                                                                                    gAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                     <  1:2170593/50‑1 (MQ=255)
                                                                                                                       aaGGGTCTTGGGTCGCCTCCCGGGCCGGCGTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTcc                              <  2:5239835/150‑1 (MQ=255)
                                                                                                                        aGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCt                             >  2:5246193/1‑150 (MQ=255)
                                                                                                                         gggTCATGGGTCGCCTCACTGGCCTGCCTGTGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                            <  1:374696/150‑1 (MQ=255)
                                                                                                                          gttctggggTCGCCTCACGGGCCGGCCGGCGCGTTCTGCGCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCg                           <  1:4583864/144‑1 (MQ=255)
                                                                                                                           gTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                          >  1:1285935/1‑150 (MQ=255)
                                                                                                                              aTGGGTCGCCTCACTGGCCTGCCAGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCATCTGGGTGCCGCGTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCGCGTCtt                       >  1:3004469/1‑150 (MQ=255)
                                                                                                                                 ggTCGCCTCTGTGGCCTGCCTGGGGGTTGTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                    <  2:4786419/150‑1 (MQ=255)
                                                                                                                                       ctcaggggCCTGCCGGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCTTCGTCTTCGCGCATCt              <  2:351903/144‑1 (MQ=255)
                                                                                                                                        tCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTg             >  2:2663158/1‑150 (MQ=255)
                                                                                                                                         aaggggccggccTGCGGGGTGTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTgg            <  2:311850/146‑1 (MQ=255)
                                                                                                                                           cTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa          <  1:1738301/150‑1 (MQ=255)
                                                                                                                                              gccggcctGCGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGAc       <  2:1866233/150‑1 (MQ=255)
                                                                                                                                                ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa     >  1:3463855/1‑150 (MQ=255)
                                                                                                                                                 ggcgtgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGaa    <  2:2171581/146‑1 (MQ=255)
                                                                                                                                                  gccggcGGGGTCGGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAg   <  1:167641/150‑1 (MQ=255)
                                                                                                                                                  gccggcGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCCGGAAGACGAAg   <  2:1285935/150‑1 (MQ=255)
                                                                                                                                                   cctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGc  >  1:3686345/1‑150 (MQ=255)
                                                                                                                                                     tgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                          >  2:2167905/1‑124 (MQ=255)
                                                                                                                                                     tgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                          <  1:2167905/124‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCACCGGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGC  >  NC_002947/1464049‑1464345

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: