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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,469,1540TA 5.2% 834.6 / 27.3 307T474S (ACC→TCC) algXalginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (157/119);  new base A (3/12);  total (166/141)
Fisher's exact test for biased strand distribution p-value = 6.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-03

Only 100 of 310 total aligned reads displayed.

TCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGG  >  NC_002947/1469005‑1469303
                                                                                                                                                      |                                                                                                                                                     
  ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                                     >  1:2762139/1‑150 (MQ=255)
  ttGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCt                                                                                                                                                     >  1:6259824/1‑150 (MQ=255)
     aGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcg                                                                                                                                                  >  2:2539048/1‑150 (MQ=255)
     aGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcg                                                                                                                                                  <  1:3500969/150‑1 (MQ=255)
      gCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcgc                                                                                                                                                 <  2:3867368/150‑1 (MQ=255)
      gCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcgc                                                                                                                                                 <  2:4863701/150‑1 (MQ=255)
       cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTc                                                                                                                                                      <  2:1923420/144‑1 (MQ=255)
       cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                                <  2:2586264/150‑1 (MQ=255)
       cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                                >  2:5321762/1‑150 (MQ=255)
       cTTTTCGATCCCATGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCa                                                                                                                                                <  1:1965385/150‑1 (MQ=255)
        ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                               >  1:688652/1‑150 (MQ=255)
         tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                              >  1:1225266/1‑150 (MQ=255)
         tttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                              >  1:293349/1‑150 (MQ=255)
          ttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGc                                                                                                                                             >  1:3104560/1‑150 (MQ=255)
          ttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCGGCGCAGGc                                                                                                                                             >  1:4716151/1‑150 (MQ=255)
            cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCAc                                                                                                                                           <  1:3411563/150‑1 (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)
                 cccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGa                                                                                                                                      >  2:433901/1‑150 (MQ=255)
                  ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACACGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCGTGCCTGCGCAGACCGCGCAGGCaggccagt                                                                                                                                     <  2:3581302/150‑8 (MQ=255)
                  ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                     >  2:2307268/1‑150 (MQ=255)
                  ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                     >  1:5415203/1‑150 (MQ=255)
                   cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAg                                                                                                                                    >  2:659871/1‑150 (MQ=255)
                    tGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgt                                                                                                                                   >  2:6523901/1‑150 (MQ=255)
                       aagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgt                                                                                                                                >  2:2463064/2‑150 (MQ=255)
                        agtaCCCCTGCGGCGTTACCAGCGCGGCGTGCGCGGCACCACCAAACACGGCCAACGCGAGCTGGGCGGGGCTGGAAATTCGTTTGAAAAGGGCCATGGGTCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgtt                                                                                                                               <  2:1097643/150‑1 (MQ=255)
                              ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                         >  1:1044967/1‑150 (MQ=255)
                               cggggggggTACCCGCGGGGGGGGCGGGGCACCACCAAACACGGCAACCGCGTGCTGGGGGGGGGGGTAATTTGTTTTGAAATGGGTAATGGGTCGCCTCCGGGCCGTGCCGGGGGGGTCTGGCAAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                        <  2:3116325/143‑1 (MQ=255)
                               cgggggcggTACCCGGGCGGCGGGCGCGCCACCACCAAACCAGCCAAACGCGTGCCGGGGGGGGCTGAAATTTGGTTTTAAATGGGTCTTGGGGCGCCTCACGGGCCGGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                        <  2:5833130/146‑1 (MQ=255)
                               cggcggcgGAACCCGCGGGGCGCGCGCGGCCCCACCAAACCTGGCCCCCGCGAGGCGGGCGGGGCTGTAATTTGGTTCGAAATGGTTCTTGGGTCGCCCCAGGGGCGGGCCGGGGCGGTGTGCACAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                        <  2:4825894/148‑1 (MQ=255)
                               cTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTTGGtttt                                                                                                                        >  1:6338286/1‑149 (MQ=255)
                                        tACCAGCGCGGCGTGCGCGGCACCACCAAACACGGCCACCGCGTGCTGGGGGGGGCTGTAAATTCGTTTGAAATGGGTCTTGGGTCGCCTCACGGGCCTGCCTGCGCGGGCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTg                                                                                                               <  2:1662587/150‑1 (MQ=255)
                                            aGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGGCGGCGCAGGCACCGCGACGTGGTTGCGTTTGCAGGGCTGGGcgt                                                                                                           >  1:189826/1‑148 (MQ=255)
                                              cgcgGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCg                                                                                                         >  1:3999310/1‑150 (MQ=255)
                                               gcggcgGGCGCGGCACCACCAAACAAGGCCAACGCGAGCTGGGGGGGGATGTAAATTCTTTTTAAATGGGTAATGGGTCGCCTCAGTGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGa                                                                                                        <  2:207827/150‑1 (MQ=255)
                                                   cgTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCtt                                                                                                    >  1:5966363/1‑150 (MQ=255)
                                                     tGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCt                                                                                                  >  1:2389158/1‑150 (MQ=255)
                                                      gcgcgcCACAACAAAACACGGCAAACGCGCGCTGGGGGGGGTTGTAATTTCTTTTCAAATGGTTCTTGGGTCGCCTCAGTGGCCGGCGTGCGGGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                                 <  1:4185058/150‑1 (MQ=255)
                                                                accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                                       >  2:4631695/1‑150 (MQ=255)
                                                                accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCtt                                                                                       >  2:6712029/1‑150 (MQ=255)
                                                                  caAACACGCCCACCGCGCGCCGGGGGGGGCGGGAATTTCTTTTGAAACGGGTCATGGGTCGCCCCACGGGCGGGCCGGGGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTc                                                                                     <  2:1352778/150‑1 (MQ=255)
                                                                        cggCCACCGCGAGGTGGGGGGGGCTTTAATTTCTTTTTAATTGGTTCTTGGGTCCCCACCCGGGCCGGCCTGCGGGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                               <  1:1288829/149‑1 (MQ=255)
                                                                        aGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                               >  2:4501092/1‑150 (MQ=255)
                                                                           ccAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTg                                                                            >  2:2450909/1‑150 (MQ=255)
                                                                             tccgcgcgCTGGGGGGGGCTGTAATTTCGTTTTAAACGGGTCTTGGGTCGCCTCACGGGCCGGCGTGCGCGGTGTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTa                                                                          <  2:6547264/148‑1 (MQ=255)
                                                                             aaCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTa                                                                          >  2:686961/1‑150 (MQ=255)
                                                                             aaCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCGGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCGTTCCGGCCCCTGTa                                                                          >  1:1390208/1‑150 (MQ=255)
                                                                              aCGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCAGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTAc                                                                         >  1:3685144/1‑150 (MQ=255)
                                                                               cgcgAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCATCTGGGTGCCGCTTTCCGGCCCCAGTAca                                                                        >  1:6546802/1‑149 (MQ=255)
                                                                                 cgcgctggggggggCTGGAAATTCTTTTGAAATGGGACATGGGCCGCCTCACTGGCCTGCCTGCGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                                      <  1:4521126/144‑1 (MQ=255)
                                                                                  gAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCg                                                                     >  2:3676706/1‑150 (MQ=255)
                                                                                        ggCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAcgcg                                                               >  2:5965585/1‑150 (MQ=255)
                                                                                         ggggggGTGGAAATGCGTTTGAACGGGGTCTGGGGTCGCCGCCCGGGCGGGCCGGCGCGGTCGGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                              <  2:6073263/148‑1 (MQ=255)
                                                                                                 ggAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGTGTTTTCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAAGCGAGCAGgcgc                                                      >  1:4313647/1‑150 (MQ=255)
                                                                                                  gAAATTCTTTTGAAAAGGGTCATGGGTCGCCTCACGGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgct                                                     <  2:5225773/150‑1 (MQ=255)
                                                                                                  caaaTTGTTTTTAAACGGGTCTGGGGGCGCCCCCGGGGCGGGCCGGCGGGGGCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgct                                                     <  1:5163524/149‑1 (MQ=255)
                                                                                                   aaaTTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCGGCGCGGTCTGCGCAGGCACCGCGAAGTTGTTGCGTTTTCAGAGCTTGGCTTAGACCATCTGGGTGCCGCGTTCCGGCCCCAGTACCACGAACGCGAGCACgcgcgg                                                    >  1:5834662/1‑148 (MQ=255)
                                                                                                   aaaTTCGTTGGAAACGGGCCCGGGGGCGCCCCCGGGGCCGGCCGGCGCGGGCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctg                                                    <  1:6706462/150‑1 (MQ=255)
                                                                                                      ttCGTTTGAAACGGGCCTGGGGTCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctggct                                                 <  2:811518/150‑1 (MQ=255)
                                                                                                       tCGTTTTAAAAGGGTCTTGGGTCGCCCCACGGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctggctg                                                <  2:4761516/150‑1 (MQ=255)
                                                                                                          tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGTGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                             <  2:2725226/150‑1 (MQ=255)
                                                                                                          tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGTGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                             <  2:6190593/150‑1 (MQ=255)
                                                                                                          tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                             >  2:1825842/1‑150 (MQ=255)
                                                                                                          tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                             >  1:6247284/1‑150 (MQ=255)
                                                                                                          tttGAAAAGGGTCATGGGTAGCCTCAATGGCCTGACTTAGCGGTCAGCGCAGGCACGGAGAAGTTCTTTCGTTTTCAGAGATTGGCTTATAACTAATTGGTGCCGCTTTCAGGCCACTGTACCACGAACGCGAGCAGTCGCTGTCTGccc                                             >  2:4979651/1‑150 (MQ=255)
                                                                                                           ttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCTAGCAGGCGCTGGCTGGccc                                            >  1:2711995/1‑150 (MQ=255)
                                                                                                            tGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa                                           >  2:2849202/1‑150 (MQ=255)
                                                                                                             gAAATGGGCCTTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAg                                          <  2:4064937/150‑1 (MQ=255)
                                                                                                              aaaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCAGTACCTCGAACGCGAGCCGGCGCTGGCTGGCCCAGt                                         >  1:4363610/1‑150 (MQ=255)
                                                                                                               ccaGGGTCCGGGGGCGCCTCCCGGGCCGGCCGGCGCGGTCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTc                                        <  1:4913717/148‑1 (MQ=255)
                                                                                                               aatgggTCTTGGGGCGCCACCGGGGCGTGCCTGCGGGGTGTGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTc                                        <  1:1551600/147‑1 (MQ=255)
                                                                                                               aaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTc                                        >  1:5606533/1‑150 (MQ=255)
                                                                                                                  gggTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                     <  2:5783112/150‑1 (MQ=255)
                                                                                                                  gggTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCGTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                     >  2:792017/1‑150 (MQ=255)
                                                                                                                    gTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                                   >  2:4610436/1‑150 (MQ=255)
                                                                                                                    gTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGt                                   >  2:1423955/1‑150 (MQ=255)
                                                                                                                     tCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTc                                  >  2:6411281/1‑150 (MQ=255)
                                                                                                                       ttGGGTCGCCTCCCGGGCCGGCCTGCGCGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCtt                                <  2:5540250/149‑1 (MQ=255)
                                                                                                                        tcgtgggCCTCCCGGGCCGGCCGGGGGGTTCGGCGCCGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTc                               <  1:2456500/144‑1 (MQ=255)
                                                                                                                          ggTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                                >  1:5681547/1‑83 (MQ=255)
                                                                                                                          ggTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                                <  2:5681547/83‑1 (MQ=255)
                                                                                                                            tCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGAACTTCTGGGTGCCGCTTTGCGGCCCCTGTACCACGAAAGCGAGCAGGCGCTGGCTGGACCAGTCCTCGTCTTCgcgc                           >  1:6618337/1‑150 (MQ=255)
                                                                                                                                   acgggCCTGCCTGCGGGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGa                    <  1:562151/147‑1 (MQ=255)
                                                                                                                                   acgggCCTGCCGGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGa                    <  2:5409288/147‑1 (MQ=255)
                                                                                                                                   aCTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGa                    >  2:3263423/1‑150 (MQ=255)
                                                                                                                                    cTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCGGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCACGTCTTCGCGCATCTGGaa                   >  2:2098506/1‑150 (MQ=255)
                                                                                                                                     cgggcgggcctgcGCGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaa                   <  2:534430/143‑1 (MQ=255)
                                                                                                                                       gcctgcctGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGAc                >  2:6294863/1‑150 (MQ=255)
                                                                                                                                         cgtgcctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTAGACCATCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAACCGACCAGGCGCTGGCTGGCCCAGTCCTCGTCTTAGCGCCTATTGAAGACg               >  2:296795/3‑150 (MQ=255)
                                                                                                                                         ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa              >  1:256568/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)
                                                                                                                                              ctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCg          >  1:167341/1‑150 (MQ=255)
                                                                                                                                              ctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCg          >  1:5492899/1‑150 (MQ=255)
                                                                                                                                              ctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCg          >  1:5296002/1‑150 (MQ=255)
                                                                                                                                                    ggTCTGCGCTGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTc    <  1:902351/150‑1 (MQ=255)
                                                                                                                                                      tCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCgg  >  1:4964652/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCTTGAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCGG  >  NC_002947/1469005‑1469303

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: