Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 721,831 A→G 4.8% L221P (CTG→CCG)  kdpB ← potassium translocating ATPase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273721,8310AG 4.8% 235.3 / 2.2 84L221P (CTG→CCG) kdpBpotassium translocating ATPase, subunit B
Reads supporting (aligned to +/- strand):  ref base A (32/46);  new base G (2/2);  total (35/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

CGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCAAAAT  >  CP009273/721688‑721978
                                                                                                                                                |                                                                                                                                                   
cgccAATAGTGGTTGGGCTCAGACCGCACAGCGGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGCCGGTTGCCAGGAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCa                                                                                                                                                >  2:1036110/1‑149 (MQ=255)
cgccAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATCGTCAGGGCAATCAGCAGAATGGTCAGGGCa                                                                                                                                                >  1:573226/1‑149 (MQ=255)
 gccAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCaa                                                                                                                                               >  2:895338/1‑149 (MQ=255)
 gccAATAGTGGTTGGGAACAGACAGACCAGCTGCGCCACCAGTAGCGTTACGCTGACTACATTGCCGCCCCACGCCGAAAACGGCCACAGCGTGGCGGTTACCAGGAAAAAGACTATAGTTAGGACAATCGACAGAATGTGCAGGGcca                                                                                                                                               >  2:534580/1‑147 (MQ=255)
      tAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCCGCCGAATGGTCCgggagacttcg                                                                                                                                          >  2:1159539/1‑141 (MQ=255)
      tAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATgtcg                                                                                                                                          >  1:70788/1‑149 (MQ=255)
      tAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCAAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGCCAAGttcg                                                                                                                                          >  1:1129759/1‑149 (MQ=255)
       aGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATttcgt                                                                                                                                         <  1:817942/149‑1 (MQ=255)
        ggggTTGGGCTCAGACAGACCAGCACGGCCACCAGTACCGTTACGCTGACGGCAGTACCGCCCCACGCGGAAAACGGCCACAGGGTGGCGGTTGCCAGTAAAAAGACGCTCGATAATCCTATCAGCAGAATGGTCAGGGCAATttcgtt                                                                                                                                        <  2:268362/147‑1 (MQ=255)
        gTGGTTGGGATCAGACCGCCCAGCAGCGCCACCGGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATttcgtt                                                                                                                                        <  2:503158/149‑1 (MQ=255)
        gTGGTTGGGACCAGACACACCACCAGCGCCACCAGCACCGGTACGCTGACTGCATTACCGCCCCCCGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATttcgtt                                                                                                                                        <  2:996393/149‑1 (MQ=255)
         tGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGGAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATttcgttc                                                                                                                                       <  1:775195/149‑1 (MQ=255)
           gctGGGAGCAGTCCGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGAATTCCCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCgg                                                                                                                                     <  2:750029/147‑1 (MQ=255)
            ttGGGATCAGACAGCCCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTCCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGc                                                                                                                                    <  1:823368/149‑1 (MQ=255)
            ttGCGATCAGACAGGCCAGCAGCGCAACCAGTCCCGCTACGTTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGc                                                                                                                                    <  2:375827/149‑1 (MQ=255)
                gATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGttt                                                                                                                                >  1:1001100/1‑149 (MQ=255)
                gATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGttt                                                                                                                                >  1:1001252/1‑149 (MQ=255)
                 aTCAGACCGGCCAGCAGGGCCACCGGTACCGTTACGCTGACTGCAGTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTGCGGCGtttt                                                                                                                               <  2:126241/149‑1 (MQ=255)
                 aTCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCTAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGtttt                                                                                                                               >  1:837889/1‑149 (MQ=255)
                  tCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTg                                                                                                                              <  1:1061983/149‑1 (MQ=255)
                   ctgactggtcatcagcGCCACCAGTATCGTTACGCTGACTGCATTACCGCCCCAGGCGGAAAACGGCCACAGCGTGGGGGTTGCCGGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGc                                                                                                                             <  1:421909/137‑1 (MQ=255)
                   cagacagaCCAGCAGCGCCACCAGTACCGTTACTCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGc                                                                                                                             <  1:495373/149‑1 (MQ=255)
                       cagacCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGGCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCg                                                                                                                         <  1:241611/149‑1 (MQ=255)
                               cagcGGCACCAGGACCGTTACGCTGACTGCATTACCGCACCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTgcgc                                                                                                                 <  1:988209/149‑1 (MQ=255)
                                  cgcCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAACGACGAAAGTCAGGGCAATCAGCAGAAGGATCAGGGCAATTTCGTGCGGCGTGTTCCGTCGCAGGGCGCCtt                                                                                                              >  2:136452/1‑149 (MQ=255)
                                     caccaGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTcca                                                                                                           <  2:1210484/149‑1 (MQ=255)
                                         aGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCGGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAt                                                                                                       >  2:395242/1‑149 (MQ=255)
                                         aGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCGGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAt                                                                                                       >  2:395227/1‑149 (MQ=255)
                                         aGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAt                                                                                                       >  1:1149436/1‑149 (MQ=255)
                                            aCCGGTACGCTGACGGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATCGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATcgc                                                                                                    <  2:37222/149‑1 (MQ=255)
                                              cGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGa                                                                                                  <  2:897021/149‑1 (MQ=255)
                                                ttgcGCTGACGGCGTTAGCGCGCTAGGGGGAAACCGGCCACAGGGTGGCGGTTGCGAGTAAAAAGACGATAGTCAGGGCAATCGGCAGAATGGTCAGGGCAATTTGGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGatc                                                                                                <  2:1080090/146‑1 (MQ=255)
                                                ttACGCTGGCTGCATCACCGCCCCACGCGGAAAACGGCCACAGCGGGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGatc                                                                                                <  2:736886/149‑1 (MQ=255)
                                                    gCTGACTGCATTACCGCCCCACGCGGAAAACGGCCCACCTGTTTCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATcc                                                                                            <  1:341818/149‑1 (MQ=255)
                                                     cTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAACGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGGTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCg                                                                                           >  1:368192/1‑149 (MQ=255)
                                                     cTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCg                                                                                           <  2:785585/149‑1 (MQ=255)
                                                     cTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGCCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTTCGTCGCTGTGCGCCTTCCACCATCGCGATCCTCCg                                                                                           >  1:1152793/1‑149 (MQ=255)
                                                      tGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGa                                                                                          <  1:799388/149‑1 (MQ=255)
                                                      tGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGa                                                                                          <  1:1174908/149‑1 (MQ=255)
                                                            cATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGa                                                                                    <  1:775871/149‑1 (MQ=255)
                                                            cATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGa                                                                                    >  1:183308/1‑149 (MQ=255)
                                                            cATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGa                                                                                    >  2:171181/1‑149 (MQ=255)
                                                             aTTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGaa                                                                                   >  2:401984/1‑149 (MQ=255)
                                                               tACCGCTCCACGAGGAAAACGGCCAAAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCCTCCACCATCGCGATCATCCGATCCAGAAAt                                                                                 <  2:55464/149‑1 (MQ=255)
                                                                 ccGCCCCACGCGGAAAACGGCGACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGt                                                                               <  1:860223/149‑1 (MQ=255)
                                                                 ccGCCCCACGCGGAAAACGGCCACAGGGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGt                                                                               <  2:921154/149‑1 (MQ=255)
                                                                  cGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGGTAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAAGGtc                                                                              >  1:657608/1‑149 (MQ=255)
                                                                  cGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGGGCCTTCCACCATCGCGATCATCCGATCCAGCAATGtc                                                                              >  1:499763/1‑149 (MQ=255)
                                                                        aCGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGccc                                                                        >  2:1145199/1‑148 (MQ=255)
                                                                        aCGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCg                                                                        >  2:1090288/1‑149 (MQ=255)
                                                                              aaaaCGGCCACAGTGTGGCGGTTGCCAGTAAAAAGACGATCGTAAGGGCAATCAGCAGAATGGTCAGGGCCATTTCGTTCGGCGTTCTGCGTCGCAGTGCGCCTCCCACCATCTCGATCATCCGATCTAGAAACGTCTCGCCGGGGTTa                                                                  <  2:182006/149‑1 (MQ=255)
                                                                              aaaaCGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCAACATCGCGTTCATCCGATCCAGAAATGTCTCGCCGGGGTTa                                                                  <  2:1191378/149‑1 (MQ=255)
                                                                               aaaCGGCCGCAGCGTGGGGGTTGCCAGTGAAAAGTCGGTGGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTCCGGCGGTTGGCGTCGCTGTGCGCCTGCCACCAGCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTaa                                                                 <  2:108923/149‑1 (MQ=255)
                                                                                aaCGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCGGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAAc                                                                <  1:395242/149‑1 (MQ=255)
                                                                                aaCGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCGGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAAc                                                                <  1:395227/149‑1 (MQ=255)
                                                                                   ggCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCt                                                             <  2:401010/149‑1 (MQ=255)
                                                                                         aGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTaccccc                                                       >  2:122014/1‑145 (MQ=255)
                                                                                             gTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAAt                                                    <  2:592682/149‑1 (MQ=255)
                                                                                                 cGGTTGCCAGTAAATAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGAGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCCCGCCGGGGTTAACGCTACACTCAATCAcc                                                <  1:755609/149‑1 (MQ=255)
                                                                                                           tAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCTGTCag                                      >  2:1177171/1‑149 (MQ=255)
                                                                                                           tAAAAAGACGAGAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGCTACAAAAACATCACTCGTAGTTTTACGCTACACTCGATCACCAGCCAGTcgg                                      >  2:566323/1‑147 (MQ=255)
                                                                                                             aaaaGACGATAGTCAGGGGAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGCTGATTCGATCCAGAACTGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCagaa                                    <  1:787686/149‑1 (MQ=255)
                                                                                                              aaaGCCGAGGGTCAGGGCGATCAGCAGAATGGTCAGGGCAATGTCGTGCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCAGCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCagaaa                                   <  2:731458/149‑1 (MQ=255)
                                                                                                              aaaGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTGACGCTACACTCCATCACCAGCCAGTCagaag                                   >  1:1063841/1‑148 (MQ=255)
                                                                                                              aaaGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCagaaa                                   <  1:865705/149‑1 (MQ=255)
                                                                                                               aaGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCagaaag                                  >  1:256027/1‑149 (MQ=255)
                                                                                                                 gACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCATGTTTCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGTCAGTCagaaagaa                                <  1:846569/149‑1 (MQ=255)
                                                                                                                  aCGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAAt                               >  1:1194075/1‑149 (MQ=255)
                                                                                                                     aTAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATAcg                            <  1:401984/149‑1 (MQ=255)
                                                                                                                     aTAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATAcg                            >  1:260330/1‑149 (MQ=255)
                                                                                                                     aTAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATAcg                            >  1:598466/1‑149 (MQ=255)
                                                                                                                     aTAGTCAGGGCAAACAGCAGAATGGTCAGGGCAATTTCGTTAGGCGCTTTGCGCCCCTGAGGGCGTACCAGCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGCCAGAAAGAACAcg                            >  2:570615/1‑149 (MQ=255)
                                                                                                                          cAGGGCAATCAGCAGAATGGGCCGGGCACTTTCGTTCGGCGTTTCGCGGCGCTGTGCGCCTTCCACCATCGCGACCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCCCCAGCCAGTCAGAAAGAATACGCGTgc                       <  2:249899/149‑1 (MQ=255)
                                                                                                                          cAGGGCAATCAGCAGAATGGGCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTgc                       <  2:274845/149‑1 (MQ=255)
                                                                                                                                 aTCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGt                >  1:1189021/1‑149 (MQ=255)
                                                                                                                                     gcgcaaTGTTCAGGGCAATTTCGTTCGGCGTTTTGAGCAGCTCTGCTCCTGCCACCATCGCGCTCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCTATCACCAGTCAGTCAGAAAGAATACGCGTGCCGCCGGTGACg            <  2:871609/145‑1 (MQ=255)
                                                                                                                                       agagTGGTCAGGGCAATTTCGTTCGGCGTTTTGCGGCCCTGGGCGCCTCCCACCATCGCGACCATCCGATCCAGCAATGTCTCGCCGGGGTTAACGCTACACTCAACCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGa          <  1:234092/149‑1 (MQ=255)
                                                                                                                                       agATTGGTCAGGGCAATCTCGTTCAGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCACGATCACCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACGCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGa          <  1:1025675/149‑1 (MQ=255)
                                                                                                                                       agAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGa          <  2:492041/149‑1 (MQ=255)
                                                                                                                                        gattGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACTCTCGATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAg         <  2:5113/146‑1 (MQ=255)
                                                                                                                                        gAATGGTCAGGGCAATGTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAg         <  1:317969/149‑1 (MQ=255)
                                                                                                                                             gTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCGGAAAGAATACGCGTGCCGCCGGTGACGGAGGCaca    >  2:920306/1‑147 (MQ=255)
                                                                                                                                              tCAGGGCAATTTCGTTCGGCGTTTTGAGTCGCTGTGCGCCCTCCACCACCGCGGTCATCCGATCCAGAAAGGGCTCGCCGGGGGTAACGCTACACCCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGCCGGAGGCaaaa   >  2:555028/1‑149 (MQ=255)
                                                                                                                                               cAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTGCCGCCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCAAAAt  <  1:538343/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGCAATCAGCAGAATGGTCAGGGCAATTTCGTTCGGCGTTTTGCGTCGCTGTGCGCCTTCCACCATCGCGATCATCCGATCCAGAAATGTCTCGCCGGGGTTAACGCTACACTCAATCACCAGCCAGTCAGAAAGAATACGCGTGCCGCCGGTGACGGAGGCAAAAT  >  CP009273/721688‑721978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: