Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,469,161 A→C 12.1% P471P (CCT→CCG algX ← alginate biosynthesis protein AlgX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,469,1610AC12.1% 638.6 / 6.7 289P471P (CCT→CCGalgXalginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base A (150/74);  new base C (4/41);  total (154/135)
Fisher's exact test for biased strand distribution p-value = 1.69e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

CTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGGTGTCGGA  >  NC_002947/1469012‑1469310
                                                                                                                                                       |                                                                                                                                                     
cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                                       <  2:2586264/150‑1 (MQ=255)
cTTTTCGATCCCATGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                                       <  1:1965385/150‑1 (MQ=255)
 ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                                      >  1:5409288/1‑150 (MQ=255)
  tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                                     >  2:5345637/1‑150 (MQ=255)
  tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                                     >  1:293349/1‑150 (MQ=255)
  tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                                     >  2:808114/1‑150 (MQ=255)
   ttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGc                                                                                                                                                    >  1:3104560/1‑150 (MQ=255)
   ttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGGGCGCGGCACCACCAAACACGGCCAACGCGAGCAGGGCGGGGCGGGAGTTTCGTTTGAAAAGGGTCATGGGTGGCCTCTCTGGCCTGCCTGCGCGGTCTGCGCgggg                                                                                                                                                    >  2:4472568/1‑149 (MQ=255)
    tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCa                                                                                                                                                   <  2:2096117/150‑1 (MQ=255)
      gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCGTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                                  <  2:5334161/150‑2 (MQ=255)
      gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                                 >  2:342649/1‑150 (MQ=255)
       aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACAAAACAAGGCCAACGCGAGCAGGGCGGGGCTGTAATTTCTTTTGAAATGGGTCTTGGGTCGCCTCACTGGCGTGCCTGCGGGGTCTGCGCAGGCACgg                                                                                                                                                <  2:3048515/150‑1 (MQ=255)
         ccccTGTTTGTAACCGTGCGGCGGTCCCCGCGGGGGGTGCGCGGCACCACCAAACAAGGCCAAGGCGAGCTGGGGGGGGGTGTAATTTCTTTTGAAATGGGTCTTGGTTCGCCTCAGTGGCGTGCGGGCGCGGTCGGCGCCGGCACGGCg                                                                                                                                              <  2:3154171/150‑1 (MQ=255)
          cccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGa                                                                                                                                             >  2:3544011/1‑150 (MQ=255)
           ccTGGTAGTACCCCTGGGGCGGTCCCAGCGGGGCGGGCGCGGCACCACCAACCCCGCCCACCGCGCGCTGGGGGGGGCTGGAATTTCGTTTGAAACGGGTCATGGGTCGCCTCACTGGCCGGCGTGCGCGGTGTGCGCAGGCACGGCGaa                                                                                                                                            <  2:4910289/150‑1 (MQ=255)
           ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACACGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCGTGCCTGCGCAGACCGCGCAGGCaggccagt                                                                                                                                            <  2:3581302/150‑8 (MQ=255)
           ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                            >  1:3500413/1‑150 (MQ=255)
            cTGGTAGTAACCCTGCGGGGGTACCTGCGCGGCGTGCGCGGCACCACAAAACAAGGCCAACGCGAGCTGGGGGGGGCTGGAATTTGTTTTTAAATGGGTAAGGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAg                                                                                                                                           <  2:6229646/150‑1 (MQ=255)
            cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACAGGCCTGCCTGCGCGGTCGGCGCAGGCACGGCGAAg                                                                                                                                           >  2:659843/1‑150 (MQ=255)
             tGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgt                                                                                                                                          >  1:796160/1‑150 (MQ=255)
                   taACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGc                                                                                                                                    >  2:4706256/1‑150 (MQ=255)
                     aCCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGt                                                                                                                                  >  1:1858559/1‑150 (MQ=255)
                       ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                                >  1:1044967/1‑150 (MQ=255)
                        cggggggggTACCCGCGGGGGGGGCGGGGCACCACCAAACACGGCAACCGCGTGCTGGGGGGGGGGGTAATTTGTTTTGAAATGGGTAATGGGTCGCCTCCGGGCCGTGCCGGGGGGGTCTGGCAAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                               <  2:3116325/143‑1 (MQ=255)
                         ggcggggGTCCCCGGGGGGGGTGCGCGGCCCCCCCACCCCGGGCCACCGCGCGCTGGGGGGGGCTGGAATTTCGTTTGAAACGGGGCTGGGGGCGCCTCCGGGGCCGGCGGGCGCGGGCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGc                                                                                                                              <  2:6324296/149‑1 (MQ=255)
                                gTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCtt                                                                                                                       >  1:1110992/1‑150 (MQ=255)
                                 tACCAGCGCGGCGTGCGCGGCACCACCAAACACGGCCACCGCGTGCTGGGGGGGGCTGTAAATTCGTTTGAAATGGGTCTTGGGTCGCCTCACGGGCCTGCCTGCGCGGGCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTg                                                                                                                      <  2:1662587/150‑1 (MQ=255)
                                 tACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTg                                                                                                                      >  1:6736265/1‑150 (MQ=255)
                                 tACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTGTGAACAGGGTCCTGGGTCGCCTCACTGGCCTGCCTCCGCGGGCTGCGCAGGCACGCCGAAGGTGTTGCGGTTGCAGAGCGTg                                                                                                                      >  2:1793137/1‑150 (MQ=255)
                                    cAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCt                                                                                                                   >  1:5472942/1‑150 (MQ=255)
                                      gcgcGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                                 >  2:4220114/1‑150 (MQ=255)
                                       cggggcgTGCGCGGCACCACCAAACACCGCCAACGCGAGCTGGGCGGGGGTGGAAATTCGTTTGAAAAGGGTCTTGGGTCGCCTCCGTGGCCTGCGTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCg                                                                                                                <  2:6825127/147‑1 (MQ=255)
                                       cgcgGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCg                                                                                                                >  1:916507/1‑150 (MQ=255)
                                        ggggcgTGCGCGCCACCACCAAACACGGCCACCGCGAGCAGGGCGGGGCTGTAATTTCGTTTTAAATGGGTATTGGGTCGCATCAGTGGCCTGCCTGGGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGa                                                                                                               <  2:207613/148‑1 (MQ=255)
                                            cgTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCtt                                                                                                           >  1:5966363/1‑150 (MQ=255)
                                              tGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCt                                                                                                         >  1:4708637/1‑150 (MQ=255)
                                                ggcgcCACCACAAAACACGCCCACCGCGCGCTGGGGGGGGCTGCAATTTCGTTTCAAATGGGTCTGGGGGCGCCTCAGTGGCCGGCCGGCGGGGTGGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                                       <  1:4883396/149‑1 (MQ=255)
                                                   gccaccacAAAACACGGCCACCGCGACCTGGGCGGGGCTGCAAATTCTTTTCAAACGGGTCTGGGGCCGCCTCACGGCCCTGCCTGCGCGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTg                                                                                                    <  1:561976/148‑1 (MQ=255)
                                                     caccacCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGAGCCGGCACGGCGAATTTGTTGCGTTTGCAGAGCATGGCTTCGAACTTCAGGGTGcc                                                                                                  >  1:4023806/1‑150 (MQ=255)
                                                         accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                                              >  1:5789885/1‑150 (MQ=255)
                                                         accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                                              >  1:6289147/1‑150 (MQ=255)
                                                         accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGGTTGCAGAGCTTGGCTTCGACCGTCTGGGTGCCGCtt                                                                                              >  1:4278403/1‑150 (MQ=255)
                                                          ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCttt                                                                                             >  2:5420936/1‑150 (MQ=255)
                                                          ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTAGCAGAGCTTGGCTTCGACCATCTGGGTGCCCCttt                                                                                             >  1:5668719/1‑150 (MQ=255)
                                                           caAACACGGCCACCGCGAGCTGGGGGGGGCTGTAATTTCGTTTGAAAAGGGTCAGGGGTCGCCCCACGGGCCGGCCTGCGGGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTc                                                                                            <  2:3765052/150‑1 (MQ=255)
                                                           caAACACGCCCAACGCGAGCCGGGGGGGGGTTTAATTTGTTTTTAAACGGTTCTTGGGTCGCCTCACTGGCCGGCGTGCGGGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTc                                                                                            <  2:1922991/150‑1 (MQ=255)
                                                           caAAAAAGGCCAACGCGAGCCGGGGGGGGCTGTAAATTCGTTTGAAAAGGGTCTTGGGTCGCCTCAGTGGCCGGCCTGCGCGGTCTGCCCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTc                                                                                            <  2:5261739/150‑1 (MQ=255)
                                                            aaccacggCCACCGCGCGCAGGGCGGGGCTGTAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTcc                                                                                           <  2:4543349/144‑1 (MQ=255)
                                                             aaaccggCCAACCGCGTGCTGGGGGGGGTTGTAATTTGGTTTTAAATGGTACCTGGGCCGCCCCATGGGCCGGCGGGCGGGGTCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTcc                                                                                           <  1:6645600/145‑1 (MQ=255)
                                                                 cgcccACCGCGCGCTGGGGGGGTCTGTAATTTCTTTTCAAAGGGTTGCGGGGCCGCCTCCGGGGCGGGCCGGCGGGGGCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                                      <  1:1288475/147‑1 (MQ=255)
                                                                    ccAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTg                                                                                   >  2:2451020/1‑150 (MQ=255)
                                                                       tccgcgcgCTGGGCGGGGCTGGAAATTGGTTTGAAATGGGTCATGGGGCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTa                                                                                 <  2:6546802/148‑1 (MQ=255)
                                                                       aCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCAGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTAc                                                                                >  1:3685144/1‑150 (MQ=255)
                                                                         cgcgAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCATCTGGGTGCCGCTTTCCGGCCCCAGTAca                                                                               >  1:6546802/1‑149 (MQ=255)
                                                                                  gggggggCTGTAAATTCGTTTTAAATGGTTCATGGGTCGCCTCACTGGCCGGCCTGCGCGCTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAcgcg                                                                      <  2:4027648/147‑1 (MQ=255)
                                                                                  ggCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAcgcg                                                                      >  2:5965585/1‑150 (MQ=255)
                                                                                   ggggggTTGTAAATTCTTTTTAAAAGGTTCTTGGGTCGCCTCACGGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                                     <  2:6652338/148‑1 (MQ=255)
                                                                                   ggggggGTGTAATTTGGTTTGAATTGGGTCTGGGGTCGCCTCTGTGGCGTGCGTGCGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                                     <  1:558500/148‑1 (MQ=255)
                                                                                   ggggggGTGGAAATGCGTTTGAACGGGGTCTGGGGTCGCCGCCCGGGCGGGCCGGCGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                                     <  2:6073263/148‑1 (MQ=255)
                                                                                   gCGGGGCTGTAATTTCGTTTCAAATGGGTCAGGGGTCGCCTCACGGGCCTGCCTGCGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                                     <  1:2739175/150‑1 (MQ=255)
                                                                                    gggggGTGTAATTTCGTTTTAAACGGGTCAGGGGTCGCCCCAGGGGCCGGCCTGCGCGGTCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAg                                                                    <  1:3959281/149‑1 (MQ=255)
                                                                                            caaaTTGTTTTTAAACGGGTCTGGGGGCGCCCCCGGGGCGGGCCGGCGGGGGCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgct                                                            <  1:5163524/149‑1 (MQ=255)
                                                                                              atttCTTTTGAAATGGGTCATGGGTCGCCTCCGGGGCCGGCGGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctgg                                                          <  2:1388482/148‑1 (MQ=255)
                                                                                                 tggTTTGAAACGGGTCCGGGGTCGCCTCCCGGGCCGCCCTGCGCGGGCCGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctggctg                                                       <  2:4761632/148‑1 (MQ=255)
                                                                                                   gTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCaggctggc                                                     >  2:517514/1‑150 (MQ=255)
                                                                                                    tttGAAAGGGGCCTTGGGGCCCCTCCGGGGCGGGCCGGCGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                                    <  2:5834662/150‑1 (MQ=255)
                                                                                                    tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                                    >  1:6247284/1‑150 (MQ=255)
                                                                                                      tGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa                                                  >  2:2849202/1‑150 (MQ=255)
                                                                                                        aatcgggTCAGGGGTCGCCCCCCGGGCCGGCCTGGGGGGTCTGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGt                                                <  1:239488/146‑1 (MQ=255)
                                                                                                        aaaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCAGTACCTCGAACGCGAGCAGGCGCTGGCTGGccccgg                                                >  1:4363248/1‑147 (MQ=255)
                                                                                                            gggacttGGGTGGCCTCCGGGGCGTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                            <  2:5783223/144‑1 (MQ=255)
                                                                                                            gggTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                            <  2:2738978/150‑1 (MQ=255)
                                                                                                            gggTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                            <  2:1110992/150‑1 (MQ=255)
                                                                                                            gggTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTATGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCCGTCCTc                                            >  2:791997/1‑150 (MQ=255)
                                                                                                              gTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                                          >  2:1423955/1‑150 (MQ=255)
                                                                                                               tCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTc                                         >  2:2505456/1‑150 (MQ=255)
                                                                                                                 ttGGGTCGCCTCACGGGCCTGCCTGCGCGGCCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCtt                                       <  2:4629047/149‑1 (MQ=255)
                                                                                                                 aTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCtt                                       >  1:1640712/1‑150 (MQ=255)
                                                                                                                    ggTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                                       <  2:5681547/83‑1 (MQ=255)
                                                                                                                      tCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCACGCACGGCGAAGTTGTTGCGTTTTCAGAGCTTGGCTTCGACCTTATGGGTGCCGCTTTCCGGCCCCAGTAACTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTAgcgc                                  >  1:150101/1‑150 (MQ=255)
                                                                                                                      tCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACCGCGAAGTTGTTGCGTTTGCAGAGCTTGGCGTCGACCTGCTGGGTGCCGCGTTCCGGCCCCAGTAACTCGAACGCGAGCAGGCGCGGGCTGGCCCAGTCCTCGTCGTCgcgc                                  >  1:133090/1‑150 (MQ=255)
                                                                                                                             aCTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGa                           >  2:3263423/1‑150 (MQ=255)
                                                                                                                              ctggcctggcTGCGGGGTGTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa                          <  2:4832571/150‑1 (MQ=255)
                                                                                                                              cTGGCGGGCCTGCGGGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa                          <  1:2522389/150‑1 (MQ=255)
                                                                                                                              cTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa                          <  2:6452086/150‑1 (MQ=255)
                                                                                                                                ggCCCGCCTGGGGGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGa                        <  2:821123/150‑1 (MQ=255)
                                                                                                                                   cgtgcctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTAGACCATCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAACCGACCAGGCGCTGGCTGGCCCAGTCCTCGTCTTAGCGCCTATTGAAGACg                      >  2:296795/3‑150 (MQ=255)
                                                                                                                                   ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa                     >  2:5559468/1‑150 (MQ=255)
                                                                                                                                     tgcctgcGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGaa                    <  1:2107196/150‑1 (MQ=255)
                                                                                                                                     tgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGac                    >  1:6094730/1‑149 (MQ=255)
                                                                                                                                     tgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGaa                    >  1:2203736/1‑150 (MQ=255)
                                                                                                                                     tgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGaa                    >  1:2705967/1‑150 (MQ=255)
                                                                                                                                        ctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCg                 >  1:5296002/1‑150 (MQ=255)
                                                                                                                                         tgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGctggctgg                                                      <  1:2214693/112‑1 (MQ=255)
                                                                                                                                         tgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGctggctgg                                                      >  2:2214693/1‑112 (MQ=255)
                                                                                                                                          ccgcgtccTGCCCAGCCCCGGCGAAGCTGTTGCGTTGGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGc               <  1:4472568/143‑1 (MQ=255)
                                                                                                                                                tCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCgg         >  1:4964652/1‑150 (MQ=255)
                                                                                                                                                 cTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGgt        <  2:6131673/150‑1 (MQ=255)
                                                                                                                                                  tGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGgtg       >  1:4583399/1‑150 (MQ=255)
                                                                                                                                                       aGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGGTGTCGGa  <  2:3631834/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGGTGTCGGA  >  NC_002947/1469012‑1469310

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: