Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 583,341 T→G 100% S33R (AGT→CGT)  ybcH ← PRK09936 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273583,3410TG96.4% 1038.1 / 12.0 255S33R (AGT→CGT) ybcHPRK09936 family protein
Reads supporting (aligned to +/- strand):  ref base T (3/7);  new base G (112/133);  total (115/140)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 255 total aligned reads displayed.

GCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTAC  >  CP009273/583196‑583487
                                                                                                                                                 |                                                                                                                                                  
 cTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAATCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtctc                                                                                                                                               <  1:4785712/150‑1 (MQ=255)
 cTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtctc                                                                                                                                               <  2:5105650/150‑1 (MQ=255)
 cTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtctc                                                                                                                                               >  1:1370/1‑150 (MQ=255)
      gCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCATTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGttt                                                                                                                                          >  2:4252059/1‑150 (MQ=255)
       cccGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGtttt                                                                                                                                         <  2:2296209/150‑1 (MQ=255)
        ccGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTtg                                                                                                                                        >  1:4910533/1‑150 (MQ=255)
         cGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTtgt                                                                                                                                       >  2:3525014/1‑150 (MQ=255)
             tAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGtt                                                                                                                                   <  1:835689/150‑1 (MQ=255)
              aaaCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTg                                                                                                                                  >  1:1534161/1‑150 (MQ=255)
              aaaCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTg                                                                                                                                  >  1:2079611/1‑150 (MQ=255)
               aaCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGc                                                                                                                                 <  1:2243967/150‑1 (MQ=255)
                      aCGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataa                                                                                                                          <  1:2643274/150‑1 (MQ=255)
                      aCGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataa                                                                                                                          <  1:2766294/150‑1 (MQ=255)
                         tGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATAc                                                                                                                       <  2:3343053/150‑1 (MQ=255)
                          gcgcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATAcc                                                                                                                      <  1:4604716/150‑1 (MQ=255)
                            gcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCtt                                                                                                                    >  1:1647603/1‑150 (MQ=255)
                             cTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCttt                                                                                                                   <  2:1075814/150‑1 (MQ=255)
                              tGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTc                                                                                                                  <  1:690824/150‑1 (MQ=255)
                              tGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTc                                                                                                                  <  1:1955766/150‑1 (MQ=255)
                              tGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTc                                                                                                                  <  2:3783139/150‑1 (MQ=255)
                                  ctggctggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCg                                                                                                              <  1:1738928/150‑1 (MQ=255)
                                    ggctggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCa                                                                                                            <  2:2015156/150‑1 (MQ=255)
                                    ggctggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCa                                                                                                            >  1:1825436/1‑150 (MQ=255)
                                     gctggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCaa                                                                                                           >  2:3637615/1‑150 (MQ=255)
                                        gggTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAgg                                                                                                        >  1:259040/1‑150 (MQ=255)
                                         ggTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGa                                                                                                       <  2:1169040/150‑1 (MQ=255)
                                          gTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaa                                                                                                      >  2:4955221/1‑150 (MQ=255)
                                           tAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaa                                                                                                     >  1:4105146/1‑150 (MQ=255)
                                            aaaTGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaaa                                                                                                    >  1:1648313/1‑150 (MQ=255)
                                            aaaTGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaaa                                                                                                    >  1:4209679/1‑150 (MQ=255)
                                            aaaTGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaaa                                                                                                    >  1:3911821/1‑150 (MQ=255)
                                              aTGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAg                                                                                                  >  1:2955935/1‑150 (MQ=255)
                                                   tCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTg                                                                                             <  1:1912203/150‑1 (MQ=255)
                                                    cGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGa                                                                                            >  1:261415/1‑150 (MQ=255)
                                                     gCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGAc                                                                                           <  1:3797044/150‑1 (MQ=255)
                                                      ccGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGAcc                                                                                          <  1:310617/150‑1 (MQ=255)
                                                          aaCGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCaaaa                                                                                      <  1:4221621/150‑1 (MQ=255)
                                                              ggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAg                                                                                  <  1:4663253/150‑1 (MQ=255)
                                                              ggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAg                                                                                  >  1:4561309/1‑150 (MQ=255)
                                                              ggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAg                                                                                  >  1:4561395/1‑150 (MQ=255)
                                                                tCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAgtg                                                                                >  2:4876643/1‑150 (MQ=255)
                                                                  cATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTc                                                                              >  1:1701288/1‑150 (MQ=255)
                                                                  cATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTc                                                                              >  2:3853380/1‑150 (MQ=255)
                                                                   aTTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCa                                                                             <  2:1932067/150‑1 (MQ=255)
                                                                    ttGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAg                                                                            <  2:1741834/150‑1 (MQ=255)
                                                                     tGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGc                                                                           >  2:2939597/1‑150 (MQ=255)
                                                                        aaaaCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCact                                                                        >  2:3803274/1‑148 (MQ=255)
                                                                               gggTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGaaaa                                                                 <  1:2405037/150‑1 (MQ=255)
                                                                                  tATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGa                                                              >  2:1615254/1‑150 (MQ=255)
                                                                                   aTCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGaa                                                             <  2:4636219/150‑1 (MQ=255)
                                                                                    tCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAAc                                                            >  1:85598/1‑150 (MQ=255)
                                                                                     cGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACt                                                           <  2:3967937/150‑1 (MQ=255)
                                                                                      gAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACtt                                                          >  1:1953520/1‑150 (MQ=255)
                                                                                      gAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACtt                                                          >  1:2209272/1‑150 (MQ=255)
                                                                                       aaGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTa                                                         >  2:2448265/1‑150 (MQ=255)
                                                                                         gCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACg                                                       <  1:3782185/150‑1 (MQ=255)
                                                                                          ccTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGc                                                      >  1:1934828/1‑150 (MQ=255)
                                                                                           cTTGCAAACGTAACTGACTCATCAGCACCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACATTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCa                                                     >  2:2277110/1‑150 (MQ=255)
                                                                                            ttGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAt                                                    <  2:1370/150‑1 (MQ=255)
                                                                                            ttGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAt                                                    <  2:1865/150‑1 (MQ=255)
                                                                                              gCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAttt                                                  <  2:2960154/150‑1 (MQ=255)
                                                                                               cAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTa                                                 >  1:573371/1‑150 (MQ=255)
                                                                                                 aaCGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTAcc                                               <  1:2486089/150‑1 (MQ=255)
                                                                                                     tAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTg                                           <  2:3761802/150‑1 (MQ=255)
                                                                                                     tAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGATGTGTGGTTGCCAGATAATACCTGTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGACAATGAACTTACGCATTTACCAGTg                                           >  2:2574939/1‑150 (MQ=255)
                                                                                                       aCTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCa                                         >  1:3254516/1‑150 (MQ=255)
                                                                                                           aCTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTcc                                     >  1:3058134/1‑150 (MQ=255)
                                                                                                            cTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCa                                    <  2:166435/150‑1 (MQ=255)
                                                                                                              catcaGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAAt                                  <  1:1013156/150‑1 (MQ=255)
                                                                                                              catcaGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAAt                                  >  1:1319198/1‑150 (MQ=255)
                                                                                                                tcaGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATgg                                >  2:4696064/1‑150 (MQ=255)
                                                                                                                tcaGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATgg                                >  2:525188/1‑150 (MQ=255)
                                                                                                                 caGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGt                               >  1:4319371/1‑150 (MQ=255)
                                                                                                                 caGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGt                               <  1:4843147/150‑1 (MQ=255)
                                                                                                                 caGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGt                               >  1:3503732/1‑150 (MQ=255)
                                                                                                                  aGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTg                              <  2:4820224/150‑1 (MQ=255)
                                                                                                                   gCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTga                             <  2:1140139/150‑1 (MQ=255)
                                                                                                                    ccccTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTgag                            <  1:4515792/150‑1 (MQ=255)
                                                                                                                      ccTGCCTCTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGaa                          <  2:1472559/150‑1 (MQ=255)
                                                                                                                       cTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGaaa                         <  1:3271851/150‑1 (MQ=255)
                                                                                                                       cTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGaaa                         <  1:5119625/150‑1 (MQ=255)
                                                                                                                       cTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTACAATGGTGAGaaa                         <  2:66510/150‑1 (MQ=255)
                                                                                                                        tGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAc                        <  2:1060539/150‑1 (MQ=255)
                                                                                                                         gCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAcg                       >  2:4152556/1‑150 (MQ=255)
                                                                                                                           cACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAcgcg                     <  1:404105/150‑1 (MQ=255)
                                                                                                                            aCTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgt                    <  1:3214939/150‑1 (MQ=255)
                                                                                                                             cTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgtt                   <  2:2021202/150‑1 (MQ=255)
                                                                                                                              tGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttg                  >  1:2154371/1‑150 (MQ=255)
                                                                                                                                 gTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGttgttg               <  2:4520174/150‑1 (MQ=255)
                                                                                                                                   aTCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAACCGCGTTGTTgcg             >  2:2814680/1‑150 (MQ=255)
                                                                                                                                      ggTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGctc          <  2:3069897/150‑1 (MQ=255)
                                                                                                                                      ggTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGctc          >  2:21534/1‑150 (MQ=255)
                                                                                                                                       gTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGctct         <  2:4462320/150‑1 (MQ=255)
                                                                                                                                         aaCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCtctc        >  2:5126558/1‑149 (MQ=255)
                                                                                                                                         aaCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTcc       <  1:4352649/150‑1 (MQ=255)
                                                                                                                                         aaCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTcc       >  2:3263506/1‑150 (MQ=255)
                                                                                                                                          aCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCg      >  1:107251/1‑150 (MQ=255)
                                                                                                                                          aCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCg      >  1:358438/1‑150 (MQ=255)
                                                                                                                                              gACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTAc  >  1:1234561/1‑150 (MQ=255)
                                                                                                                                              gACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTAc  >  2:859083/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTAC  >  CP009273/583196‑583487

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