Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,281,581 G→T 43.0% 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,5810GT43.0% 259.1 / inf 269D48Y (GAT→TAT) nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (60/93);  new base T (67/49);  total (127/142)
Fisher's exact test for biased strand distribution p-value = 3.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

GTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTAT  >  CP009273/4281434‑4281730
                                                                                                                                                   |                                                                                                                                                     
gTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCAGTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGAt                                                                                                                                                     >  2:1375739/1‑150 (MQ=255)
   aaCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtc                                                                                                                                                  <  2:4395075/150‑1 (MQ=255)
    aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtct                                                                                                                                                 <  2:4165492/150‑1 (MQ=255)
     cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc                                                                                                                                                <  1:1000560/150‑1 (MQ=255)
     cGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc                                                                                                                                                <  1:1633993/150‑1 (MQ=255)
      gtggtgACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGAGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCg                                                                                                                                               <  1:2835119/150‑1 (MQ=255)
            aCCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCAtt                                                                                                                                         <  1:2394981/150‑1 (MQ=255)
              cccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTa                                                                                                                                       >  1:705119/1‑150 (MQ=255)
              cccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTa                                                                                                                                       >  2:453094/1‑150 (MQ=255)
              cccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTa                                                                                                                                       >  2:453532/1‑150 (MQ=255)
              cccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTa                                                                                                                                       >  2:453102/1‑150 (MQ=255)
              cccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTACTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTa                                                                                                                                       >  1:2986285/1‑150 (MQ=255)
               ccGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTAt                                                                                                                                      <  2:1500996/150‑1 (MQ=255)
                 gTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCGGGCGCAGTTCGCGCTCCTGCTGCTCGCCGTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGt                                                                                                                                    >  1:3905690/1‑150 (MQ=255)
                   tGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTcg                                                                                                                                  <  2:3620957/150‑1 (MQ=255)
                   tGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTcg                                                                                                                                  <  2:1052390/150‑1 (MQ=255)
                    gACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTcgc                                                                                                                                 >  2:4169873/1‑150 (MQ=255)
                     aCGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTcgcg                                                                                                                                >  1:540016/1‑150 (MQ=255)
                       ggCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCa                                                                                                                              >  2:3275127/1‑150 (MQ=255)
                       ggCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCa                                                                                                                              >  1:776604/1‑150 (MQ=255)
                       ggCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCa                                                                                                                              >  1:1506779/1‑150 (MQ=255)
                          cTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACa                                                                                                                           <  2:3056333/150‑1 (MQ=255)
                            tCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATa                                                                                                                         <  1:3762882/150‑1 (MQ=255)
                            tCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATa                                                                                                                         <  1:224981/150‑1 (MQ=255)
                               cgGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCt                                                                                                                      >  1:4228652/1‑150 (MQ=255)
                                gggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCta                                                                                                                     <  1:2972310/150‑1 (MQ=255)
                                 ggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtat                                                                                                                    <  2:1404004/150‑1 (MQ=255)
                                 ggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtat                                                                                                                    <  2:2266837/150‑1 (MQ=255)
                                 ggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtat                                                                                                                    <  2:1404156/150‑1 (MQ=255)
                                 ggAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCtat                                                                                                                    >  2:2951502/1‑150 (MQ=255)
                                    gTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGc                                                                                                                 <  1:132346/150‑1 (MQ=255)
                                     tGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAgtc                                                                                                                <  1:3270455/150‑3 (MQ=255)
                                       cGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCctg                                                                                                              <  2:1422151/150‑1 (MQ=255)
                                       cGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCctg                                                                                                              <  2:2466084/150‑1 (MQ=255)
                                       cGGTGGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCctg                                                                                                              <  2:3961835/150‑1 (MQ=255)
                                          ttGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCctgctg                                                                                                           >  2:268755/1‑150 (MQ=255)
                                          ttGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCctgctg                                                                                                           <  1:854586/150‑1 (MQ=255)
                                          ttGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCctgctg                                                                                                           <  1:854562/150‑1 (MQ=255)
                                             cgcTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCg                                                                                                        <  2:3716144/150‑1 (MQ=255)
                                              gcTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGt                                                                                                       <  2:1588925/150‑1 (MQ=255)
                                                  gCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGa                                                                                                   <  2:9144/150‑1 (MQ=255)
                                                  gCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGa                                                                                                   <  2:9240/150‑1 (MQ=255)
                                                    cTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAAc                                                                                                 >  2:2793831/1‑150 (MQ=255)
                                                      gCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTc                                                                                               >  1:4508992/1‑150 (MQ=255)
                                                      gCCATCATGCGACTCACTTGCATCGGCATTCAGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTAGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGGCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTc                                                                                               >  1:32156/1‑150 (MQ=255)
                                                       ccATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa                                                                                              <  1:2266431/150‑1 (MQ=255)
                                                           gatgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCt                                                                                          <  2:3798670/150‑1 (MQ=255)
                                                           gatgCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCt                                                                                          <  2:3797737/150‑1 (MQ=255)
                                                                gACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAg                                                                                     >  2:1905289/1‑150 (MQ=255)
                                                                gACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAg                                                                                     >  1:283626/1‑150 (MQ=255)
                                                                  cTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCg                                                                                   <  1:1481603/150‑1 (MQ=255)
                                                                    cACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGgt                                                                                 >  1:370626/1‑150 (MQ=255)
                                                                      cTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGgtgt                                                                               <  2:1711629/150‑1 (MQ=255)
                                                                        tGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTgg                                                                             >  2:2175538/1‑150 (MQ=255)
                                                                         gCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTggg                                                                            >  1:1377052/1‑150 (MQ=255)
                                                                           aTCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggc                                                                          >  2:2497774/1‑150 (MQ=255)
                                                                           aTCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggc                                                                          <  2:803586/150‑1 (MQ=255)
                                                                           aTCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggc                                                                          >  1:4252194/1‑150 (MQ=255)
                                                                            tCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGACCTCAAGCTGGCAGCGGTGTGGggcg                                                                         >  2:3640886/1‑150 (MQ=255)
                                                                                 aTTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCAt                                                                    >  1:2554993/1‑150 (MQ=255)
                                                                                 aTTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCAt                                                                    >  1:2144675/1‑150 (MQ=255)
                                                                                   tCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGa                                                                  >  2:4274998/1‑150 (MQ=255)
                                                                                      ggCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAgg                                                               >  2:222578/1‑150 (MQ=255)
                                                                                        cgcAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGtt                                                             >  2:3923991/1‑150 (MQ=255)
                                                                                           aGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGGTCgc                                                          >  1:2950298/1‑150 (MQ=255)
                                                                                           aGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgc                                                          >  2:3455673/1‑150 (MQ=255)
                                                                                               cgcgCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgct                                                      <  1:1905645/150‑1 (MQ=255)
                                                                                               cgcgCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgct                                                      <  2:588119/150‑1 (MQ=255)
                                                                                                gcgcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgctg                                                     >  1:630665/1‑150 (MQ=255)
                                                                                                  gcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCCTCAGCTATTTCTCGGTCATTTATGTCGCGCAACATCGATATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgctgct                                                   >  1:2932925/1‑150 (MQ=255)
                                                                                                  gcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgctgct                                                   <  1:3379618/150‑1 (MQ=255)
                                                                                                  gcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgctgct                                                   >  2:4084835/1‑150 (MQ=255)
                                                                                                  gcTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCgctgctgct                                                   <  2:1468237/150‑1 (MQ=255)
                                                                                                   cTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCtt                                                  >  1:218598/1‑150 (MQ=255)
                                                                                                    tCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCttt                                                 >  2:3537744/1‑150 (MQ=255)
                                                                                                    tCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCttt                                                 >  1:891542/1‑150 (MQ=255)
                                                                                                    tCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCttt                                                 <  1:1268641/150‑1 (MQ=255)
                                                                                                    tCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCttt                                                 <  1:4221129/150‑1 (MQ=255)
                                                                                                    tCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCttt                                                 >  1:2635241/1‑150 (MQ=255)
                                                                                                        gctgctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGt                                             <  2:2621932/150‑1 (MQ=255)
                                                                                                         ctgctCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtg                                            >  1:3731610/1‑150 (MQ=255)
                                                                                                               gCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgctg                                      >  1:990188/1‑150 (MQ=255)
                                                                                                               gCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGtgctgctg                                      >  2:2457690/1‑150 (MQ=255)
                                                                                                                    tGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTc                                 <  1:1096568/150‑1 (MQ=255)
                                                                                                                     ggCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTcc                                <  2:2073315/150‑1 (MQ=255)
                                                                                                                     ggCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTcc                                <  2:850545/150‑1 (MQ=255)
                                                                                                                      gCGTACTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCCACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCg                               <  2:664504/150‑1 (MQ=255)
                                                                                                                          aCTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTg                           >  2:124898/1‑150 (MQ=255)
                                                                                                                           cTGACGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTgg                          <  1:4250476/150‑1 (MQ=255)
                                                                                                                             gACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAg                        >  1:1228847/1‑150 (MQ=255)
                                                                                                                              aCGTATTGTTTTCTCATCAGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgc                       >  1:3182368/1‑150 (MQ=255)
                                                                                                                                 tATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAgcgcg                    >  1:2091936/1‑150 (MQ=255)
                                                                                                                                   ttgttgTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCt                  <  2:3161977/150‑1 (MQ=255)
                                                                                                                                           tcatcaGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCt          >  2:4308026/1‑150 (MQ=255)
                                                                                                                                            catcaGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTg         <  2:80170/150‑1 (MQ=255)
                                                                                                                                            catcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTg         <  2:4214803/150‑1 (MQ=255)
                                                                                                                                            catcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTg         <  2:4215037/150‑1 (MQ=255)
                                                                                                                                            catcaGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTg         <  2:1913443/150‑1 (MQ=255)
                                                                                                                                              tcaGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGc       <  1:2901748/150‑1 (MQ=255)
                                                                                                                                              tcaGCTATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGc       <  1:4288076/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCGCCTGGAGCGCGCTGTTTGCCTGGCATTAT  >  CP009273/4281434‑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: