Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,046,172 A→T 37.5% K256N (AAA→AAT hemN → coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,046,1720AT37.5% 270.7 / inf 275K256N (AAA→AAThemNcoproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Reads supporting (aligned to +/- strand):  ref base A (90/82);  new base T (47/56);  total (137/138)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 283 total aligned reads displayed.

AACATCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTG  >  CP009273/4046023‑4046319
                                                                                                                                                     |                                                                                                                                                   
aaCATCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCaaa                                                                                                                                                     <  1:3762290/150‑1 (MQ=255)
aaCATCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCaaa                                                                                                                                                     >  1:5045449/1‑150 (MQ=255)
 aCATCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAg                                                                                                                                                    >  1:4058347/1‑150 (MQ=255)
  cATCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGa                                                                                                                                                   >  2:3793759/1‑150 (MQ=255)
    tCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAatgatg                                                                                                                                                 <  1:479630/150‑1 (MQ=255)
       aCCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTg                                                                                                                                              >  2:2132113/1‑150 (MQ=255)
           gATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCt                                                                                                                                          >  1:4311683/1‑150 (MQ=255)
            cctgaCGGCCCGCCGACACAGACGCCGGAGAGTGTCGCCTGTACCCTGACACGTGTGGCGGAGCTGCACCCCGATAGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTg                                                                                                                                         <  2:3629842/146‑1 (MQ=255)
            aTTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTg                                                                                                                                         <  2:4760434/150‑1 (MQ=255)
            aTTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTg                                                                                                                                         <  1:2333492/150‑1 (MQ=255)
               tACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCg                                                                                                                                      >  2:3289826/1‑150 (MQ=255)
                  ggCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGt                                                                                                                                   <  2:5086528/150‑1 (MQ=255)
                  ggCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGt                                                                                                                                   >  1:1354694/1‑150 (MQ=255)
                       gCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCgca                                                                                                                              <  2:2715423/150‑1 (MQ=255)
                         cGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACTTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCgcagc                                                                                                                            >  1:2808893/1‑150 (MQ=255)
                          gAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCgcagca                                                                                                                           >  1:2156587/1‑150 (MQ=255)
                          gAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCgcagca                                                                                                                           >  1:1306705/1‑150 (MQ=255)
                              cAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCaaaaa                                                                                                                       >  2:4876668/1‑150 (MQ=255)
                               aGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAAc                                                                                                                      <  2:1836077/150‑1 (MQ=255)
                               aGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAAc                                                                                                                      <  2:4359772/150‑1 (MQ=255)
                                gACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACt                                                                                                                     <  2:4955516/150‑1 (MQ=255)
                                gACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACt                                                                                                                     >  2:3830841/1‑150 (MQ=255)
                                gACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACt                                                                                                                     <  2:4814240/150‑1 (MQ=255)
                                  cGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCg                                                                                                                   <  2:2904854/150‑1 (MQ=255)
                                  cGCCGGAGAGTATCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGAAAGACTGAGAGTCATAAACAACGCGCATCTGCCCACCATTTTAGCAGCTCAGCGCAAAATCAAAGATGCAGACCTGCCGAGTCCGCAGCAAAAACTCg                                                                                                                   >  2:1701955/1‑150 (MQ=255)
                                    ccGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGat                                                                                                                 <  2:3159939/150‑1 (MQ=255)
                                      ggAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGatat                                                                                                               <  1:2239236/150‑1 (MQ=255)
                                       gagaGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATc                                                                                                              <  1:1122612/150‑1 (MQ=255)
                                         gagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCt                                                                                                            <  2:1796716/150‑1 (MQ=255)
                                          agTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTg                                                                                                           <  2:2982880/150‑1 (MQ=255)
                                           gTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGc                                                                                                          >  1:1222524/1‑150 (MQ=255)
                                            tttCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCa                                                                                                         <  2:3784100/150‑1 (MQ=255)
                                               cGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGa                                                                                                      <  1:463000/150‑1 (MQ=255)
                                               cGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGa                                                                                                      >  2:2829470/1‑150 (MQ=255)
                                               cGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGa                                                                                                      >  2:3299934/1‑150 (MQ=255)
                                                  cTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAc                                                                                                   <  1:3793313/150‑1 (MQ=255)
                                                  cTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAAc                                                                                                   >  1:3116387/1‑150 (MQ=255)
                                                    ttACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCa                                                                                                 <  2:94629/150‑1 (MQ=255)
                                                       cccTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCg                                                                                              <  2:3030276/150‑1 (MQ=255)
                                                            aaaCGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTc                                                                                         <  1:1657243/150‑1 (MQ=255)
                                                              aCGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCt                                                                                       <  2:1834885/150‑1 (MQ=255)
                                                                       ggAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATc                                                                              <  2:445634/150‑1 (MQ=255)
                                                                       ggAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATc                                                                              >  1:1393041/1‑150 (MQ=255)
                                                                         aGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCgg                                                                            >  1:3971870/1‑150 (MQ=255)
                                                                             gAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTa                                                                        <  1:2790807/150‑1 (MQ=255)
                                                                              aaCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTAt                                                                       >  2:1362811/1‑150 (MQ=255)
                                                                              aaCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTAt                                                                       >  2:1100045/1‑150 (MQ=255)
                                                                              aaCCACGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTAt                                                                       >  2:1675326/1‑150 (MQ=255)
                                                                                ccccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCa                                                                     >  1:260288/1‑150 (MQ=255)
                                                                                 cccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAg                                                                    >  1:4242649/1‑150 (MQ=255)
                                                                                  ccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGt                                                                   <  2:2347484/150‑1 (MQ=255)
                                                                                  ccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGt                                                                   <  1:1017497/150‑1 (MQ=255)
                                                                                   cGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGtt                                                                  <  1:410910/150‑1 (MQ=255)
                                                                                    gATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGttt                                                                 <  2:2387134/150‑1 (MQ=255)
                                                                                     aTCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTa                                                                >  2:1934960/1‑150 (MQ=255)
                                                                                     aTCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTa                                                                >  1:1482914/1‑150 (MQ=255)
                                                                                       cGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATc                                                              >  1:3543147/1‑150 (MQ=255)
                                                                                       cGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATc                                                              <  2:629940/150‑1 (MQ=255)
                                                                                         tCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCgg                                                            >  1:1593197/1‑150 (MQ=255)
                                                                                          cTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGt                                                           >  1:3981455/1‑150 (MQ=255)
                                                                                            gAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTAt                                                         <  1:5126329/150‑1 (MQ=255)
                                                                                            gAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATAACCTGCAGGAAACCAACGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTAt                                                         >  2:257439/1‑150 (MQ=255)
                                                                                               tgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGa                                                      >  1:86164/1‑150 (MQ=255)
                                                                                               tgtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGAGATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGa                                                      <  1:3022739/150‑1 (MQ=255)
                                                                                                gtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGAt                                                     <  1:4805998/150‑1 (MQ=255)
                                                                                                gtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGAt                                                     <  1:2707185/150‑1 (MQ=255)
                                                                                                    ttAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACt                                                 <  2:1294832/150‑1 (MQ=255)
                                                                                                        cTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGc                                             <  1:1986055/150‑1 (MQ=255)
                                                                                                         tACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGcc                                            >  2:3973622/1‑150 (MQ=255)
                                                                                                         tACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGcc                                            >  2:1173931/1‑150 (MQ=255)
                                                                                                              gcATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTcc                                       >  2:1544505/1‑150 (MQ=255)
                                                                                                                  cTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGAt                                   >  2:1738984/1‑150 (MQ=255)
                                                                                                                  cTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGAt                                   <  1:2327197/150‑1 (MQ=255)
                                                                                                                   tGCCGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATg                                  <  2:4281755/150‑1 (MQ=255)
                                                                                                                     ccGACCATTTTTGCTGCTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGAc                                <  1:3160932/150‑1 (MQ=255)
                                                                                                                          cATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGca                           <  2:229403/150‑2 (MQ=255)
                                                                                                                          cATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCt                           <  2:2881740/150‑1 (MQ=255)
                                                                                                                            tttttGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTgg                         <  2:4570325/150‑1 (MQ=255)
                                                                                                                             ttttGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGc                        >  1:963494/1‑150 (MQ=255)
                                                                                                                              tttGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCg                       >  1:1845325/1‑150 (MQ=255)
                                                                                                                                tgctgcTCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGt                     >  1:2296737/1‑150 (MQ=255)
                                                                                                                                tgctgcTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGt                     >  2:2487188/1‑150 (MQ=255)
                                                                                                                                 gctgctCAGCGCAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTg                    <  1:2471610/150‑1 (MQ=255)
                                                                                                                                  ctgctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTgg                   <  2:2470963/150‑1 (MQ=255)
                                                                                                                                   tgctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGc                  >  2:2339841/1‑150 (MQ=255)
                                                                                                                                   tgctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGc                  >  1:4881306/1‑150 (MQ=255)
                                                                                                                                   tgctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGc                  <  1:775348/150‑1 (MQ=255)
                                                                                                                                    gctCAGCGCACAATCACAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATCTCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGcc                 <  2:1736902/150‑1 (MQ=255)
                                                                                                                                      tCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCa               >  1:3562486/1‑150 (MQ=255)
                                                                                                                                      tCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCa               >  2:3573638/1‑150 (MQ=255)
                                                                                                                                       cAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAg              >  1:1654826/1‑150 (MQ=255)
                                                                                                                                        aGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGc             >  1:4518910/1‑150 (MQ=255)
                                                                                                                                          cgcAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGt           <  2:4086524/150‑1 (MQ=255)
                                                                                                                                          cgcAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGt           <  1:3423356/150‑1 (MQ=255)
                                                                                                                                            cAAAATCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGa         >  1:1739632/1‑150 (MQ=255)
                                                                                                                                             aaaaTCAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGaa        >  1:880459/1‑150 (MQ=255)
                                                                                                                                                aTCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGc     >  1:1008568/1‑150 (MQ=255)
                                                                                                                                                  cAATGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGt   <  1:4848942/150‑1 (MQ=255)
                                                                                                                                                  cAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGt   >  2:3097596/1‑150 (MQ=255)
                                                                                                                                                   aaaGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtg  >  1:2430111/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AACATCGACCTGATTTACGGCCTGCCGAAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTG  >  CP009273/4046023‑4046319

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: