Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,325,982 G→T 100% T226T (ACG→ACT)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,325,9820GT99.2% 810.0 / ‑2.3 263T226T (ACG→ACT) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (131/128);  total (133/130)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 263 total aligned reads displayed.

TATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAC  >  CP009273/1325835‑1326127
                                                                                                                                                   |                                                                                                                                                 
tatGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTAc                                                                                                                                                 >  2:5027428/1‑150 (MQ=255)
   ggTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCAc                                                                                                                                              >  2:270177/1‑150 (MQ=255)
      ttCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGcc                                                                                                                                           >  1:3462580/1‑150 (MQ=255)
      ttCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCAGGTGGTCGAGCGGGAGCGTGAAATTAAAGCGTTAGTGCCGGAAGAGTTCTGGGAAGTCGATTCCAGCACAACCACGcc                                                                                                                                           >  1:4919190/1‑150 (MQ=255)
        cgccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCAt                                                                                                                                         >  2:3989768/1‑150 (MQ=255)
         gccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATc                                                                                                                                        <  1:4149225/150‑1 (MQ=255)
         gccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATc                                                                                                                                        <  1:4281190/150‑1 (MQ=255)
          ccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGACGTCGATGCCAGCACTACCACGCCATCt                                                                                                                                       >  1:1120875/1‑150 (MQ=255)
           cgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGTAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTg                                                                                                                                      >  2:3898176/1‑150 (MQ=255)
            gctgctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTgg                                                                                                                                     >  2:2525330/1‑150 (MQ=255)
               gctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGa                                                                                                                                  >  2:1197978/1‑150 (MQ=255)
                 tATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAg                                                                                                                                >  1:1140995/1‑150 (MQ=255)
                  aTGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTAGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGc                                                                                                                               <  1:5283049/150‑1 (MQ=255)
                      aaaaaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTg                                                                                                                           >  2:2831790/1‑150 (MQ=255)
                       aaaaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTgg                                                                                                                          >  2:513472/1‑150 (MQ=255)
                          aGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGt                                                                                                                       >  2:573214/1‑150 (MQ=255)
                          aGATCGCTCGGGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGt                                                                                                                       <  1:36415/150‑1 (MQ=255)
                           gATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGATCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTAGGAAGTAGATGCCAGCACTACCACCCCATCTGGTGAAGCGTTGGCGtt                                                                                                                      >  1:168829/1‑150 (MQ=255)
                              cGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACa                                                                                                                   >  1:3221365/1‑150 (MQ=255)
                               gCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAg                                                                                                                  >  2:5009334/1‑150 (MQ=255)
                                cTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCGGGTGAAGCGTTTGCGTTACAgg                                                                                                                 >  1:4099444/1‑150 (MQ=255)
                                  cGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTg                                                                                                               <  1:3390162/150‑1 (MQ=255)
                                          gTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtcatca                                                                                                       <  2:3171193/150‑1 (MQ=255)
                                          gTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtcatca                                                                                                       <  1:2712218/150‑1 (MQ=255)
                                          gTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtcatca                                                                                                       <  2:2047036/150‑1 (MQ=255)
                                          gTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtcatca                                                                                                       >  1:2127340/1‑150 (MQ=255)
                                           tCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAg                                                                                                      <  1:2466328/150‑1 (MQ=255)
                                            cTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGa                                                                                                     >  1:2423028/1‑150 (MQ=255)
                                              gCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAAc                                                                                                   <  1:1004553/150‑1 (MQ=255)
                                                 ggTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAATCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGAc                                                                                                >  1:2317132/1‑150 (MQ=255)
                                                 ggTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGAc                                                                                                >  1:2791038/1‑150 (MQ=255)
                                                 ggTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGAc                                                                                                <  1:3252539/150‑1 (MQ=255)
                                                   tCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACaa                                                                                              <  1:2546551/150‑1 (MQ=255)
                                                         gCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGtt                                                                                        <  2:1910813/150‑1 (MQ=255)
                                                           aGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtcc                                                                                      <  1:3784577/150‑1 (MQ=255)
                                                            gTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtccg                                                                                     <  2:1029332/150‑1 (MQ=255)
                                                             tCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtccgt                                                                                    >  1:5157254/1‑150 (MQ=255)
                                                             tCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtccgt                                                                                    >  1:1056156/1‑150 (MQ=255)
                                                              cGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtccgtc                                                                                   <  1:2998452/150‑1 (MQ=255)
                                                                gTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtccgtccg                                                                                 <  1:462046/150‑1 (MQ=255)
                                                                 tGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCgg                                                                                <  2:3198829/150‑1 (MQ=255)
                                                                 tGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCgg                                                                                >  2:5208232/1‑150 (MQ=255)
                                                                     ggTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaa                                                                            >  2:3334802/1‑150 (MQ=255)
                                                                     ggTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaa                                                                            <  2:4191958/150‑1 (MQ=255)
                                                                     ggTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaa                                                                            >  2:4922387/1‑150 (MQ=255)
                                                                          gCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAg                                                                       >  1:3485488/1‑150 (MQ=255)
                                                                            cTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGaa                                                                     >  1:4901629/1‑150 (MQ=255)
                                                                            cTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGaa                                                                     >  2:1603839/1‑150 (MQ=255)
                                                                             tggtggTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAAc                                                                    >  1:224844/1‑150 (MQ=255)
                                                                              ggtggtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACa                                                                   <  1:2877658/150‑1 (MQ=255)
                                                                               gtggtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACaa                                                                  <  1:307842/150‑1 (MQ=255)
                                                                                    cGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCa                                                             <  2:906684/150‑1 (MQ=255)
                                                                                     gagcgtgagcgtGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAg                                                            >  2:2184280/1‑150 (MQ=255)
                                                                                          tgagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGc                                                       >  2:1451106/1‑150 (MQ=255)
                                                                                           gagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCg                                                      <  2:4371272/150‑1 (MQ=255)
                                                                                              cgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCAATACCACGCCATCTGGTGAATCGTTTGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAATAACAAACTCAGTCTGCGGTa                                                   >  1:4006673/1‑150 (MQ=255)
                                                                                                tgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAg                                                 <  1:3107038/150‑1 (MQ=255)
                                                                                                 gaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGt                                                <  2:3655540/150‑1 (MQ=255)
                                                                                                 gaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGt                                                >  1:4946411/1‑150 (MQ=255)
                                                                                                 gaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGt                                                <  1:5052865/150‑1 (MQ=255)
                                                                                                 gaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGt                                                >  2:3727918/1‑150 (MQ=255)
                                                                                                  aaaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTc                                               <  2:2976798/150‑1 (MQ=255)
                                                                                                   aaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTct                                              >  1:2385832/1‑150 (MQ=255)
                                                                                                    aTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctg                                             >  2:3788021/1‑150 (MQ=255)
                                                                                                    aTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctg                                             >  2:4478787/1‑150 (MQ=255)
                                                                                                    aTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctg                                             <  2:1142450/150‑1 (MQ=255)
                                                                                                     ttAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgc                                            >  2:698102/1‑150 (MQ=255)
                                                                                                      tAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTGGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgct                                           <  1:1825735/150‑1 (MQ=255)
                                                                                                       aaaGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgctg                                          <  2:2902400/150‑1 (MQ=255)
                                                                                                       aaaGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgctg                                          <  2:5280867/150‑1 (MQ=255)
                                                                                                        aaGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTgg                                         <  1:4338253/150‑1 (MQ=255)
                                                                                                            gTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGaaaa                                     <  1:3037986/150‑1 (MQ=255)
                                                                                                             ttCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAATAAACTCAGGCTGCGGTAAGTCTGCTGGaaaaa                                    <  2:1551974/150‑1 (MQ=255)
                                                                                                              tCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAg                                   >  1:616743/1‑150 (MQ=255)
                                                                                                              tCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGACAGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAg                                   >  2:2376041/1‑150 (MQ=255)
                                                                                                                    cGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCt                             <  1:1502746/150‑1 (MQ=255)
                                                                                                                     ggAAGGGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTa                            >  2:776132/1‑150 (MQ=255)
                                                                                                                       aaGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACa                          >  2:3537475/1‑150 (MQ=255)
                                                                                                                       aaGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACa                          >  1:4859512/1‑150 (MQ=255)
                                                                                                                          agTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCg                       >  1:376230/1‑150 (MQ=255)
                                                                                                                          agTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCg                       >  2:2632425/1‑150 (MQ=255)
                                                                                                                          agTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCg                       >  1:1376765/1‑150 (MQ=255)
                                                                                                                          agTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCg                       >  1:4232269/1‑150 (MQ=255)
                                                                                                                             tCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGc                    <  2:3824832/150‑1 (MQ=255)
                                                                                                                             tCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGc                    <  1:3559109/150‑1 (MQ=255)
                                                                                                                              cTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCt                   <  1:3810465/150‑1 (MQ=255)
                                                                                                                              cTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCt                   <  2:1728615/150‑1 (MQ=255)
                                                                                                                                  gAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGaa               >  1:1983494/1‑150 (MQ=255)
                                                                                                                                  gAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGaa               <  1:1315814/150‑1 (MQ=255)
                                                                                                                                  gAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGaa               >  2:3628951/1‑150 (MQ=255)
                                                                                                                                   aaGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAAc              <  2:1057483/150‑1 (MQ=255)
                                                                                                                                    aGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACg             >  2:2566413/1‑150 (MQ=255)
                                                                                                                                      tCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTg           <  2:4602761/150‑1 (MQ=255)
                                                                                                                                           gCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGAc      >  2:1942539/1‑150 (MQ=255)
                                                                                                                                            ccAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACa     >  2:4786851/1‑150 (MQ=255)
                                                                                                                                              aGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACaaa   >  2:1883172/1‑150 (MQ=255)
                                                                                                                                              aGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACaaa   >  2:1010338/1‑150 (MQ=255)
                                                                                                                                               gCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACTAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAc  <  1:75951/150‑1 (MQ=255)
                                                                                                                                               gCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAc  <  1:1452737/150‑1 (MQ=255)
                                                                                                                                               gCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAc  <  2:1321887/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAC  >  CP009273/1325835‑1326127

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