Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,281,581 G→T 60.9% D48Y (GAT→TAT)  nrfE → heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,281,5810GT60.9% 120.5 / inf 323D48Y (GAT→TAT) nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (55/71);  new base T (102/94);  total (158/165)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 324 total aligned reads displayed.

TCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTAT  >  CP009273/4281435‑4281730
                                                                                                                                                  |                                                                                                                                                     
     gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCg                                                                                                                                               <  2:100022/150‑1 (MQ=255)
     gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCg                                                                                                                                               <  1:403284/150‑1 (MQ=255)
     gtgtTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCg                                                                                                                                               >  2:1358930/1‑150 (MQ=255)
        ttGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTc                                                                                                                                            >  1:307493/1‑150 (MQ=255)
        ttGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTc                                                                                                                                            <  2:4579405/150‑1 (MQ=255)
        ttGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTc                                                                                                                                            <  1:3902317/150‑1 (MQ=255)
           aCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt                                                                                                                                         <  2:2055969/150‑1 (MQ=255)
           aCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt                                                                                                                                         <  2:4065539/150‑1 (MQ=255)
             cccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTa                                                                                                                                       >  1:3654765/1‑150 (MQ=255)
              ccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTAt                                                                                                                                      <  2:4569956/150‑1 (MQ=255)
                 ttGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTc                                                                                                                                   >  2:1476707/1‑150 (MQ=255)
                  tGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTcg                                                                                                                                  >  1:3990212/1‑150 (MQ=255)
                  tGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTcg                                                                                                                                  >  1:1784627/1‑150 (MQ=255)
                  tGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTcg                                                                                                                                  >  1:3661110/1‑150 (MQ=255)
                    aCGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTcgcg                                                                                                                                <  2:4292767/150‑1 (MQ=255)
                      ggCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCa                                                                                                                              >  1:203027/1‑150 (MQ=255)
                      ggCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCa                                                                                                                              >  1:2896532/1‑150 (MQ=255)
                        ccTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAAc                                                                                                                            <  2:3269621/150‑1 (MQ=255)
                        ccTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAAc                                                                                                                            <  2:2014645/150‑1 (MQ=255)
                          ttCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCTTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACAt                                                                                                                          <  1:1988717/150‑1 (MQ=255)
                            cgcgGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAg                                                                                                                        <  2:3017647/150‑1 (MQ=255)
                             gcgGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGc                                                                                                                       >  1:4377122/1‑150 (MQ=255)
                              cgGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCt                                                                                                                      >  2:1614657/1‑150 (MQ=255)
                                 gAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCGTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtata                                                                                                                   >  1:2148505/1‑150 (MQ=255)
                                           gcgcTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTc                                                                                                         <  1:875622/150‑1 (MQ=255)
                                            cgcTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCg                                                                                                        <  2:2838999/150‑1 (MQ=255)
                                                 gCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGa                                                                                                   <  2:3796748/150‑1 (MQ=255)
                                                 gCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGTGCAACATAGCTATAGCCTGCTGTCGTGGGa                                                                                                   <  1:4290793/150‑1 (MQ=255)
                                                 gCCTGCCATGATGCAACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGa                                                                                                   >  2:288315/1‑150 (MQ=255)
                                                  ccTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGaa                                                                                                  <  1:500504/150‑1 (MQ=255)
                                                     gCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTc                                                                                               >  1:1068715/1‑150 (MQ=255)
                                                       cATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCaa                                                                                             <  2:3686846/150‑1 (MQ=255)
                                                          gttgCGACTCACTTGCATCGGCATTCTGGCGCTGTTCGCGCTCCTGCTGCTAGCCTTTGGCGTACTGACTTATTGTTTTCTCATCAGAGATTTCTCGTTCATTTTTGTCTCGCAACATAGCTTTAGCATGATGTCGTTGGTACTcttgct                                                                                          <  2:478376/148‑6 (MQ=255)
                                                          gatgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCt                                                                                          >  2:3216555/1‑150 (MQ=255)
                                                          gatgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCt                                                                                          >  2:2428418/1‑150 (MQ=255)
                                                          gatgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCt                                                                                          <  1:3431760/150‑1 (MQ=255)
                                                          gatgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCCATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCt                                                                                          >  2:3817624/1‑150 (MQ=255)
                                                           atgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTg                                                                                         >  2:2268846/1‑150 (MQ=255)
                                                              cGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCa                                                                                      <  1:3732683/150‑1 (MQ=255)
                                                              cGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCa                                                                                      <  2:3920217/150‑1 (MQ=255)
                                                              cGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCa                                                                                      <  2:3377822/150‑1 (MQ=255)
                                                              cGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCa                                                                                      >  1:1876856/1‑150 (MQ=255)
                                                                aCTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGc                                                                                    >  2:785935/1‑150 (MQ=255)
                                                                 cTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCg                                                                                   >  1:3370651/1‑150 (MQ=255)
                                                                 cTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCg                                                                                   <  2:2540831/150‑1 (MQ=255)
                                                                  tCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCgg                                                                                  <  2:4537693/150‑1 (MQ=255)
                                                                     cTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGgtgt                                                                               <  1:4420625/150‑1 (MQ=255)
                                                                       tGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTgg                                                                             <  1:3435465/150‑1 (MQ=255)
                                                                       tGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTgg                                                                             <  1:1109733/150‑1 (MQ=255)
                                                                        gCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTggg                                                                            >  1:2494548/1‑150 (MQ=255)
                                                                          aTCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggc                                                                          >  1:2359340/1‑150 (MQ=255)
                                                                          aTCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggc                                                                          >  1:3405408/1‑150 (MQ=255)
                                                                            cGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggcgg                                                                        <  1:776098/150‑1 (MQ=255)
                                                                                  tCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGa                                                                  >  1:4448741/1‑150 (MQ=255)
                                                                                  tCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGa                                                                  <  2:2958351/150‑1 (MQ=255)
                                                                                    tGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAg                                                                <  2:3119295/150‑1 (MQ=255)
                                                                                    tGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAg                                                                <  2:2212963/150‑1 (MQ=255)
                                                                                         cAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCg                                                           >  2:149893/1‑150 (MQ=255)
                                                                                         cAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCg                                                           <  2:837547/150‑1 (MQ=255)
                                                                                           gTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgct                                                         <  2:2417961/150‑1 (MQ=255)
                                                                                           gTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgct                                                         <  2:1753440/150‑1 (MQ=255)
                                                                                                 gcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgctgct                                                   <  1:2262654/150‑1 (MQ=255)
                                                                                                 gcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgctgct                                                   >  2:3053446/1‑150 (MQ=255)
                                                                                                  cTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCtt                                                  <  2:3390279/150‑1 (MQ=255)
                                                                                                     ctgctgCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTgg                                               <  2:1114013/150‑1 (MQ=255)
                                                                                                       gctgctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGt                                             >  1:2754044/1‑150 (MQ=255)
                                                                                                        ctgctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtg                                            <  2:3714315/150‑1 (MQ=255)
                                                                                                        ctgctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGATCGCTGCTGCTTTGGGtg                                            >  2:2346833/1‑150 (MQ=255)
                                                                                                         tgctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgc                                           >  1:3029834/1‑150 (MQ=255)
                                                                                                         tgctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgc                                           >  2:941639/1‑150 (MQ=255)
                                                                                                          gctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgct                                          <  1:1741801/150‑1 (MQ=255)
                                                                                                           ctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctg                                         >  2:1538412/1‑150 (MQ=255)
                                                                                                           ctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATGTATGTCGCGCCACATAGCGATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCCTGAAGGTTCGCTGCTGCGTTGGGagctg                                         >  2:1777553/1‑150 (MQ=255)
                                                                                                            tCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAACCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgc                                        <  1:2125414/150‑1 (MQ=255)
                                                                                                            tCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgc                                        <  2:1576619/150‑1 (MQ=255)
                                                                                                             cGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgct                                       <  1:3198110/150‑1 (MQ=255)
                                                                                                               ccTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgctgc                                     >  1:1231820/1‑150 (MQ=255)
                                                                                                                       gTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGcc                             >  2:4455351/1‑150 (MQ=255)
                                                                                                                        tACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCt                            >  2:2136205/1‑150 (MQ=255)
                                                                                                                        tACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCt                            <  1:1526500/150‑1 (MQ=255)
                                                                                                                        tACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCt                            <  2:4030094/150‑1 (MQ=255)
                                                                                                                           tGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGa                         >  1:3519729/1‑150 (MQ=255)
                                                                                                                              cGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcg                      <  1:1040632/150‑1 (MQ=255)
                                                                                                                              cGTATTGTTTTCTCATCAGCTATGTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcg                      <  2:2485643/150‑1 (MQ=255)
                                                                                                                              cGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcg                      >  1:575606/1‑150 (MQ=255)
                                                                                                                              cGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcg                      <  2:4349503/150‑1 (MQ=255)
                                                                                                                               gTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcgc                     <  1:2596720/150‑1 (MQ=255)
                                                                                                                                tATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcgcg                    >  2:3720368/1‑150 (MQ=255)
                                                                                                                                  ttGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCt                  <  1:2850824/150‑1 (MQ=255)
                                                                                                                                    gTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGt                >  1:2331147/1‑150 (MQ=255)
                                                                                                                                     ttttCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGtt               <  1:3570922/150‑1 (MQ=255)
                                                                                                                                      tttCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGttt              <  2:1217013/150‑1 (MQ=255)
                                                                                                                                       ttCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTg             <  2:3956806/150‑1 (MQ=255)
                                                                                                                                        tctcATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGc            >  2:2669443/1‑150 (MQ=255)
                                                                                                                                            atcaGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTgg        >  2:1881328/1‑150 (MQ=255)
                                                                                                                                            atcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTgg        >  2:638806/1‑150 (MQ=255)
                                                                                                                                               aGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCAt     >  1:1692181/1‑150 (MQ=255)
                                                                                                                                                gCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCAtt    <  2:1247088/150‑1 (MQ=255)
                                                                                                                                                 cGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTa   <  2:2263819/150‑1 (MQ=255)
                                                                                                                                                 cGATTTCTCGGTCATTTATGTCGCGCAACATAGCTAGAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTa   <  2:4276173/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTAT  >  CP009273/4281435‑4281730

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: