Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,443,135 C→G 6.5% intergenic (‑66/+110) ghrB ← / ← PP_1262 2‑ketoaldonate reductase/hydroxypyruvate/glyoxylate reductase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,443,1350CG 6.5% 338.1 / 2.2 117intergenic (‑66/+110)ghrB/PP_12622‑ketoaldonate reductase/hydroxypyruvate/glyoxylate reductase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (53/49);  new base G (2/7);  total (60/57)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

GGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGTTTTCAGGTCTGC  >  NC_002947/1442986‑1443283
                                                                                                                                                     |                                                                                                                                                    
gggCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTAc                                                                                                                                                      <  1:4204897/150‑1 (MQ=255)
agattcttgccgcgcatgaggcagacctgaaaacccatcaggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTAc                                                                                                                                                      >  2:2275725‑M1/126‑150 (MQ=255)
  gCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACgg                                                                                                                                                    <  2:4408579/150‑1 (MQ=255)
  gCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACgg                                                                                                                                                    <  1:2210681/150‑1 (MQ=255)
  gCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGCCTGTTCCGGCACCTTCGCGGGCAAGCCCGCTCCTAcaa                                                                                                                                                    <  1:2491155/150‑3 (MQ=39)
   cAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACggg                                                                                                                                                   <  1:2044671/150‑1 (MQ=255)
    aGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACgggg                                                                                                                                                  >  1:1615727/1‑150 (MQ=255)
    aGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACgggg                                                                                                                                                  <  1:5407054/150‑1 (MQ=255)
        tCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCCACCCCCCACCAAACGGGGGgc                                                                                                                                              >  1:254107/1‑150 (MQ=255)
        tCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGgc                                                                                                                                              <  2:1615727/150‑1 (MQ=255)
            gTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGGAAGCCCCCGCCGACCGGGGGGCGCgg                                                                                                                                          >  1:4645530/1‑150 (MQ=255)
            gTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCCAGCCCCCACCAAACGGGGGGGGCgg                                                                                                                                          >  1:2522432/1‑150 (MQ=255)
            gTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCgg                                                                                                                                          >  2:3323515/1‑150 (MQ=255)
            gTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTATTTTTCATCGGGCAACTTCTCTTCAAGGCGGGGTGGGGCTTTGCTGGTCCCTTCGCGGGCAATCCCCCACCAAAGGGGGGGCGCgg                                                                                                                                          >  1:3617309/1‑150 (MQ=255)
            gTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACCGCTGAAGGCCCGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCCAGCCCCCGCCCACCGGGGGGCGCgg                                                                                                                                          >  1:3675455/1‑150 (MQ=255)
             ttGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGt                                                                                                                                         <  1:2227096/150‑1 (MQ=255)
              tGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTg                                                                                                                                        <  1:2878084/150‑1 (MQ=255)
                   tCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCACCTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGGGACCCCCCCCCGACGGGGGGGCGCggggggggg                                                                                                                                   >  1:3402995/1‑143 (MQ=255)
                   tCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCGTg                                                                                                                                   >  2:3826011/1‑150 (MQ=255)
                    gcagacctgaaaacccatcaggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGt                                                                                                                                  >  1:2426736‑M1/106‑150 (MQ=255)
                    cTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGt                                                                                                                                  <  2:2505029/150‑1 (MQ=255)
                    cTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGt                                                                                                                                  <  2:1406649/150‑1 (MQ=255)
                    cTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGt                                                                                                                                  <  2:5200423/150‑1 (MQ=255)
                    cTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGt                                                                                                                                  <  2:5200427/150‑1 (MQ=255)
                    cTTGCTGCAGGCGTTCGGCCATGGCCGAGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGt                                                                                                                                  <  2:1228935/150‑1 (MQ=255)
                     ttGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTa                                                                                                                                 <  1:825275/150‑1 (MQ=255)
                     ttGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTa                                                                                                                                 <  1:3261866/150‑1 (MQ=255)
                     ttGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCGGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTa                                                                                                                                 <  1:1062040/150‑1 (MQ=255)
                     ttGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCGGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTa                                                                                                                                 <  1:81633/150‑1 (MQ=255)
                     gtgctgcAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTa                                                                                                                                 <  2:3787233/149‑1 (MQ=255)
                     gtgctgcAGGCGTTCGGCCATGGCCGGGGGGATCCGGCTGAAGGCCAGGACGGTCTTTTTCCTCGGGCAACGGCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTa                                                                                                                                 <  1:4122117/149‑1 (MQ=255)
                        acctgaaaacccatcaggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGa                                                                                                                              >  1:9293‑M1/102‑150 (MQ=255)
                          gcAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGc                                                                                                                            <  1:273402/150‑1 (MQ=255)
                          gcAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGc                                                                                                                            <  1:273681/150‑1 (MQ=255)
                          ctgaaaacccatcaggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGc                                                                                                                            >  2:3377519‑M1/100‑150 (MQ=255)
                           tgaaaacccatcaggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCg                                                                                                                           >  1:483016‑M1/99‑150 (MQ=255)
                            gaaaacccatcaggtgctgatcgccaactaaacagcttctgcaatgcgggggctgacaatggtatgttgctaagctttgatttcgcctgttcaggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTTTAGGAGCgg                                                                                                                          >  1:4109364‑M1/98‑150 (MQ=255)
                            gaaaacccatcaggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCgg                                                                                                                          >  1:1848328‑M1/98‑150 (MQ=255)
                            gaaaacccatcaggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCgg                                                                                                                          >  2:117001‑M1/98‑150 (MQ=255)
                             aaaacccatcaggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGTGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCggg                                                                                                                         >  1:474661‑M1/97‑150 (MQ=255)
                             aaaacccatcaggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCggg                                                                                                                         >  1:654282‑M1/97‑150 (MQ=255)
                              gCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGc                                                                                                                        >  2:1738544/1‑150 (MQ=255)
                                gTTCGGCCATGGCCGGGGTGATACGGCTGAAGTCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACAAGGCACCGCGCCCCCCCGTAGGAGCGGGCtt                                                                                                                      <  2:4136211/150‑1 (MQ=25)
                                  tCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGc                                                                                                                    <  2:1995631/150‑1 (MQ=255)
                                   cGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGcc                                                                                                                   >  2:3504585/1‑150 (MQ=255)
                                    ggCCATGGCCGGGGTGCTCCGGCGGACGGCCAGGACGGTCTTTTTCCGCGGGCAACGTCTCTGCACGGCGGGGTGGGGCTGGGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGccc                                                                                                                  <  2:3675455/150‑1 (MQ=255)
                                    ggCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCGTGTAGGAGCGGGCTTGccc                                                                                                                  >  2:557097/1‑150 (MQ=255)
                                        ggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCgg                                                                                                                          >  1:5365396‑M1/86‑138 (MQ=255)
                                        ggtgctgatcgccaactaaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCgg                                                                                                                          <  2:5365396‑M1/53‑1 (MQ=255)
                                           gCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGt                                                                                                                                         <  1:2783686/120‑1 (MQ=255)
                                           gCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGt                                                                                                                                         >  2:2783686/1‑120 (MQ=255)
                                             cGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTg                                                                                                         >  2:2522157/1‑150 (MQ=255)
                                                    aTACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGACGGTGCCGGCAc                                                                                                  >  2:5090069/1‑150 (MQ=255)
                                                     tACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACa                                                                                                 >  2:3859575/1‑150 (MQ=255)
                                                      aCGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCGTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAAccg                                                                                                >  2:1599051/1‑148 (MQ=255)
                                                       cGGCTGCCGGCCTGGCGGGTCTTTTTCAGCGGGCACGTTCTGTGCTGGGGGGGGTGGGGCTGTGGGGGGCCCGGCGCGGGCACGCCCGCGCCTCGGGGGGGGGGGGGTGCCTTGTCGGCGCGGGCTTGCCCGCGAAGGTGCCGGAACAgg                                                                                               <  2:5260671/150‑1 (MQ=255)
                                                        taaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCggg                                                                                                                         >  2:2446241‑M1/70‑123 (MQ=255)
                                                        taaaaagctgctgcaatgcgggggcggacaatggtaggttgctaaccattgatttcgcctgttccggcaCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCggg                                                                                                                         <  1:2446241‑M1/54‑1 (MQ=255)
                                                         gCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCCCACCCCCGCCAACCGGGGGGGGGGGTGGCGTGGTGGGGCGGGGGTGGCCGCGCAGGGGGCCGGACAGGCg                                                                                             >  1:1845/1‑150 (MQ=255)
                                                           tGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAACCCCCCACCAACGGGGGGGCCCGGTGCCGTGGAGGAGCCGGCGTGCCCCCGCAGGTGCCGGCACAGGCGaa                                                                                           >  1:3910879/1‑150 (MQ=255)
                                                           tGAAGCCCTGGAGGGTCTTTTTCATCGGGAAACTTCTCTGAACGGGGGGGTGGGGCTGTGCTGGCCCCTTCGCGGGAAAGCCCGCTCCTAGGGGGGGGCGCGGTGCCTTTTAGTAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGaa                                                                                           <  2:3286378/150‑1 (MQ=255)
                                                                 cccggcgggTCTTTTACAGCGGCCAAGTTCGCTGCACGGCGGGGGGGGGCTGGGCTGGTCCGGGCGCGGGCACGCCCGCTCCTAGGGGGGGGCGCGGTGCCTGGTAGGAGGGGGCTGCCCCGCGAAGGTGCCGGAACAGGCGAAATCAAt                                                                                     <  2:4347443/143‑1 (MQ=255)
                                                                  ccggccGGTTTTTTTCTGGGGGCACCTTCGCGGCCCGGGGGGGGGGGGCGGTGCGGGGCCCTTCGGGGTCACCCCCGCTCCCCGGGGGGGGGGCGGTGCCTTGGCGGCGCGGGGTGGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATg                                                                                    <  2:727522/145‑1 (MQ=255)
                                                                  cAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCGTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATg                                                                                    >  2:1621715/1‑150 (MQ=255)
                                                                      cgggTCTTTTGCTGCGGGCACGTGCTGGGCACGGGGGGGGGGGGGTGGGCTGGCCCGGGCGCGGGCACGCCGGCGCGTCGGGGGGGGGGGGGGGCGTTGTGGGCGGGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTa                                                                                <  2:3402995/148‑1 (MQ=255)
                                                                        ggTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCGCCGACGGGGGGGCGCGGTGGCGGGGAGGGGCGGGCTTGCCCGCGAAGGTGGCGGGACAGGCGCAATCAATGGTTAGc                                                                              >  1:5252275/1‑150 (MQ=255)
                                                                         gacTTTTTCATCGGGCAACTTCTCTTCAATGAGGGGAGGGGATGTGATGGTCCCTTATCGGGCCAACCCGCTCCAAAGGGGGGGCGCGGTGAGGTGTAGGACGGGGCTTTCAAGCCTAGGGGCCCGACCAATCGAAATAAATAGTTAGCa                                                                             >  2:1836806/3‑150 (MQ=255)
                                                                            tttttCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCGCCTACCGGGGGGCGCGGTGCCTGGTTGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAAcc                                                                          >  2:1365883/1‑150 (MQ=255)
                                                                             ttttCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCGTGTAGGAGCGGGCTTGCCCGCGAAGGTTCCGGAACAGGCGAAATCAATGGTTAGAAACCt                                                                         >  2:1630808/1‑150 (MQ=255)
                                                                             ttttCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAACCCCGCTCCTACGGGGGGGGGCGGTGGCGTGTAGGGGCGGGCTTGGCCGCGAAGGTGCCGGGACACGCGGAATAAATGGTTAGCAAcca                                                                         >  1:4378500/1‑149 (MQ=255)
                                                                                  aTCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGGGCCGGGTCGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACCGGCGCAATCAATGGTTAGCAACCTACCag                                                                    >  2:2567462/1‑149 (MQ=255)
                                                                                      cgCACCTTCTCTGCCTGGGGGGGTGGGGGTGTGGTGGTCCCTTGGCGGGCAAGCCCCCGCCGCGGGGGGGGCGGGGTGCCTTGTCGGTGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTc                                                                <  2:2378256/149‑1 (MQ=255)
                                                                                       gaaccTTCGCTGCTCGGGGGGGTGGGGCTGTGGTGGTCCTTTCGCGGGCACGCCCGCCCCTCGGGGGGGGCGCGGTGCCTTTTAGGCGGGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTcc                                                               <  2:4299832/146‑1 (MQ=255)
                                                                                            ttCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCCCTCCTACGGGGGGGCGCGGTGCCGTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATAAATGGTTAGCAACCTACCATTGTCCGcccc                                                          >  1:5436669/1‑150 (MQ=255)
                                                                                             tctcTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGGAAGCCCGCTGCTACGGGGGGGCGGGGTGGCTTGTTGGTGGGGGGGTGGCCGCGGAGGTGGCGGTAAAGGCGAAATCAATTGTTAGGAACCTAGCATTGTGCGgcccc                                                         >  2:2290448/1‑150 (MQ=255)
                                                                                             tctcTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCAACCATTGTCCGccccc                                                         >  1:3548247/1‑150 (MQ=255)
                                                                                             tctcTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGAAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGccccc                                                         >  1:806869/1‑150 (MQ=255)
                                                                                                  gCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCCAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGGAGGAGCGGGCTTGCCCGCGCAGGGGCCGGAACAGGCGAAATCAATGGTTAGCAACCGACCCTTGTCCGCCCCCGCAtt                                                    >  2:3329343/1‑150 (MQ=255)
                                                                                                    aaGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGGAAACCCCCTCCAACGGGGGGGCGCGGTGGCGGGTAGGTGCGGGCTTGGCCGCGGAGGTGCCGGGACAGGCGGAATAAAAGGTTTTCAACCAACCAATTTGCGccccccccatgc                                                  >  1:5195396/1‑143 (MQ=255)
                                                                                                                   tgtgCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCGTGTGGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCCATGGTTAGCCACCTACCCTTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGtt                                   >  2:3353817/1‑150 (MQ=39)
                                                                                                                             taggaGCGGGCTTGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGc                         <  1:3152482‑M1/145‑1 (MQ=255)
                                                                                                                             CCTTCGCGGGCAAGCCCGCTCCAAAGGGGGGGCGCGGTGCCGTGTAGGAGCGGGCTTGCCCCCGAAGGTGCCGGAACACGCGAAATCAATGGTTAGCAACCAACCATTGTCCCCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGc                         <  2:3835466‑M1/150‑1 (MQ=255)
                                                                                                                              CTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCa                        <  2:3575800‑M1/150‑1 (MQ=255)
                                                                                                                                      GCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATgg                <  1:162111‑M1/150‑1 (MQ=255)
                                                                                                                                       cacgcccgcccCGCCGGGGGGGCGCGGTGCCTGGTTGGAGGGGGCTGGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATggg               <  2:4645530/147‑1 (MQ=255)
                                                                                                                                       CAAGCCCGCTCCTACAAGGCACCGCGCCCCCCCGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATggg               <  2:5178978‑M1/150‑1 (MQ=255)
                                                                                                                                        cccCCCGCTCCTCGGGGGGGGCGCGGGGCCTTTTCGCCGGGGGCTCCCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGt              >  1:3323515‑M1/4‑150 (MQ=255)
                                                                                                                                           CCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATg                 <  1:2833811‑M1/144‑1 (MQ=255)
                                                                                                                                           CCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGtttt           <  2:316940‑M1/150‑1 (MQ=255)
                                                                                                                                           CCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGtttt           <  1:5061261‑M1/150‑1 (MQ=255)
                                                                                                                                            CCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGTTTTc          <  1:1612105‑M1/150‑1 (MQ=255)
                                                                                                                                             CGCTCCTCCGGGGGGGCGCGGTGCCTTGTAGGAGGGGGGTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGTTTTCa         >  2:3291906‑M1/1‑150 (MQ=255)
                                                                                                                                             CGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGTTTTCa         <  1:3765252‑M1/150‑1 (MQ=255)
                                                                                                                                               CTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGTTTTCAgg       <  1:4285959‑M1/150‑1 (MQ=255)
                                                                                                                                                   TACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTa                                      <  1:1137060‑M1/115‑1 (MQ=255)
                                                                                                                                                   TACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTa                                      >  2:1137060‑M1/1‑115 (MQ=255)
                                                                                                                                                   TACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGTTTTCAGGTCTg   <  2:2543441‑M1/150‑1 (MQ=255)
                                                                                                                                                   TACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGTTTTCAGGTCTg   <  2:870307‑M1/150‑1 (MQ=255)
                                                                                                                                                    cCGGGGGGGCGCGGGGCCTTGTCGGAGCGGGCTGGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGTTTTCAGGTCTGc  >  2:5195396‑M1/2‑150 (MQ=255)
                                                                                                                                                     CGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTa                                      >  1:4036939‑M1/1‑113 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGCAGGATCACGTTGAAGTCTTGCTGCAGGCGTTCGGCCATGGCCGGGGTGATACGGCTGAAGGCCAGGACGGTCTTTTTCATCGGGCAACTTCTCTGCAAGGCGGGGTGGGGCTGTGCTGGTCCCTTCGCGGGCAAGCCCGCTCCTACGGGGGGGCGCGGTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGGTGCCGGAACAGGCGAAATCAATGGTTAGCAACCTACCATTGTCCGCCCCCGCATTGCAGCAGCTTTTTAGTTGGCGATCAGCACCTGATGGGTTTTCAGGTCTGC  >  NC_002947/1442986‑1443283

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