Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,035,421 G→T 91.1% intergenic (‑62/‑8) mobA ← / → yihD molybdopterin‑guanine dinucleotide synthase/DUF1040 protein YihD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,035,4210GT91.1% 535.8 / 42.6 199intergenic (‑62/‑8)mobA/yihDmolybdopterin‑guanine dinucleotide synthase/DUF1040 protein YihD
Reads supporting (aligned to +/- strand):  ref base G (8/10);  new base T (78/103);  total (86/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

AAGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGACGG  >  CP009273/4035276‑4035567
                                                                                                                                                 |                                                                                                                                                  
taGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCa                                                                                                                                                  <  2:1228293/147‑1 (MQ=255)
     aTCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACt                                                                                                                                             <  2:3042370/148‑1 (MQ=255)
          ttATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGaa                                                                                                                                        >  1:184246/1‑148 (MQ=255)
           tATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGaaa                                                                                                                                       <  1:2950879/148‑1 (MQ=255)
           tATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGaaa                                                                                                                                       <  2:2350929/148‑1 (MQ=255)
            aTCTACGCCGCCCATTCGTCTGGCTTTACCGCATGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAAt                                                                                                                                      <  1:3241865/148‑1 (MQ=255)
                aCGGCGGCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTaa                                                                                                                                  <  2:2218461/148‑1 (MQ=255)
                 cgccgcCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTaaa                                                                                                                                 <  1:561052/148‑1 (MQ=255)
                 cgccgcCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGAGTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTaaa                                                                                                                                 <  2:3203049/148‑1 (MQ=255)
                   ccgccCATTCGTCTGGCTTTACCGCCTGACAGCACAACTCCTGATATCGTCGTCATCAGATTCAACGATATCGCCTATTTTAATGTGGGATTTACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATTAAATGTAAACg                                                                                                                               >  1:2112016/1‑148 (MQ=255)
                    cgccCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGt                                                                                                                              <  1:415526/148‑1 (MQ=255)
                         aTTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCtgaa                                                                                                                         >  1:806160/1‑148 (MQ=255)
                           tCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaatg                                                                                                                       <  2:1694409/148‑1 (MQ=255)
                           tCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaatg                                                                                                                       >  1:1850367/1‑148 (MQ=255)
                           tCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaatg                                                                                                                       <  2:2390640/148‑1 (MQ=255)
                           tCGTCTGGCTTTACCGCCTGCCAGCACAACGACTTTTATCGTCGTCATAAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTTCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATTAAATGTAAACGTCtgaatg                                                                                                                       <  1:2830424/148‑1 (MQ=255)
                            cGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaatga                                                                                                                      <  1:2411342/148‑1 (MQ=255)
                                 ggCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTa                                                                                                                 <  2:1980686/148‑1 (MQ=255)
                                  gCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTAt                                                                                                                <  2:2277699/148‑1 (MQ=255)
                                  gCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTAt                                                                                                                <  1:2753130/148‑1 (MQ=255)
                                  gCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGGTAt                                                                                                                >  2:1247346/1‑148 (MQ=255)
                                    tttACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTg                                                                                                              >  2:1009927/1‑148 (MQ=255)
                                       aCCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAAc                                                                                                           <  2:1911445/148‑1 (MQ=255)
                                         cGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActc                                                                                                         <  2:1208410/148‑1 (MQ=255)
                                         cGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTAAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActc                                                                                                         <  1:2723079/148‑1 (MQ=255)
                                          gCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActcc                                                                                                        <  1:2075528/148‑1 (MQ=255)
                                          gCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActcc                                                                                                        >  1:2699980/1‑148 (MQ=255)
                                          gCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActcc                                                                                                        <  1:2345704/148‑1 (MQ=255)
                                           ccTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActcct                                                                                                       <  2:505342/148‑1 (MQ=255)
                                                cAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTccagc                                                                                                  <  1:772040/148‑1 (MQ=255)
                                                cAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTccagc                                                                                                  <  1:1329740/148‑1 (MQ=255)
                                                  gCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTccagcca                                                                                                >  1:676588/1‑148 (MQ=255)
                                                   cacaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTccagccag                                                                                               >  1:232311/1‑148 (MQ=255)
                                                    acaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCcagccagc                                                                                              <  1:1287467/148‑1 (MQ=255)
                                                    acaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCcagccagc                                                                                              <  1:2936788/148‑1 (MQ=255)
                                                    acaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCcagccagc                                                                                              <  2:3294356/148‑1 (MQ=255)
                                                    acaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCcagccagc                                                                                              <  1:3128689/148‑1 (MQ=255)
                                                    acaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCcagccagc                                                                                              <  1:1475386/148‑1 (MQ=255)
                                                    acaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCcagccagc                                                                                              <  1:2841062/148‑1 (MQ=255)
                                                    acaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCcagccagc                                                                                              <  2:3159083/148‑1 (MQ=255)
                                                     caACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCagccagcc                                                                                             <  1:920824/148‑1 (MQ=255)
                                                     caACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCagccagcc                                                                                             <  1:862039/148‑1 (MQ=255)
                                                      aaCGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCt                                                                                            >  2:1061263/1‑148 (MQ=255)
                                                       aCGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTg                                                                                           >  2:2208536/1‑148 (MQ=255)
                                                           cTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCaa                                                                                       <  1:1897853/148‑1 (MQ=255)
                                                             gTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCaaaa                                                                                     <  2:78483/148‑1 (MQ=255)
                                                               tATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAag                                                                                   <  1:452583/148‑1 (MQ=255)
                                                                aTCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAaga                                                                                  >  1:2954639/1‑148 (MQ=255)
                                                                aTCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAaga                                                                                  >  1:2183588/1‑148 (MQ=255)
                                                                 tcgtcgTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAagag                                                                                 <  2:947134/148‑1 (MQ=255)
                                                                  cgtcgtCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGc                                                                                <  2:2329667/148‑1 (MQ=255)
                                                                  cgtcgtCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGc                                                                                <  2:2330039/148‑1 (MQ=255)
                                                                  cgtcgtCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGc                                                                                <  1:1446201/148‑1 (MQ=255)
                                                                        catcaGATTCACCGAGATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGAcc                                                                          <  2:1813368/148‑1 (MQ=255)
                                                                         atcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCt                                                                         <  1:604546/148‑1 (MQ=255)
                                                                          tcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCtt                                                                        <  2:661225/148‑1 (MQ=255)
                                                                           caGATTCACCGTTATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTATGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTa                                                                       <  1:1070621/148‑1 (MQ=255)
                                                                             gATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAc                                                                     <  1:984136/148‑1 (MQ=255)
                                                                             gATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAc                                                                     >  2:2660351/1‑148 (MQ=255)
                                                                              aTTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAcc                                                                    >  1:201113/1‑148 (MQ=255)
                                                                                   ccGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACctgctg                                                               >  2:2004281/1‑148 (MQ=255)
                                                                                        atCTCCTCTATAATTGTGGGATTGACCCTGCTAACTAGTAAGACTAAAGAGTAAGGATCATCACTATTAAATTTAAAACTCTTAATGAATAAATTTTACTCCTCAAGCCAGCCTGGCAAAAATAGCCAGACCTAAACCTTCTGAAAtt                                                          >  2:866259/1‑148 (MQ=255)
                                                                                        atCGCCTCTTTTATTGTGGGATTGACCCTGCGAACGTGTCTGTCTCAAGAGTAAGGATCATCACTCTGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCGGGCAAAAAGAGCCAGACCGTAACCTGCTGCAAtt                                                          >  2:2521839/1‑148 (MQ=255)
                                                                                         tCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAttt                                                         <  2:149563/148‑1 (MQ=255)
                                                                                         tCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAttt                                                         <  2:1133622/148‑1 (MQ=255)
                                                                                             ctctTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGc                                                     <  2:448648/148‑1 (MQ=255)
                                                                                               ctTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAg                                                   >  1:1898508/1‑148 (MQ=255)
                                                                                                  ttATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGaaa                                                >  1:1146168/1‑148 (MQ=255)
                                                                                                  ttATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGaaa                                                >  1:2097010/1‑148 (MQ=255)
                                                                                                     ttGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTg                                             <  1:1547555/148‑1 (MQ=255)
                                                                                                     ttGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTg                                             <  1:1649715/148‑1 (MQ=255)
                                                                                                     ttGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTg                                             <  1:355934/148‑1 (MQ=255)
                                                                                                      tgtgGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTgg                                            <  1:1043577/148‑1 (MQ=255)
                                                                                                      tgtgGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTgg                                            >  2:1963720/1‑148 (MQ=255)
                                                                                                      tgtgGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTgg                                            <  1:2356062/148‑1 (MQ=255)
                                                                                                          ggATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGaa                                        <  1:2129965/148‑1 (MQ=255)
                                                                                                           gATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGaaa                                       <  1:1489633/148‑1 (MQ=255)
                                                                                                             ttGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAaga                                     <  2:3234244/148‑1 (MQ=255)
                                                                                                                 cccTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCa                                 >  2:844368/1‑148 (MQ=255)
                                                                                                                  ccTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTGTGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAg                                >  1:3145711/1‑148 (MQ=255)
                                                                                                                    tGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGt                              >  2:2539990/1‑148 (MQ=255)
                                                                                                                      cTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGttt                            <  1:1973010/148‑1 (MQ=255)
                                                                                                                        aaCGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTg                          <  1:2244420/148‑1 (MQ=255)
                                                                                                                         aCGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGa                         <  1:1624088/148‑1 (MQ=255)
                                                                                                                             gtctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGc                     >  2:2075139/1‑148 (MQ=255)
                                                                                                                             gtctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGc                     >  2:1106168/1‑148 (MQ=255)
                                                                                                                               ctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGa                   >  2:2833123/1‑148 (MQ=255)
                                                                                                                                 gtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAAc                 <  2:107906/148‑1 (MQ=255)
                                                                                                                                  tctcAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACt                >  1:2462724/1‑148 (MQ=255)
                                                                                                                                   ctcAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTg               >  1:627478/1‑148 (MQ=255)
                                                                                                                                   ctcAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTg               >  2:2444396/1‑148 (MQ=255)
                                                                                                                                   ctcAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTg               >  2:1285625/1‑148 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa          >  2:2773506/1‑148 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa          >  1:678773/1‑148 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa          >  2:1787788/1‑148 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa          >  2:1788531/1‑148 (MQ=255)
                                                                                                                                         agTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGAt         <  2:1170807/148‑1 (MQ=255)
                                                                                                                                              ggATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGAc    <  2:2699980/148‑1 (MQ=255)
                                                                                                                                              ggATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGAc    >  1:792324/1‑148 (MQ=255)
                                                                                                                                                aGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGACgg  <  2:1898508/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AAGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGACGG  >  CP009273/4035276‑4035567

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: