Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,881,078 G→T 34.7% G198G (GGC→GGA cysH ← phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,881,0780GT34.7% 329.6 / 188.7 216G198G (GGC→GGAcysHphosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent
Reads supporting (aligned to +/- strand):  ref base G (78/63);  new base T (38/37);  total (116/100)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 220 total aligned reads displayed.

CAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGTTTGC  >  CP009273/2880931‑2881225
                                                                                                                                                   |                                                                                                                                                   
cAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAgg                                                                                                                                                     >  2:75561/1‑148 (MQ=255)
cAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAgg                                                                                                                                                     >  1:1446421/1‑148 (MQ=255)
  atatTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGcc                                                                                                                                                   >  1:1886431/1‑148 (MQ=255)
   tatTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCa                                                                                                                                                  >  2:184687/1‑148 (MQ=255)
     tttACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATg                                                                                                                                                <  1:1944557/148‑1 (MQ=255)
      ttACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGt                                                                                                                                               >  1:1602835/1‑148 (MQ=255)
       tACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGtt                                                                                                                                              <  1:277004/148‑1 (MQ=255)
       tACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCAGAAGTGGTGAGATTTCAGGCCATGtt                                                                                                                                              >  2:3024002/1‑148 (MQ=255)
        aCCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGttt                                                                                                                                             >  1:3267668/1‑148 (MQ=255)
        aCCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGttt                                                                                                                                             <  1:1367539/148‑1 (MQ=255)
         cccTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGtttt                                                                                                                                            <  1:2018615/148‑1 (MQ=255)
         cccTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGtttt                                                                                                                                            <  1:346743/148‑1 (MQ=255)
          ccTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGttttt                                                                                                                                           <  2:1286744/148‑1 (MQ=255)
          acTTCGTGTAACCCACATTCCATTTTTAAGCCAAAGAAACGCGTTTCTTCTTGCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATTATGTGTGATATTTCAGGCCATGttttt                                                                                                                                           <  2:2382269/147‑1 (MQ=255)
           cTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTg                                                                                                                                          >  1:740825/1‑148 (MQ=255)
            ttCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGc                                                                                                                                         <  1:851482/148‑1 (MQ=255)
             tCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCa                                                                                                                                        >  1:2981877/1‑148 (MQ=255)
             tCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCa                                                                                                                                        <  2:923322/148‑1 (MQ=255)
                 gtAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTa                                                                                                                                    >  2:3322471/1‑148 (MQ=255)
                  tAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTAc                                                                                                                                   >  2:1208716/1‑148 (MQ=255)
                  tAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTAc                                                                                                                                   >  2:1803559/1‑148 (MQ=255)
                  tAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTAc                                                                                                                                   >  2:1803526/1‑148 (MQ=255)
                     cccACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGa                                                                                                                                <  2:1596151/148‑1 (MQ=255)
                       cacaTTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATa                                                                                                                              <  1:57515/148‑1 (MQ=255)
                        acaTTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTAAGTCCATGTTTTTGCAGGTACTGATaa                                                                                                                             >  2:1921800/1‑148 (MQ=255)
                         caTTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATaaa                                                                                                                            >  1:1240618/1‑148 (MQ=255)
                            tCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAg                                                                                                                         <  1:1921256/148‑1 (MQ=255)
                             cccTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGt                                                                                                                        <  1:231561/148‑1 (MQ=255)
                               cTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTc                                                                                                                      >  2:3377821/1‑148 (MQ=255)
                                       gCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATcc                                                                                                              <  1:2004292/148‑1 (MQ=255)
                                           aaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGt                                                                                                          >  2:1027882/1‑148 (MQ=255)
                                             gAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCg                                                                                                        <  1:187078/148‑1 (MQ=255)
                                              aaaCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGa                                                                                                       <  1:868096/148‑1 (MQ=255)
                                                 cgcgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATaa                                                                                                    <  2:826158/148‑1 (MQ=255)
                                                 cgcgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATaa                                                                                                    <  2:2979576/148‑1 (MQ=255)
                                                  gcgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAAt                                                                                                   <  2:2702504/148‑1 (MQ=255)
                                                      ttcttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCTATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgc                                                                                               >  1:2269013/1‑148 (MQ=255)
                                                      ttcttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgc                                                                                               >  1:2107795/1‑148 (MQ=255)
                                                      ttcttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgc                                                                                               <  2:2721910/148‑1 (MQ=255)
                                                       tcgtcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgca                                                                                              <  2:2929287/145‑1 (MQ=255)
                                                           cttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCAc                                                                                          <  2:725922/148‑1 (MQ=255)
                                                            ttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACt                                                                                         <  1:3515447/148‑1 (MQ=255)
                                                             tcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACtt                                                                                        >  2:1566877/1‑148 (MQ=255)
                                                             tcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACtt                                                                                        >  2:2813348/1‑148 (MQ=255)
                                                                  cATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAAt                                                                                   <  2:2596366/148‑1 (MQ=255)
                                                                  cATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAAt                                                                                   <  2:703957/148‑1 (MQ=255)
                                                                  cATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAAt                                                                                   <  1:2739692/148‑1 (MQ=255)
                                                                  cATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAAt                                                                                   <  1:1205504/148‑1 (MQ=255)
                                                                  cATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGAAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAAt                                                                                   >  1:1629780/1‑148 (MQ=255)
                                                                       cGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGcc                                                                              <  2:3124768/148‑1 (MQ=255)
                                                                         ggTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCAc                                                                            >  2:3291588/1‑148 (MQ=255)
                                                                           ttCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGc                                                                          >  1:2971430/1‑148 (MQ=255)
                                                                            tCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCt                                                                         >  1:1130978/1‑148 (MQ=255)
                                                                            tCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCt                                                                         >  1:2630767/1‑148 (MQ=255)
                                                                                aTTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAAt                                                                     >  2:342661/1‑148 (MQ=255)
                                                                                  ttACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTg                                                                   <  1:1049584/148‑1 (MQ=255)
                                                                                  ttACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTg                                                                   >  1:1764011/1‑148 (MQ=255)
                                                                                  ctACGGGTCGTATGGGTATCGCCCACCGATAAATATCCGTCCTCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTg                                                                   <  1:3367185/147‑1 (MQ=255)
                                                                                   tACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGc                                                                  <  1:1430862/148‑1 (MQ=255)
                                                                                       ggTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGc                                                              >  1:682573/1‑148 (MQ=255)
                                                                                         ttGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCAc                                                            <  2:1223973/148‑1 (MQ=255)
                                                                                         ttGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCAc                                                            >  1:2839133/1‑148 (MQ=255)
                                                                                          tGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCAcc                                                           <  1:1765848/148‑1 (MQ=255)
                                                                                           gTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCg                                                          <  2:2553478/148‑1 (MQ=255)
                                                                                              tGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTa                                                       <  2:2803048/148‑1 (MQ=255)
                                                                                                ggTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTaaa                                                     <  1:367468/148‑1 (MQ=255)
                                                                                                 gTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAAt                                                    <  1:2740023/148‑1 (MQ=255)
                                                                                                   aTCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTg                                                  >  2:3519997/1‑148 (MQ=255)
                                                                                                         cACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGa                                            >  1:3055915/1‑148 (MQ=255)
                                                                                                         cACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGa                                            >  1:3055991/1‑148 (MQ=255)
                                                                                                           ccGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACt                                          >  2:2225377/1‑148 (MQ=255)
                                                                                                           ccGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACt                                          >  2:2862569/1‑148 (MQ=255)
                                                                                                               tAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCg                                      >  1:1047950/1‑148 (MQ=255)
                                                                                                                aaaTATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCgg                                     >  2:153410/1‑148 (MQ=255)
                                                                                                                  atatCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGAt                                   >  2:1273778/1‑148 (MQ=255)
                                                                                                                  atatCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGAt                                   >  2:1486018/1‑148 (MQ=255)
                                                                                                                   tatCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGAtt                                  <  2:2105981/148‑1 (MQ=255)
                                                                                                                   tatCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGAtt                                  <  2:1446421/148‑1 (MQ=255)
                                                                                                                    atCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTg                                 >  2:1648928/1‑148 (MQ=255)
                                                                                                                    atCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTg                                 <  1:744469/148‑1 (MQ=255)
                                                                                                                       cTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTc                              <  2:2489261/148‑1 (MQ=255)
                                                                                                                        ttCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCa                             <  2:1398642/148‑1 (MQ=255)
                                                                                                                            tCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGgc                         >  2:399984/1‑148 (MQ=255)
                                                                                                                               cATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCa                      <  2:3200674/148‑1 (MQ=255)
                                                                                                                                 tAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAgg                    >  2:3006128/1‑148 (MQ=255)
                                                                                                                                  aaTGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGc                   <  1:1501726/148‑1 (MQ=255)
                                                                                                                                  aaTGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGc                   <  1:3401462/148‑1 (MQ=255)
                                                                                                                                   aTGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGcc                  >  2:2769611/1‑148 (MQ=255)
                                                                                                                                     gggTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAg                <  2:3267668/148‑1 (MQ=255)
                                                                                                                                      ggTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGc               >  2:3255860/1‑148 (MQ=255)
                                                                                                                                      ggTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGAACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGc               >  2:2453748/1‑148 (MQ=255)
                                                                                                                                          atatTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAAc           <  1:184687/148‑1 (MQ=255)
                                                                                                                                          atatTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAAc           <  1:554387/148‑1 (MQ=255)
                                                                                                                                           tatTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAAcc          <  2:1763994/148‑1 (MQ=255)
                                                                                                                                            atTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCa         >  2:2099864/1‑148 (MQ=255)
                                                                                                                                              ttCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAgg       <  1:399984/148‑1 (MQ=255)
                                                                                                                                              ttCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAgg       >  2:2212037/1‑148 (MQ=255)
                                                                                                                                              ttCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAgg       >  1:47514/1‑148 (MQ=255)
                                                                                                                                               tCAGTCCATTTTTTTTCAGTTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAATTTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGt      <  1:2134113/148‑1 (MQ=255)
                                                                                                                                                   gCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGTTTgc  >  1:2413467/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGTTTGC  >  CP009273/2880931‑2881225

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: