Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,469,152 G→A 5.5% T474T (ACC→ACT) ‡ algX ← alginate biosynthesis protein AlgX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,469,1520GA 5.5% 671.5 / 18.6 219T474T (ACC→ACT) ‡algXalginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base G (111/93);  new base A (2/10);  total (113/106)
Fisher's exact test for biased strand distribution p-value = 1.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 231 total aligned reads displayed.

GGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGT  >  NC_002947/1469003‑1469300
                                                                                                                                                     |                                                                                                                                                    
ggTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGcgcg                                                                                                                                                      >  2:2217836/1‑150 (MQ=255)
ggTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGcgcg                                                                                                                                                      <  2:3567687/150‑1 (MQ=255)
 gTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGACATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCgg                                                                                                                                                     >  2:3542646/1‑150 (MQ=255)
 gTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCgg                                                                                                                                                     >  2:4273499/1‑150 (MQ=255)
    ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                                  >  1:4111699/1‑150 (MQ=255)
    ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                                  >  1:3971245/1‑150 (MQ=255)
      gAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                                <  2:4992457/150‑1 (MQ=255)
        gCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcgc                                                                                                                                              <  2:796840/150‑1 (MQ=255)
         cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAACCAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                             <  1:1536189/150‑1 (MQ=255)
          ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                            <  1:990879/150‑1 (MQ=255)
             tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCa                                                                                                                                         <  1:1476497/150‑1 (MQ=255)
             tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCa                                                                                                                                         <  1:4477452/150‑1 (MQ=255)
             tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCa                                                                                                                                         <  2:5518957/150‑1 (MQ=255)
              cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCag                                                                                                                                        <  2:3748891/150‑2 (MQ=255)
              cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                        <  1:2945028/150‑2 (MQ=255)
              cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                        <  1:3361975/150‑2 (MQ=255)
              cGATCCCCTGGTAGGAACCCTGCGGCGGGACCAGCGCGGCGTGCGCGGCACCACCACACAAGGCCACCGCGAGCAGGGCGGGGCTGGAACTTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                        <  2:3402219/150‑2 (MQ=255)
               gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                       >  2:3691894/1‑150 (MQ=255)
                aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                      >  1:3545992/1‑150 (MQ=255)
                aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                      >  1:4993592/1‑150 (MQ=255)
                 tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGc                                                                                                                                     >  2:360739/1‑150 (MQ=255)
                  ccccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCGTGCCTGCGCAGACCGCGCAGGCaggcca                                                                                                                                    <  1:3091624/150‑6 (MQ=255)
                   cccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCgg                                                                                                                                   >  1:2598016/1‑149 (MQ=255)
                    ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTTGAAATTCTTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                  <  2:404583/150‑1 (MQ=255)
                    ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAACTTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                  >  1:2215712/1‑150 (MQ=255)
                    ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                  >  2:4591955/1‑150 (MQ=255)
                     gtGGTTGTAACCCTGGGGGGGTACCAGCGCGGGGTGCGCGGCACCACAAAAAATGGCAAACGGTTGCTGGGGGGGGCTGTAAATTGTTTTTAAATGGGTCTTGGGTCGCCTCACTGGCCTGCCTGCGGGGTCTGCGCAGGCACGGCGAAg                                                                                                                                 <  2:5050829/149‑1 (MQ=255)
                     cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAg                                                                                                                                 >  1:4752656/1‑150 (MQ=255)
                     cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGAGCGGGCATCGCAGGCACCGCGAAg                                                                                                                                 >  1:1072621/1‑150 (MQ=255)
                     cTGGGAGTAACCCTGCGGCGGTACCAGAGCGGCGTGCGCGGCACAACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTAGTTTGAAAAGTGTCATGGGTCGCCTCACTGGCCTGCCAGCGCGGTCAGCGCCGGCACGGCGCAg                                                                                                                                 >  1:1184229/1‑150 (MQ=255)
                        gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAggtg                                                                                                                              >  1:2096372/1‑150 (MQ=255)
                           gtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGACGttgttg                                                                                                                           >  2:1930004/1‑150 (MQ=255)
                             aaCCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCg                                                                                                                         >  2:4742371/1‑150 (MQ=255)
                              aCCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGACGTTGTTGCGt                                                                                                                        >  2:167745/1‑150 (MQ=255)
                              aCCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAACTCGTTTGAAAAGGGTCATTGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGt                                                                                                                        >  2:3969577/1‑150 (MQ=255)
                                ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                      >  1:837916/1‑150 (MQ=255)
                                ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                      >  1:588456/1‑150 (MQ=255)
                                ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                      >  1:2012640/1‑150 (MQ=255)
                                 cTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                     >  1:5392094/1‑150 (MQ=255)
                                      gcggTACCCGCGCGGCGTGCGCGCCACCACAAAACACGCCCACCGCGAGCTGGGGGGGGCTTCAATTTCTTTTCAAATGGGTCTTGGGCCGCCAAACGGGCCGGCCTGCGCGGTGTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCagag                                                                                                                <  2:2565384/150‑1 (MQ=255)
                                        ggTACCTGCGGGGCTTGCGCGGCACCACCAAACACGCCCAACGCGTGCTGGGGGGGGGTGTAATTTCTTTTGAAATGGGTCTTGGTTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCt                                                                                                              <  2:4899661/150‑1 (MQ=255)
                                        ggTACCAGCGCGGCGTGCGCGCCACCACCAAACACGGCCAACGCGAGCAGGGCGGGGCTGTAAATTCTTTTGAAACGGGTCATGGGTCGCCTCACGGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCt                                                                                                              <  2:4529995/150‑1 (MQ=255)
                                         ttACCCGCGGGGCGTGGGCGGCACCACCAAACTCGCCAAAGGCTCGCAGGGGGGGGCTGTAATTTCGTTTTAAATGGTTCTTGGGTCGCCCCCGGGGCGGGCGGGCGGGGTTGGGGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCtt                                                                                                             <  2:1004241/149‑1 (MQ=255)
                                         gTACCCGGGCGGCGTGCGGGCCCCCCCCAACCAAGGCCCCGGCGTGCGGGGGGGGGCTGTAATTTGTTTTGAACCGGGTATGGGGGGGCCCCCGGGGCGGGCGGGCGCGTCCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCtt                                                                                                             <  2:1004075/150‑1 (MQ=255)
                                         gTACCCGGGCGGCGGGCGCGCCCCCACCACCCACGGCCACCGCGCGCGGGGGGGGGCTGTAAATTCGTTTCAAACGGTGCTGGGGTCGCCCCAGGGGCCGGCCGGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCtt                                                                                                             <  2:2622910/150‑1 (MQ=255)
                                          tACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGcgcg                                                                                                                                                      >  2:1111723/1‑108 (MQ=255)
                                              aGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                       >  1:4640484/1‑150 (MQ=255)
                                               gcgcGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCACCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                       >  2:1311367/1‑150 (MQ=255)
                                                 gcggcgTGCGCGGCACCACAAACCAAGGCCACCGCGAGCCGGGGGGGGTGGGAAATTCTTTTGAAACGGGTCATGGGTCGCCTCACGGGCCGGCGTGCGGGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGa                                                                                                     <  2:2882855/150‑1 (MQ=255)
                                                        gcgcgcCACCCCCAACCCCGCCCACCGCGCGCGGGGGGGGGGTGTAATTTCTTTTCAAATGGGCCTGGGGCCGCCGCCGGGGCCGGCGGGGGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                              <  1:1047809/150‑1 (MQ=255)
                                                        gcgcgcCACCACCAAACAAGGCAAACGCGAGCCGGGGGGGGCTGTAAATTCTTTTTAAAAGGGTCATGGGTCGCCTCACTGGCCGGCCTGCGCGTTTTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                              <  2:5377174/150‑1 (MQ=255)
                                                        gcgcGGCACCACCAAACACGGCCAACGCGAGCAGGGCGGGGCTGGAATTTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                              <  2:5108812/150‑1 (MQ=255)
                                                           cgccaccacCAAACACGGCCACCGCGAGCAGGGGGGGGCTGGAAATTCGTTGGAACAGGGTCTTGGGTCGCCCCACGGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGt                                                                                           <  2:1294216/147‑1 (MQ=255)
                                                             gCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGGTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGc                                                                                         >  1:3830643/1‑150 (MQ=255)
                                                              caccacCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGcc                                                                                        >  2:4058385/1‑150 (MQ=255)
                                                                    caAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCAGGGTGCCGCTTTc                                                                                  >  1:2930833/1‑150 (MQ=255)
                                                                     aaacccggCCAACGCGAGCCGGGGGGGGCTGTAATTTGTTTTTAAATGGGTCAGGGGTCGCCTCACTGGCCTGCGGGGGGGGTCGGCGCTGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTcc                                                                                 <  1:1813070/144‑1 (MQ=255)
                                                                     aaaCAAGGCCAACGCGAGCTGGGGGGGGCTGTAAATTGGTTTTAAATGGGTCTTGGGCCGCCTCACTGGCCTGCGTGGGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTcc                                                                                 <  1:2746613/150‑1 (MQ=255)
                                                                     aaaCAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTcc                                                                                 >  2:4833503/1‑150 (MQ=255)
                                                                         acgcccACCGCGCGCGGGGGGGGGCTGGAACTTCGTTTTAAACGGTCCATGGGGCGCCTCAGTGGCCGGCCGGCGGGGGGGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGcc                                                                             <  2:2011954/146‑1 (MQ=255)
                                                                            gCCAACGCGAGCAGGGCGGGGCGGGAAATTCGTTTGAAATGGGTCATGGGGCGCCTCACGGGCCTGCGTGCGCGGTCTGGGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCt                                                                          <  2:3197601/150‑1 (MQ=255)
                                                                              cAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGt                                                                        >  1:4827303/1‑150 (MQ=255)
                                                                               aaCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCAGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCAGGGTGCCGCTTTCCGGCCCCTGTa                                                                       >  2:1732322/1‑150 (MQ=255)
                                                                                   ggcgctggggggggCTGGAATTTCTTTTTAAATGGGTACTGGGTGGCCAAACTGGCGTGCGGGCGCGGCCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                                   <  2:60688/144‑1 (MQ=255)
                                                                                    gAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCg                                                                  >  2:4330028/1‑150 (MQ=255)
                                                                                        aGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCATCTGGGTGCCGCTTTACGGCCCCTGTACCTCGAAcg                                                              >  1:2464659/1‑150 (MQ=255)
                                                                                           ggggggTTGTAATTTCGTTTAAAACGGGACCTGGGGCCCCTCACGGGCCGCCCTGCGGGTCCGGCGCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                           <  1:2574791/148‑1 (MQ=255)
                                                                                              gggCTGGAAATTCGTTTGAAATGGGTCTTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCa                                                        <  2:3702971/150‑1 (MQ=255)
                                                                                               ggCTGTAAATTGTTTTTAAATGGGTCTTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAg                                                       <  1:4675386/150‑1 (MQ=255)
                                                                                                gCTGTAATTTCTTTTTAAATGGGCCTGGGGCCGCCCCCCGGGCCTGCGGGCGGGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAgg                                                      <  1:224457/150‑1 (MQ=255)
                                                                                                   ggAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAAGCGAGCAGgcgc                                                   >  1:3630850/1‑150 (MQ=255)
                                                                                                    gAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCACGCACCGCGAAGTTGTTTCGGTTTCAGAGCATGGCGTCGAACTTCTTGGTGCCGCTTTCCGGCCCCAGTAACACGAAACCGAGCAGgcgca                                                  >  1:195689/1‑149 (MQ=255)
                                                                                                          cttttCAAACGGGTCAGGGGTCGCCCCACTGGCCGGCCTGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGctggctgg                                            <  2:1401637/148‑1 (MQ=255)
                                                                                                           ttttGAAATGGGTCTGGGGTCGCCGCACGGGCCTGCGGGGGCGGTGGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCtggctggc                                           <  1:3914299/149‑1 (MQ=255)
                                                                                                            tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGTTGTTGTTGGGTTTGGTGAGGTTGGCTTCGAGGTTGTGGGTGGCGGTTTGCGGCCCCTGTTCCTCGTAGGCGAGGAGGGGGTGGGTGGcc                                          >  2:119831/1‑150 (MQ=255)
                                                                                                            tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                          >  1:3258714/1‑150 (MQ=255)
                                                                                                               gAAATGGGTCCTGGGCCGCCCCCCGGGCCGGCGTGCGGGGTCTGCGCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAg                                       <  2:1404874/150‑1 (MQ=255)
                                                                                                               gAAACGGGTCTTGGGTCGCCTCACTGGCCGGCCTGCGCGGGCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAg                                       <  2:2948185/150‑1 (MQ=255)
                                                                                                               gAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAg                                       >  2:124161/1‑150 (MQ=255)
                                                                                                                aaaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGt                                      >  1:4528101/1‑150 (MQ=255)
                                                                                                                   tgggTCTTGGGTCGCCTCAGTGGCCGGCGGGCGGGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCt                                   <  1:578704/149‑1 (MQ=255)
                                                                                                                      gTCCTGGGTCGCCTCACTGGCCTGCCTGGGCGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                                <  2:3830643/150‑1 (MQ=255)
                                                                                                                         aTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCAGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCtt                             >  1:1373072/1‑150 (MQ=255)
                                                                                                                          tGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTACCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTc                            <  2:4262952/150‑1 (MQ=255)
                                                                                                                           gggTCGCCTCAGGGGCGTGCGTGCGGGGTCGGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcg                           <  1:4273499/150‑1 (MQ=255)
                                                                                                                           gggTCCCCACACGGGCCTGCCTGCGCGGACTGCCCAGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcg                           <  1:4163860/150‑1 (MQ=255)
                                                                                                                            ggTCGCCTCACTGGCCTGCCAGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                          >  2:1659475/1‑150 (MQ=255)
                                                                                                                                gCCTCAGGGGCCTGCGTGGGCGGTCTGGGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCAt                      <  2:3771158/150‑1 (MQ=255)
                                                                                                                                    cACTGCCCCCCCTCCGCGGTCTCCCAAGCCACGGCGAAGTTGTTGCGATTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTgg                  <  2:588456/150‑1 (MQ=255)
                                                                                                                                      cgggCCTGCCGGCGCTGGATGCGAAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACTCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa                <  2:4452857/148‑1 (MQ=255)
                                                                                                                                      cTGGCCTGCCTGCGCGGTCTGTGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa                <  1:5109724/150‑1 (MQ=255)
                                                                                                                                      cTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa                <  2:4113997/150‑1 (MQ=255)
                                                                                                                                         gcctgcctGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGAc             >  1:2365505/1‑150 (MQ=255)
                                                                                                                                           tttcccgcgcGATCAGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa           >  2:5226261/7‑150 (MQ=255)
                                                                                                                                           ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa           >  1:5109989/1‑150 (MQ=255)
                                                                                                                                            tgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGATGACGaa          <  2:1513458/150‑1 (MQ=255)
                                                                                                                                            tgcctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGaa          >  2:1144195/1‑150 (MQ=255)
                                                                                                                                                 gcgcGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGc     >  2:1848194/1‑150 (MQ=255)
                                                                                                                                                 gcgcGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGc     >  2:5429929/1‑150 (MQ=255)
                                                                                                                                                    cgGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGt  >  2:2460911/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGT  >  NC_002947/1469003‑1469300

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: