Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,557,402 C→G 6.8% G271G (GGG→GGC PP_1368 ← glycosyltransferase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,557,4020CG 6.8% 452.0 / 16.2 144G271G (GGG→GGCPP_1368glycosyltransferase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (69/65);  new base G (6/4);  total (75/69)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 146 total aligned reads displayed.

GTCCCAGCAGGCACACAAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTCCCCAGTGGATGCAGGTC  >  NC_002947/1557260‑1557549
                                                                                                                                              |                                                                                                                                                   
gTCCCAGCAGGCACACAAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGgcgc                                                                                                                                              >  1:554606/1‑150 (MQ=255)
  cccAGCAGGCACACAAACACCGGCCGCGCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGcgccg                                                                                                                                            <  2:1161566/150‑1 (MQ=255)
     agcagGCACACAAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGcgccgaac                                                                                                                                         <  2:741365/150‑4 (MQ=255)
        agGCACACAAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGgc                                                                                                                                      >  1:1887494/1‑150 (MQ=255)
           cacacaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGgcagc                                                                                                                                   <  1:1364930/150‑1 (MQ=255)
            acacaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGgcagca                                                                                                                                  >  1:239993/1‑150 (MQ=255)
            acacaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGgcagca                                                                                                                                  >  1:1815998/1‑150 (MQ=255)
            acacaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGgcagca                                                                                                                                  >  1:2063857/1‑150 (MQ=255)
            acacaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGgcagca                                                                                                                                  >  2:546688/1‑150 (MQ=255)
             cacaAACACCGGCCCCCCAGCCAGCAGATAAAAGTCATACGAACACGCCCGCCATGCTAGGATTGTCGCATCTTGGGTGGAGCTGCCCACCGAGGGTTGCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGcagcag                                                                                                                                 >  2:2080634/1‑150 (MQ=255)
              acaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCGCCGGGgccgccgaa                                                                                                                                       >  2:2313529/1‑140 (MQ=255)
              acaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCGCCGGGgccgccgaa                                                                                                                                       <  1:2313529/143‑4 (MQ=255)
              acaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGcgccgc                                                                                                                                           <  1:977326/139‑1 (MQ=255)
              acaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGa                                                                                                                                >  1:232603/1‑150 (MQ=255)
              acaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGa                                                                                                                                >  1:438994/1‑150 (MQ=255)
              acaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGa                                                                                                                                >  1:119217/1‑150 (MQ=255)
              acaAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGa                                                                                                                                >  1:179027/1‑150 (MQ=255)
                aaaCACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCa                                                                                                                                         <  2:16294/139‑1 (MQ=255)
                aaaCACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACt                                                                                                                              >  2:661323/1‑150 (MQ=255)
                 aaCACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACAACCCGGCCTATGAGGATTGCCGCAGCTGCTGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTg                                                                                                                             >  1:1856103/1‑150 (MQ=255)
                  acacCGGCCCCCCAGCCAGTAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGa                                                                                                                            >  2:819443/1‑150 (MQ=255)
                  acacCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCGCCGGGgccgccgaa                                                                                                                                       >  1:1216872/1‑136 (MQ=255)
                  acacCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGa                                                                                                                            >  2:2377451/1‑150 (MQ=255)
                    acCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAAt                                                                                                                          >  2:1699984/1‑150 (MQ=255)
                     ccGGCCCCCCAGCCAGCAGATAAAAGTAATACGTCACCCTTCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAAtt                                                                                                                         <  1:1933127/150‑1 (MQ=255)
                        gCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGc                                                                                                                      <  2:232603/150‑1 (MQ=255)
                           ccccAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCgg                                                                                                                   >  2:175051/1‑150 (MQ=255)
                           ccccAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCgg                                                                                                                   >  2:1482624/1‑150 (MQ=255)
                            cccAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAAAAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGaa                                                                                                                           >  2:1306982/1‑141 (MQ=255)
                                   agcagATAAAAGTAATACGTACCCGCCCGCCATATTAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGc                                                                                                           <  1:11847/150‑1 (MQ=255)
                                   agcagATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGc                                                                                                           <  1:1482624/150‑1 (MQ=255)
                                   agcagATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGc                                                                                                           <  2:438994/150‑1 (MQ=255)
                                       gATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCCGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACa                                                                                                       <  2:2085661/150‑1 (MQ=255)
                                         tAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAGGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGgcagc                                                                                                                                   >  1:1262168/1‑120 (MQ=255)
                                         tAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAGGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGgcagc                                                                                                                                   <  2:1262168/120‑1 (MQ=255)
                                                 acaCGTAACCGCCCGCCATATCACGATTGCCGCAGCTGCGGTCGAGCTGGCCACCAAGGGTGTCAACAGGCTGGCCGGCGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCAc                                                                                             <  2:357731/148‑1 (MQ=255)
                                                 aTACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGTGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCAc                                                                                             <  2:1373999/150‑1 (MQ=255)
                                                 aTACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCAc                                                                                             <  2:2251299/150‑1 (MQ=255)
                                                     gTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGgcagca                                                                                                                                  >  1:2228629/1‑109 (MQ=255)
                                                     gTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCCAGGCCAGCTCGGCGGCGCGGCCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACAAgg                                                                                         >  2:1059498/1‑150 (MQ=255)
                                                        aCCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACATGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGa                                                                                                                                >  2:418739/1‑108 (MQ=255)
                                                          cgcccgccATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCt                                                                                    <  2:1905460/150‑1 (MQ=255)
                                                           gcccgccATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTg                                                                                   <  1:2377451/150‑1 (MQ=255)
                                                           gcccgccATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTg                                                                                   <  2:1815998/150‑1 (MQ=255)
                                                            cccgccATATTAGGATTGCCGCAGCCGCGGTCGAGCTGCCCACCGATGGTGTCAACAGGCTGGCCGACGTCCTTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTgg                                                                                  <  2:1145894/150‑1 (MQ=255)
                                                            cccgccATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTgg                                                                                  <  1:661323/150‑1 (MQ=255)
                                                               gccATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGaa                                                                               >  1:409295/1‑150 (MQ=255)
                                                                ccATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAAt                                                                              >  2:414245/1‑150 (MQ=255)
                                                                  atacCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTTGGAATcc                                                                            <  1:1299883/150‑1 (MQ=255)
                                                                     tCAGGATTGCCGCAGGTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGATGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGt                                                                         <  2:334709/150‑1 (MQ=255)
                                                                     tCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTgg                                                                                  >  1:1317258/1‑141 (MQ=255)
                                                                     tCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGt                                                                         >  2:1453398/1‑150 (MQ=255)
                                                                        ggATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCa                                                                      <  2:1887494/150‑1 (MQ=255)
                                                                         gATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAg                                                                     <  1:1729081/150‑1 (MQ=255)
                                                                         gATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAACCCAGTGCAg                                                                     <  2:1774894/150‑1 (MQ=255)
                                                                          aTTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCa                                                                           <  1:407930/143‑1 (MQ=255)
                                                                          aTTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCa                                                                           >  2:407930/1‑143 (MQ=255)
                                                                          aTTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAgg                                                                    >  1:1700527/1‑150 (MQ=255)
                                                                           ttGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGt                                                                   <  1:819443/150‑1 (MQ=255)
                                                                               cGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTg                                                                                   <  1:717909/130‑1 (MQ=255)
                                                                               cGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTg                                                                                   >  2:717909/1‑130 (MQ=255)
                                                                               cGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTAc                                                               >  1:1102054/1‑150 (MQ=255)
                                                                                 cAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTAccc                                                             <  2:1820742/150‑1 (MQ=255)
                                                                                      gCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAg                                                        <  2:8819/150‑1 (MQ=255)
                                                                                          tCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGcccc                                                    >  1:1237990/1‑150 (MQ=255)
                                                                                             aGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCgca                                                 <  2:658946/150‑1 (MQ=255)
                                                                                                   ccaccaAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAg                                                                          <  2:1370863/119‑1 (MQ=255)
                                                                                                   ccaccaAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAg                                                                          >  1:1370863/1‑119 (MQ=255)
                                                                                                     accaAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCt                                                                                                               >  2:555168/1‑80 (MQ=255)
                                                                                                     accaAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGc                                                                       <  1:1499940/120‑1 (MQ=255)
                                                                                                     accaAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAgg                                         <  2:1346365/150‑1 (MQ=255)
                                                                                                     accaAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAgg                                         >  1:196061/1‑150 (MQ=255)
                                                                                                      ccaAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTgg                                                                                  <  2:1959417/108‑1 (MQ=255)
                                                                                                      ccaAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTgg                                                                                  >  1:1959417/1‑108 (MQ=255)
                                                                                                       caAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTg                                                                                   <  2:314622/106‑1 (MQ=255)
                                                                                                       caAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTg                                                                                   >  1:314622/1‑106 (MQ=255)
                                                                                                        aaGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTAc                                      >  1:324280/1‑150 (MQ=255)
                                                                                                          gggTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACaa                                    <  1:2034459/150‑1 (MQ=255)
                                                                                                          gggTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACaa                                    <  1:649259/150‑1 (MQ=255)
                                                                                                                 aaCAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGa                                                                                                             >  2:1656025/1‑70 (MQ=255)
                                                                                                                 aaCAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGa                                                                                                             <  1:1656025/70‑1 (MQ=255)
                                                                                                                 aaCAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCa                                                                           <  1:1641331/104‑1 (MQ=255)
                                                                                                                      gCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAgcg                        >  1:2234192/1‑150 (MQ=255)
                                                                                                                        tGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCa                      <  2:2234192/150‑1 (MQ=255)
                                                                                                                        tGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCa                      >  1:133285/1‑150 (MQ=255)
                                                                                                                        tGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCa                      <  2:133285/150‑1 (MQ=255)
                                                                                                                        tGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCACGTCTACCCCCAAGACCCGCAGCACCAGGTACAAGACGCTGTAGCGCa                      >  2:679087/1‑150 (MQ=255)
                                                                                                                         ggCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTgg                                                                                  <  1:1406076/89‑1 (MQ=255)
                                                                                                                           ccGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTc                   <  1:799690/150‑1 (MQ=255)
                                                                                                                                gTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTCCCCAg              >  2:1554395/1‑150 (MQ=255)
                                                                                                                                  cAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTCCCCAGTg            <  2:2063857/150‑1 (MQ=255)
                                                                                                                                     ttCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTCCCCAGTGGAt         >  1:107454/1‑150 (MQ=255)
                                                                                                                                      ctggcggcGCCGGGGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTCCCCAGTGGATg        >  2:1595813/3‑150 (MQ=255)
                                                                                                                                       tggcggcGCCGGGGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAg                                          <  2:682140/114‑1 (MQ=255)
                                                                                                                                       tggcggcGCCGGGGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAg                                          >  1:682140/2‑115 (MQ=255)
                                                                                                                                       tggcggcGCCGGGGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTCCCCAGTGGATGc       >  2:1163312/2‑150 (MQ=255)
                                                                                                                                        ggcggcCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTCCCCAGTGGATGCa      >  1:1038481/1‑150 (MQ=255)
                                                                                                                                         gcggcCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTCCCCAGTGGATGCAg     >  1:404199/1‑150 (MQ=255)
                                                                                                                                            gcCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTCCCCAGTGGATGCAGGTc  <  2:2365716/150‑1 (MQ=255)
                                                                                                                                              cccGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTa                                       >  2:201704/1‑111 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GTCCCAGCAGGCACACAAACACCGGCCCCCCAGCCAGCAGATAAAAGTAATACGTAACCGCCCGCCATATCAGGATTGCCGCAGCTGCGGTCGAGCTGCCCACCAAGGGTGTCAACAGGCTGGCCGACGTCAGTTCGGCGGCCCCGGCGCCGCCAGGCAGCAGACTGAATTGGCCGGCGCTGAGCGACAGCATCTGCACCAGGAAGCTGGGAATCCAGTGCAGGTCTACCCCCAAGCCCCGCAGCACCAGGTACAAGACGCTGTAGCGCAGTCCCCAGTGGATGCAGGTC  >  NC_002947/1557260‑1557549

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: