Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 627,794 G→A 5.3% D53N (GAC→AAC)  ybdD → DUF466 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273627,7940GA 5.3% 1162.1 / 32.7 321D53N (GAC→AAC) ybdDDUF466 family protein
Reads supporting (aligned to +/- strand):  ref base G (174/130);  new base A (7/10);  total (181/140)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 321 total aligned reads displayed.

TCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAGG  >  CP009273/627647‑627941
                                                                                                                                                   |                                                                                                                                                   
tCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGAc                                                                                                                                                   <  1:1419746/150‑1 (MQ=255)
tCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACATACGAAGAGTTTTTCCGTGAGCGGCAGGAc                                                                                                                                                   <  2:2170355/150‑1 (MQ=255)
   cTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGAcgcg                                                                                                                                                >  2:4163141/1‑150 (MQ=255)
    tGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTTTGCCTGATTACGACAACTATGTCGAACATATGCGGGTTTACCATCCCGATCAAAAGCCGATGACCAACGAAGAGTTTTTCCGGGAGCGGTAGGAcgcgc                                                                                                                                               >  2:76201/1‑150 (MQ=255)
        aaaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTa                                                                                                                                           >  2:245001/1‑150 (MQ=255)
         aaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTAc                                                                                                                                          >  1:2879922/1‑150 (MQ=255)
              cGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTgg                                                                                                                                     >  1:480025/1‑150 (MQ=255)
               ggAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGa                                                                                                                                    >  1:1754922/1‑150 (MQ=255)
               ggAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGa                                                                                                                                    >  1:2924269/1‑150 (MQ=255)
                   aaTATGTAGGTCAGGCGACGACGCAGATGATTCCTATGCCTGATAACCACAACAAACACGAACATAAGCGGCATAAACATCCCCATCAAACGCCGAAGACCTACGACGACTAGATCCGTGATCGGCAGGACGCGCGCTACGGTGGAAAAg                                                                                                                                >  2:3822546/1‑150 (MQ=255)
                        ttAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAggcggc                                                                                                                           <  1:1812063/150‑1 (MQ=255)
                        ttAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAggcggc                                                                                                                           <  1:705496/150‑1 (MQ=255)
                          aGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGgcgc                                                                                                                         >  1:570251/1‑150 (MQ=255)
                               aGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCCACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCgct                                                                                                                    <  2:2298351/150‑1 (MQ=255)
                                   ggcgAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCgctgctg                                                                                                                <  2:1054963/150‑1 (MQ=255)
                                   ggcgAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCgctgctg                                                                                                                >  2:4922817/1‑150 (MQ=255)
                                   ggcgAAGCTGATGATTGGGCTGCCTGATTACGACAACTATGTCGAACATATGCGGGGTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCgctgctg                                                                                                                <  2:3673140/150‑1 (MQ=255)
                                          cTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCAACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTCAAtt                                                                                                         >  2:3238259/1‑150 (MQ=255)
                                           tgatgaTTTGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAAttt                                                                                                        <  2:1967671/150‑1 (MQ=255)
                                           tgatgaTTGGTCTGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAAttt                                                                                                        <  2:3898306/150‑1 (MQ=255)
                                           tgatgaTTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAAttt                                                                                                        <  2:1569294/150‑1 (MQ=255)
                                            gatgatTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTg                                                                                                       >  1:1250885/1‑150 (MQ=255)
                                              tgatTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTa                                                                                                     >  2:648110/1‑150 (MQ=255)
                                               gatTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGACAAGGCGGCGCGCGCTGCTGCTAAAGTTGTAc                                                                                                    >  2:4535428/1‑150 (MQ=255)
                                               gatTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTAc                                                                                                    <  1:4474589/150‑1 (MQ=255)
                                                 ttGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGt                                                                                                  >  2:2085850/1‑150 (MQ=255)
                                                 ttGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTActt                                                                                                  <  2:1702694/150‑3 (MQ=255)
                                                  tGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTc                                                                                                 >  1:3850411/1‑150 (MQ=255)
                                                         ccTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCttt                                                                                          <  2:3493148/150‑1 (MQ=255)
                                                         ccTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCttt                                                                                          <  1:3253614/150‑1 (MQ=255)
                                                           tGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTaa                                                                                        >  1:4600780/1‑150 (MQ=255)
                                                              ttACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACg                                                                                     <  2:3577684/150‑1 (MQ=255)
                                                               tACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTATACGa                                                                                    >  2:4127437/1‑150 (MQ=255)
                                                               tACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa                                                                                    >  2:1786667/1‑150 (MQ=255)
                                                                aCGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGAt                                                                                   <  2:1287142/150‑1 (MQ=255)
                                                                 cGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGAtt                                                                                  >  1:4495391/1‑150 (MQ=255)
                                                                 cGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGAtt                                                                                  <  1:181041/150‑1 (MQ=255)
                                                                 cGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGAtt                                                                                  <  2:112773/150‑1 (MQ=255)
                                                                    cAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCa                                                                               <  2:3786871/150‑1 (MQ=255)
                                                                      aCTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACt                                                                             >  2:2228899/1‑150 (MQ=255)
                                                                      aCTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGAGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCAACt                                                                             >  2:5019409/1‑150 (MQ=255)
                                                                       cTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACtt                                                                            >  2:1686111/1‑150 (MQ=255)
                                                                         aTGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACtttt                                                                          >  2:5046749/1‑150 (MQ=255)
                                                                           gTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACtgtttt                                                                        >  2:3526471/1‑150 (MQ=255)
                                                                             cGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGACAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGa                                                                      <  2:1226271/150‑1 (MQ=255)
                                                                             cGAACATATGCGGGGTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGa                                                                      >  2:2544879/1‑150 (MQ=255)
                                                                               aaCATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAAc                                                                    <  1:1634621/150‑1 (MQ=255)
                                                                                aCATTTGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACg                                                                   <  2:3179417/150‑1 (MQ=255)
                                                                                aCATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACg                                                                   <  2:2009685/150‑1 (MQ=255)
                                                                                  atatGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCt                                                                 <  1:2429637/150‑1 (MQ=255)
                                                                                   tatGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTg                                                                <  1:2112491/150‑1 (MQ=255)
                                                                                     ttcGGGTTTACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTAAGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCTGGCTTTAAACGATTCCACTTTTTTGAACGCTGCa                                                              <  1:4678046/148‑1 (MQ=255)
                                                                                       cGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACg                                                            >  1:1088769/1‑150 (MQ=255)
                                                                                       cGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACg                                                            <  1:457632/150‑1 (MQ=255)
                                                                                        gggTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGc                                                           >  1:3159432/1‑150 (MQ=255)
                                                                                          gTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCaa                                                         >  2:2291433/1‑150 (MQ=255)
                                                                                           ttAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAc                                                        >  2:1492796/1‑150 (MQ=255)
                                                                                           ttAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAc                                                        <  1:2131773/150‑1 (MQ=255)
                                                                                             aaCCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGt                                                      >  1:4093072/1‑150 (MQ=255)
                                                                                                cATCCCGATCAACCGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATc                                                   <  2:2754737/150‑1 (MQ=255)
                                                                                                cATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATc                                                   <  1:159174/150‑1 (MQ=255)
                                                                                                 aTCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCt                                                  >  1:5027040/1‑150 (MQ=255)
                                                                                                    ccGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGgt                                               >  1:4177186/1‑150 (MQ=255)
                                                                                                     cGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTGTGAACGCTGCACGCGACGTTTCTGgtg                                              >  2:4034729/1‑150 (MQ=255)
                                                                                                       aTCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTc                                            <  1:3858540/150‑1 (MQ=255)
                                                                                                       aTCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTc                                            >  1:4225005/1‑150 (MQ=255)
                                                                                                         cAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAg                                          >  2:1758563/1‑150 (MQ=255)
                                                                                                          aaaCGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGc                                         >  1:3205827/1‑150 (MQ=255)
                                                                                                          aaaCGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGGATCTGGTGTCAGc                                         >  2:4291556/1‑150 (MQ=255)
                                                                                                           aaCGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCg                                        <  2:682117/150‑1 (MQ=255)
                                                                                                             cGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTg                                      <  1:3781638/150‑1 (MQ=255)
                                                                                                               ccGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGAc                                    >  1:2156526/1‑150 (MQ=255)
                                                                                                                  aTGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTgg                                 >  2:162090/1‑150 (MQ=255)
                                                                                                                   tGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGTTGGAAATGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGc                                >  2:341754/1‑150 (MQ=255)
                                                                                                                   tGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGc                                >  2:1167093/1‑150 (MQ=255)
                                                                                                                    gACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCa                               >  1:5110977/1‑150 (MQ=255)
                                                                                                                    gACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCa                               >  1:3008650/1‑150 (MQ=255)
                                                                                                                     aCCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAg                              >  1:1441663/1‑150 (MQ=255)
                                                                                                                     aCCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAg                              <  2:29039/150‑1 (MQ=255)
                                                                                                                     aCCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAg                              >  1:3430557/1‑150 (MQ=255)
                                                                                                                       cTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGt                            >  1:2018026/1‑150 (MQ=255)
                                                                                                                         aCGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaa                          >  1:352005/1‑150 (MQ=255)
                                                                                                                              gagTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaatgaat                     <  1:3364097/150‑1 (MQ=255)
                                                                                                                              gagTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATACCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaatgaat                     <  1:2620443/150‑1 (MQ=255)
                                                                                                                               agTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaatgaatg                    >  1:1341521/1‑150 (MQ=255)
                                                                                                                                gTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATgg                   >  2:4490521/1‑150 (MQ=255)
                                                                                                                                 tttttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTGCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGa                  <  1:3751478/150‑1 (MQ=255)
                                                                                                                                 tttttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGa                  <  1:2229025/150‑1 (MQ=255)
                                                                                                                                  ttttCCGTGAGCGGCTGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTTAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGAc                 <  2:3409074/150‑1 (MQ=255)
                                                                                                                                    ttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGATGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTc               <  1:5227643/150‑1 (MQ=255)
                                                                                                                                      ccGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCa             >  2:2873337/1‑150 (MQ=255)
                                                                                                                                      ccGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCa             >  2:221389/1‑150 (MQ=255)
                                                                                                                                       cGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCAc            <  2:2364081/150‑1 (MQ=255)
                                                                                                                                        gTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCAcc           <  1:4015301/150‑1 (MQ=255)
                                                                                                                                          gAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCgg         <  2:2658655/150‑1 (MQ=255)
                                                                                                                                          gAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCgg         <  2:5214992/150‑1 (MQ=255)
                                                                                                                                           aGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGa        <  2:3623545/150‑1 (MQ=255)
                                                                                                                                           aGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGa        >  2:4794273/1‑150 (MQ=255)
                                                                                                                                            gCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAACCGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGAc       <  2:955115/150‑1 (MQ=255)
                                                                                                                                               gCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCa    <  2:4619852/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAGG  >  CP009273/627647‑627941

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