Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,881,078 G→T 35.0% 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,0780GT35.0% 249.7 / 136.9 162G198G (GGC→GGAcysHphosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent
Reads supporting (aligned to +/- strand):  ref base G (57/48);  new base T (32/25);  total (89/73)
Fisher's exact test for biased strand distribution p-value = 8.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

CAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGTTTGC  >  CP009273/2880931‑2881225
                                                                                                                                                   |                                                                                                                                                   
cAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTCTCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAgg                                                                                                                                                     >  2:284892/1‑148 (MQ=255)
   tatTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCa                                                                                                                                                  <  2:606148/148‑1 (MQ=255)
    atTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCAt                                                                                                                                                 <  1:1106850/148‑1 (MQ=255)
     tttACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATg                                                                                                                                                <  2:1096950/148‑1 (MQ=255)
       tACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGtt                                                                                                                                              >  1:1080351/1‑148 (MQ=255)
       tACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGtt                                                                                                                                              <  2:1659507/148‑1 (MQ=255)
         cccTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGtttt                                                                                                                                            <  1:3082353/148‑1 (MQ=255)
         cccTTCGTGTAACCCACATTCCCTTTTAAAGCCAAAGAAACGCGTTACATCTTCCGACATGCCGGGTAACCATTTACGGGTTGTATGGGTATCGCCCACCCATAAAAATCCTTCATCCCATAATGGGTGATATTTCAGGCCaatttat                                                                                                                                            >  1:882263/1‑142 (MQ=255)
         cccTTCGGGAAACCCACAGTCCCTTTTTAAGCCAAAGAAACGCGTTTCGTCTTCCGCCATGCCGGGTTCCCATGTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGtttt                                                                                                                                            <  1:3082344/148‑1 (MQ=255)
            ttCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGc                                                                                                                                         <  1:3291863/148‑1 (MQ=255)
            ttCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGc                                                                                                                                         >  2:23512/1‑148 (MQ=255)
                tgtAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGt                                                                                                                                     >  1:2575016/1‑148 (MQ=255)
                tgtAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTACTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGt                                                                                                                                     <  2:3205324/148‑1 (MQ=255)
                tgtAACCCACATTCCCTTTTTAAGCCAAAGAAACGCGTGTCTTCTTCCGCCATGCCGGGTTCCCTTTTACGGGTTGTATGGGTATAGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGAAgtt                                                                                                                                     <  2:3208889/148‑3 (MQ=255)
                    aCCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTg                                                                                                                                 <  1:2549827/148‑1 (MQ=255)
                    aCCCACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTg                                                                                                                                 <  1:580722/148‑1 (MQ=255)
                     cccACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGa                                                                                                                                <  1:3022186/148‑1 (MQ=255)
                     cccACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGa                                                                                                                                >  2:2833487/1‑148 (MQ=255)
                     cccACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGa                                                                                                                                <  1:1490614/148‑1 (MQ=255)
                     cccACATTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGa                                                                                                                                <  2:199601/148‑1 (MQ=255)
                         caTTCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATaaa                                                                                                                            >  1:3229364/1‑148 (MQ=255)
                            tCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAg                                                                                                                         >  2:1517503/1‑148 (MQ=255)
                            tCCCTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAg                                                                                                                         <  1:861936/148‑1 (MQ=255)
                               cTTTTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTc                                                                                                                      <  1:1120770/148‑1 (MQ=255)
                                 ttttAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCgg                                                                                                                    >  1:801184/1‑148 (MQ=255)
                                     aaGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTAt                                                                                                                <  2:1836319/148‑1 (MQ=255)
                                       gCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATcc                                                                                                              >  2:913212/1‑148 (MQ=255)
                                          aaaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAg                                                                                                           <  1:782803/148‑1 (MQ=255)
                                          aaaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAg                                                                                                           <  1:782785/148‑1 (MQ=255)
                                            aGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTc                                                                                                         >  2:331078/1‑148 (MQ=255)
                                            aGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTc                                                                                                         >  1:1219894/1‑148 (MQ=255)
                                             gAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCg                                                                                                        >  1:718367/1‑148 (MQ=255)
                                                aCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATa                                                                                                     >  1:1808441/1‑148 (MQ=255)
                                                aCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATa                                                                                                     >  2:251208/1‑148 (MQ=255)
                                                aCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATa                                                                                                     >  1:1196796/1‑148 (MQ=255)
                                                aCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATa                                                                                                     >  1:1196121/1‑148 (MQ=255)
                                                 cgcgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATaa                                                                                                    <  1:2395167/148‑1 (MQ=255)
                                                   cgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGTGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATc                                                                                                  <  1:117271/148‑1 (MQ=255)
                                                   cgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATc                                                                                                  >  1:1427018/1‑148 (MQ=255)
                                                   cgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATc                                                                                                  >  1:2005258/1‑148 (MQ=255)
                                                   cgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATc                                                                                                  >  1:337982/1‑148 (MQ=255)
                                                   cgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATc                                                                                                  >  1:1478892/1‑148 (MQ=255)
                                                      ttcttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgc                                                                                               >  1:1887001/1‑148 (MQ=255)
                                                      ttcttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgc                                                                                               >  2:1286916/1‑148 (MQ=255)
                                                       tcttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgca                                                                                              <  2:1152505/148‑1 (MQ=255)
                                                       tcttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgca                                                                                              <  2:1152425/148‑1 (MQ=255)
                                                        cttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgcag                                                                                             >  1:588198/1‑148 (MQ=255)
                                                         ttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCGATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgcagc                                                                                            >  2:625206/1‑148 (MQ=255)
                                                         ttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgcagc                                                                                            >  2:2822466/1‑148 (MQ=255)
                                                          tcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgcagca                                                                                           <  1:551055/148‑1 (MQ=255)
                                                            ttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACt                                                                                         >  2:487513/1‑148 (MQ=255)
                                                             tcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACtt                                                                                        >  1:1345391/1‑148 (MQ=255)
                                                                gCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTaa                                                                                     <  1:3264926/148‑1 (MQ=255)
                                                                gCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTaa                                                                                     <  2:564686/148‑1 (MQ=255)
                                                                gCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTaa                                                                                     <  1:2623575/148‑1 (MQ=255)
                                                                 ccATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTaaa                                                                                    >  2:3268127/1‑148 (MQ=255)
                                                                  cATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAAt                                                                                   >  1:764066/1‑148 (MQ=255)
                                                                     gCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACg                                                                                <  2:2191868/148‑1 (MQ=255)
                                                                     gCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACg                                                                                >  2:1956930/1‑148 (MQ=255)
                                                                      ccGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGc                                                                               <  2:2956683/148‑1 (MQ=255)
                                                                         ggTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCAc                                                                            >  2:2276564/1‑148 (MQ=255)
                                                                          gTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACg                                                                           <  1:1819401/148‑1 (MQ=255)
                                                                          gTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACg                                                                           <  1:1289372/148‑1 (MQ=255)
                                                                          gTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACg                                                                           >  1:2501392/1‑148 (MQ=255)
                                                                           ttCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGc                                                                          >  2:2397496/1‑148 (MQ=255)
                                                                                aTTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAAt                                                                     <  1:2091493/148‑1 (MQ=255)
                                                                                  ttACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTg                                                                   >  1:378586/1‑148 (MQ=255)
                                                                                  ttACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTg                                                                   >  1:2376705/1‑148 (MQ=255)
                                                                                  ttACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTg                                                                   <  1:1200831/148‑1 (MQ=255)
                                                                                    aCGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGcc                                                                 >  2:2449552/1‑148 (MQ=255)
                                                                                        gTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCa                                                             <  1:1317990/148‑1 (MQ=255)
                                                                                               gggTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTaa                                                      >  1:904202/1‑148 (MQ=255)
                                                                                                  tATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAAtt                                                   >  1:704289/1‑148 (MQ=255)
                                                                                                   aTCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTg                                                  <  1:2519366/148‑1 (MQ=255)
                                                                                                    tCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTgg                                                 <  2:1802932/148‑1 (MQ=255)
                                                                                                    tCGCCAACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCTTGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCTATAATCGGCTGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTgg                                                 <  2:882263/148‑1 (MQ=255)
                                                                                                         cACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACTCCACGCTGAATTGCCAGACCCGGTCAATTGGCACGa                                            >  2:1980524/1‑148 (MQ=255)
                                                                                                          aCCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACATTAAATACGCCACGCTGAATAGCCAGCAGCGGTAAAAAGTCACGAc                                           >  2:2200667/1‑148 (MQ=255)
                                                                                                               tAAATATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCg                                      <  1:755005/148‑1 (MQ=255)
                                                                                                                aaaTATCCTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCgg                                     >  2:2817949/1‑148 (MQ=255)
                                                                                                                 aaTATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGa                                    >  2:2121346/1‑148 (MQ=255)
                                                                                                                 aaTATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGa                                    >  2:456769/1‑148 (MQ=255)
                                                                                                                  atatCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGAt                                   <  2:1407771/148‑1 (MQ=255)
                                                                                                                      ccTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGtt                               >  2:472801/1‑148 (MQ=255)
                                                                                                                       cTTCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTc                              <  2:1903767/148‑1 (MQ=255)
                                                                                                                        ttCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCa                             >  1:2423129/1‑148 (MQ=255)
                                                                                                                         tCATCCCATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCAc                            >  1:354779/1‑148 (MQ=255)
                                                                                                                             cccATAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGgcg                        >  2:256895/1‑148 (MQ=255)
                                                                                                                                aTAATGGGTGATATTTCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAg                     >  2:1726014/1‑148 (MQ=255)
                                                                                                                                         gATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCaaa            >  2:1823563/1‑148 (MQ=255)
                                                                                                                                          atatTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAAc           <  1:284892/148‑1 (MQ=255)
                                                                                                                                          atatTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAAc           <  1:284872/148‑1 (MQ=255)
                                                                                                                                            atTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCa         >  1:1441814/1‑148 (MQ=255)
                                                                                                                                              ttCAGTCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAgg       >  1:12972/1‑148 (MQ=255)
                                                                                                                                              ttCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAgg       >  2:330736/1‑148 (MQ=255)
                                                                                                                                                cAGTCCATGTTTTTGCAGGGACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGtt     <  2:1741088/148‑1 (MQ=255)
                                                                                                                                                  ggCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGTTTg   >  1:730573/1‑148 (MQ=255)
                                                                                                                                                   gCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGTTTgc  >  1:1054559/1‑148 (MQ=255)
                                                                                                                                                   gCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGTTTgc  >  1:1700393/1‑148 (MQ=255)
                                                                                                                                                   gCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATACGCCACGCTGAATTGCCAGCACCGGTAAATTGGCACGACTGCCGGATTGTTCACGGCGCAGGCCAGCAAACCAGGTTTgc  >  1:1054524/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: