Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,120,216 C→A 100% D389Y (GAT→TAT)  glcD ← glycolate oxidase subunit, FAD‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,120,2160CA97.1% 1093.4 / 17.8 315D389Y (GAT→TAT) glcDglycolate oxidase subunit, FAD‑linked
Reads supporting (aligned to +/- strand):  ref base C (7/2);  new base A (162/144);  total (169/146)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 318 total aligned reads displayed.

GATTTTTTCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAACG  >  CP009273/3120067‑3120362
                                                                                                                                                     |                                                                                                                                                  
   tttttCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATcc                                                                                                                                                 <  1:4499241/150‑1 (MQ=255)
     tttCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCgg                                                                                                                                               >  1:2638395/1‑150 (MQ=255)
     tttCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCgg                                                                                                                                               <  2:3585961/150‑1 (MQ=255)
      ttCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGc                                                                                                                                              >  1:433001/1‑150 (MQ=255)
       tctcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCa                                                                                                                                             <  2:4282089/150‑1 (MQ=255)
        ctcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCAt                                                                                                                                            >  1:2097989/1‑150 (MQ=255)
         tcGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATg                                                                                                                                           <  2:4239014/150‑1 (MQ=255)
          cGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGa                                                                                                                                          >  2:391337/1‑150 (MQ=255)
          cGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGa                                                                                                                                          <  1:3190087/150‑1 (MQ=255)
            cccGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGaaa                                                                                                                                        >  1:3408589/1‑150 (MQ=255)
             ccGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAg                                                                                                                                       <  2:2574168/150‑1 (MQ=255)
                aTGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACg                                                                                                                                    <  2:4753976/150‑1 (MQ=255)
                aTGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACg                                                                                                                                    <  1:3281145/150‑1 (MQ=255)
                  gCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGtt                                                                                                                                  <  1:225450/150‑1 (MQ=255)
                     aTGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAAGTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGc                                                                                                                               <  2:4097242/150‑1 (MQ=255)
                      tGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCa                                                                                                                              <  1:4606778/150‑1 (MQ=255)
                      tGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCa                                                                                                                              <  2:1306994/150‑1 (MQ=255)
                       gTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCaa                                                                                                                             >  1:1948423/1‑150 (MQ=255)
                         tCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAaca                                                                                                                           <  1:1266708/150‑1 (MQ=255)
                          cGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAacac                                                                                                                          >  2:1206629/1‑150 (MQ=255)
                          cGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAacac                                                                                                                          <  2:1024009/150‑1 (MQ=255)
                             cACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCTAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTa                                                                                                                       <  2:1392807/150‑1 (MQ=255)
                             cACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTa                                                                                                                       <  2:398548/150‑1 (MQ=255)
                               cTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTaaa                                                                                                                     >  2:1881859/1‑150 (MQ=255)
                                tGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAAt                                                                                                                    <  2:3684127/150‑1 (MQ=255)
                                tGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAAt                                                                                                                    >  1:4420362/1‑150 (MQ=255)
                                 gATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATc                                                                                                                   <  2:3360894/150‑1 (MQ=255)
                                    gctgcCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCata                                                                                                                >  1:3251346/1‑150 (MQ=255)
                                         cgccAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATtgct                                                                                                           <  2:253259/150‑1 (MQ=255)
                                           ccAACTACAACGCAGAGTTCGAGGATCATCCCGCCCAGCTATTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGTATTAACGGGTGCATGATGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATtgctgc                                                                                                         >  1:3868640/1‑150 (MQ=255)
                                               cTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTTCTGCGATa                                                                                                     >  2:788778/1‑150 (MQ=255)
                                                   aaCGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACg                                                                                                 <  2:1826562/150‑1 (MQ=255)
                                                     cGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACggg                                                                                               <  2:2221254/150‑1 (MQ=255)
                                                      gCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGc                                                                                              >  2:1707030/1‑150 (MQ=255)
                                                           gTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGc                                                                                         <  2:2966180/150‑1 (MQ=255)
                                                              cGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCtt                                                                                      >  2:1528348/1‑150 (MQ=255)
                                                              cGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCtt                                                                                      <  2:1423002/150‑1 (MQ=255)
                                                               gAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTc                                                                                     <  1:4479152/150‑1 (MQ=255)
                                                                  gATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAg                                                                                  <  2:1257968/150‑1 (MQ=255)
                                                                     cTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTAc                                                                               <  1:830434/150‑1 (MQ=255)
                                                                     cTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTAc                                                                               >  2:2471680/1‑150 (MQ=255)
                                                                      ttCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACg                                                                              >  1:3112253/1‑150 (MQ=255)
                                                                      ttCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACg                                                                              >  1:564548/1‑150 (MQ=255)
                                                                        cccGCCCAGCTCGTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGcc                                                                            <  1:2976885/150‑1 (MQ=255)
                                                                          cGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcag                                                                          <  2:433029/150‑1 (MQ=255)
                                                                          cGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcag                                                                          <  2:877480/150‑1 (MQ=255)
                                                                             ccAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcaggca                                                                       >  2:3369882/1‑150 (MQ=255)
                                                                              cAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCcaggcag                                                                      >  2:4185771/1‑150 (MQ=255)
                                                                                gCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAggg                                                                    >  1:4288664/1‑150 (MQ=255)
                                                                                 ctctTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAgggg                                                                   <  2:714400/150‑2 (MQ=255)
                                                                                  tctTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGgcg                                                                  <  1:1481876/150‑1 (MQ=255)
                                                                                  tctACCGCGCGGGCAAATTAACCGGGTTCGTTGGCATCGAAAAGGATTATCTGGTGCATGTAGCCAAATCCGGCATGAAAGAAGTTGGCAACACGTAAATCATATTTCTGCTATAAACGCGCAATGCCTACCAGTACGACAGGCAGGgcg                                                                  >  2:3255519/1‑150 (MQ=255)
                                                                                     tCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGa                                                               <  2:3705221/150‑1 (MQ=255)
                                                                                             ggCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGAt                                                       <  2:1948423/150‑1 (MQ=255)
                                                                                              gCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCTTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATg                                                      <  1:4045408/150‑1 (MQ=255)
                                                                                              gCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATg                                                      <  1:1842333/150‑1 (MQ=255)
                                                                                              gCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATg                                                      <  1:2002589/150‑1 (MQ=255)
                                                                                                aaaTTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGt                                                    >  2:570712/1‑150 (MQ=255)
                                                                                                aaaTTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGt                                                    >  2:3412564/1‑150 (MQ=255)
                                                                                                  aTTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGc                                                  >  1:987713/1‑150 (MQ=255)
                                                                                                      aCCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGccatc                                              >  1:1322886/1‑150 (MQ=255)
                                                                                                           gTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGc                                         <  1:2471680/150‑1 (MQ=255)
                                                                                                                 tGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCagtagt                                   >  2:3288835/1‑150 (MQ=255)
                                                                                                                    cATCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAAt                                >  1:801769/1‑150 (MQ=255)
                                                                                                                     aTCGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATc                               >  2:4608737/1‑150 (MQ=255)
                                                                                                                      tCGAAAAGGATTAACGGGAGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATcc                              <  2:2100287/150‑1 (MQ=255)
                                                                                                                       cGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCg                             >  2:3418639/1‑150 (MQ=255)
                                                                                                                       cGAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCg                             >  1:1209346/1‑150 (MQ=255)
                                                                                                                        gAAAAGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCgg                            >  2:1561348/1‑150 (MQ=255)
                                                                                                                         aaaaGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCggg                           >  1:1157715/1‑150 (MQ=255)
                                                                                                                         aaaaGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCggg                           >  2:2378976/1‑150 (MQ=255)
                                                                                                                           aaGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGga                         >  2:4511701/1‑150 (MQ=255)
                                                                                                                            aGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGgag                        >  1:3025423/1‑150 (MQ=255)
                                                                                                                            aGGATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGgag                        <  2:3505184/150‑1 (MQ=255)
                                                                                                                              gATTAACGGGTGCATGTTGCCATATGCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGAt                      <  2:3891089/150‑1 (MQ=255)
                                                                                                                              gATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGAt                      <  2:2275529/150‑1 (MQ=255)
                                                                                                                              gATTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGAt                      <  1:4392671/150‑1 (MQ=255)
                                                                                                                               aTTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATa                     >  1:2854084/1‑150 (MQ=255)
                                                                                                                               aTTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATa                     >  2:298434/1‑150 (MQ=255)
                                                                                                                               aTTAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATa                     <  2:4138755/150‑1 (MQ=255)
                                                                                                                                 tAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACg                   <  2:2872296/150‑1 (MQ=255)
                                                                                                                                 tAACGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACg                   >  1:1703295/1‑150 (MQ=255)
                                                                                                                                  aaCGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGt                  >  2:2422791/1‑150 (MQ=255)
                                                                                                                                    cGGGTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTcc                <  1:2550256/150‑1 (MQ=255)
                                                                                                                                     gggTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCt               <  1:296134/150‑1 (MQ=255)
                                                                                                                                       gTGCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAc             <  2:2910661/150‑1 (MQ=255)
                                                                                                                                         gCATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccg           >  1:179890/1‑150 (MQ=255)
                                                                                                                                          cATGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccgc          >  1:72120/1‑150 (MQ=255)
                                                                                                                                           aTGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccgcc         >  1:139437/1‑150 (MQ=255)
                                                                                                                                            tGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGTAGTAGTAATCCGGGGAGATACGTCCTAccgccg        <  2:672816/150‑1 (MQ=255)
                                                                                                                                            tGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccgccg        <  2:381636/150‑1 (MQ=255)
                                                                                                                                            tGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccgccg        >  2:2184250/1‑150 (MQ=255)
                                                                                                                                            tGTTGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTAccgccg        <  2:3176468/150‑1 (MQ=255)
                                                                                                                                              ttGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCggg      <  1:2396919/150‑1 (MQ=255)
                                                                                                                                               tGCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGa     <  1:3588002/150‑1 (MQ=255)
                                                                                                                                                gCCATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGaa    >  1:3798835/1‑150 (MQ=255)
                                                                                                                                                 ccATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAAc   <  2:86515/150‑1 (MQ=255)
                                                                                                                                                 ccATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAAc   <  2:3209630/150‑1 (MQ=255)
                                                                                                                                                  cATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAACg  >  1:2789406/1‑150 (MQ=255)
                                                                                                                                                  cATATCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAACg  >  1:3438012/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GATTTTTTCTCGCCCGATGCCATGTTCGCCACTGATGCTGCCGCCAACTTCAACGCAGAGTTCGAGGATCTTCCCGCCCAGCTCTTCCGCGCGGGCAAATTCACCGGGTTCGTTGGCATCGAAAAGGATTAACGGGTGCATGTTGCCATCTCCGGCATGAAAGACGTTGGCAACACGTAAATCATATTGCTGCGATAAACGGGCAATGCCTTCCAGTACGCCAGGCAGGGCGCGACGCGGGATGGTGCCATCCATGCAGTAGTAATCCGGGGAGATACGTCCTACCGCCGGGAACG  >  CP009273/3120067‑3120362

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: ‑