Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,923 A→T 5.1% intergenic (+25/+84) PP_5256 → / ← amaA cyclic nucleotide‑binding protein/L‑pipecolate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,9230AT 5.1% 755.2 / 9.8 233intergenic (+25/+84)PP_5256/amaAcyclic nucleotide‑binding protein/L‑pipecolate oxidase
Reads supporting (aligned to +/- strand):  ref base A (107/113);  new base T (10/3);  total (117/116)
Fisher's exact test for biased strand distribution p-value = 8.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

GTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGCGAAAAACCGGACTCATAGAACGCTG  >  NC_002947/6005774‑6006071
                                                                                                                                                       |                                                                                                                                                     
gaggcggacagcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTa                                                                                                                                                       <  2:6873931‑M1/16‑1 (MQ=255)
gaggcggacagcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTa                                                                                                                                                       <  2:6525364‑M1/16‑1 (MQ=255)
 tGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTaa                                                                                                                                                      <  1:3881425/150‑1 (MQ=255)
  ggcggacagcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTaaa                                                                                                                                                     <  2:2092237‑M1/18‑1 (MQ=255)
  ggCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTaaa                                                                                                                                                     >  2:1397978/1‑150 (MQ=255)
   gcggacagcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAAc                                                                                                                                                    >  2:4165896‑M1/132‑150 (MQ=255)
   gCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAAc                                                                                                                                                    <  1:667053/150‑1 (MQ=255)
   gCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAAc                                                                                                                                                    >  2:8184095/1‑150 (MQ=255)
    cggacagcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAAcc                                                                                                                                                   >  1:7126750‑M1/131‑150 (MQ=255)
      aaGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCACATCGCGGGTAAAaccc                                                                                                                                                 >  1:4258936/1‑149 (MQ=255)
       acagcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGc                                                                                                                                                >  1:5384942‑M1/128‑150 (MQ=255)
       acagcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGc                                                                                                                                                >  1:3345323‑M1/128‑150 (MQ=255)
        gTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCTGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCt                                                                                                                                               <  1:3626156/150‑1 (MQ=255)
        gTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTATTAGCGGGTAAAccccca                                                                                                                                               >  2:4876652/1‑147 (MQ=255)
        cagcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCt                                                                                                                                               >  2:8759561‑M1/127‑150 (MQ=255)
        cagcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCt                                                                                                                                               >  1:1857711‑M1/127‑150 (MQ=255)
         agcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTc                                                                                                                                              >  2:7602990‑M1/126‑150 (MQ=255)
         agcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTc                                                                                                                                              >  1:2912297‑M1/126‑150 (MQ=255)
         agcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTc                                                                                                                                              >  1:1364327‑M1/126‑150 (MQ=255)
          gcctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                                                                             >  2:438434‑M1/125‑150 (MQ=255)
          gCAGACGGCGTTGAGCACATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGGGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                                                                             <  2:8455217/150‑4 (MQ=255)
          gCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                                                                             <  2:2683642/150‑1 (MQ=255)
          gCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                                                                             <  1:7305070/150‑1 (MQ=255)
          gCAGACGGCGTTGAGCAAAGTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                                                                             <  1:5949990/150‑4 (MQ=255)
            ctgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACACCctccaa                                                                                                                                           <  1:4876652‑M1/28‑3 (MQ=255)
            aGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAAaccccaccaa                                                                                                                                           >  2:575578/1‑143 (MQ=255)
            aGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                           <  1:7927760/150‑1 (MQ=255)
             tgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                          >  1:8521876‑M1/122‑150 (MQ=255)
             tgccgattccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                          >  2:7433513‑M1/122‑150 (MQ=255)
             gACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                          <  2:3205380/150‑1 (MQ=255)
               cGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                        <  1:819122/150‑1 (MQ=255)
               cGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                        <  1:5261939/150‑1 (MQ=255)
               cGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                        <  1:1315317/150‑1 (MQ=255)
               cGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                        <  2:587956/150‑1 (MQ=255)
                ggCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAag                                                                                                                                       <  2:4993873/150‑1 (MQ=255)
                   gTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAagtag                                                                                                                                    <  1:8244046/150‑1 (MQ=255)
                   gTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAagtag                                                                                                                                    <  1:2562248/150‑1 (MQ=255)
                    ttGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAagtagt                                                                                                                                   <  1:6492677/150‑1 (MQ=255)
                    tccgttttcgccgatgagcggtgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAAACCCCACCAacaagttca                                                                                                                                   <  1:1502541‑M1/36‑4 (MQ=255)
                     tGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTg                                                                                                                                  <  2:6284900/150‑1 (MQ=255)
                       aGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGcc                                                                                                                                >  2:7804738/1‑149 (MQ=255)
                         cAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCAcc                                                                                                                              <  2:1695640/150‑1 (MQ=255)
                             ttGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCt                                                                                                                          <  1:5583271/150‑1 (MQ=255)
                             ttGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCt                                                                                                                          <  2:6181794/150‑1 (MQ=255)
                             ttGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCt                                                                                                                          <  2:2787838/150‑1 (MQ=255)
                                 aGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCAgaga                                                                                                                      >  2:7109807/1‑150 (MQ=255)
                                  gCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGagag                                                                                                                     <  2:3072294/150‑1 (MQ=255)
                                  gCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTAGCGGGTTAAACCCCAACTACAACTAGTGGACCCCCGGggag                                                                                                                     >  2:1843596/1‑150 (MQ=255)
                                   cGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAgg                                                                                                                    <  2:2629799/150‑1 (MQ=255)
                                     gTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCAAATCGCGGGTAAAACCCCACCGACAAGAAGTGCACCCCCAGAGAGGGc                                                                                                                  >  2:5312038/1‑150 (MQ=255)
                                        ccGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCtgt                                                                                                               <  2:1678325/150‑1 (MQ=255)
                                        ccGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCtgt                                                                                                               <  2:6829376/150‑1 (MQ=255)
                                        ccGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCtgt                                                                                                               <  2:5282891/150‑1 (MQ=255)
                                         tgtgagcgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTAGTGCACCGCCTGAGAGGGCtgtg                                                                                                              >  2:6194626‑M1/94‑150 (MQ=255)
                                         cggcggCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCtgtg                                                                                                              <  1:2038130/150‑1 (MQ=255)
                                         cggcggCAGAGCCGGGAGTCGTGGGTGTTGCAGGAGCCGGCTTCAATCAAGAAGGGTGGCGTCCAGAGTTGGTTGCACAAGTGAGCCTGTGCCCGCCGCGTCGCGGGGAAACGCGCGCCTACAAGTAGTGCACCCCCGGGGAGGGCggtg                                                                                                              >  2:1502541/1‑150 (MQ=255)
                                          ggcggcAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTgtgg                                                                                                             <  2:1651409/150‑1 (MQ=255)
                                            cggcAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTgtggtg                                                                                                           <  1:1022842/150‑1 (MQ=255)
                                               cgcgccttcgctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAAc                                                                                                        >  2:888417‑M1/88‑150 (MQ=255)
                                               cAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAAc                                                                                                        <  2:1872264/150‑1 (MQ=255)
                                                  agCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTg                                                                                                     <  2:3641654/150‑1 (MQ=255)
                                                   gCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGt                                                                                                    <  2:7303259/150‑1 (MQ=255)
                                                     cGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAg                                                                                                  <  1:4478087/150‑1 (MQ=255)
                                                     cGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAg                                                                                                  <  2:3139268/150‑1 (MQ=255)
                                                       ggAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGa                                                                                                <  1:8797784/150‑1 (MQ=255)
                                                         ctgcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGc                                                                                              >  1:141631‑M1/78‑150 (MQ=255)
                                                           gcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCgg                                                                                           >  2:2508331‑M1/76‑150 (MQ=255)
                                                           gcggacagcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCgg                                                                                           >  1:2164762‑M1/76‑150 (MQ=255)
                                                            cGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCggg                                                                                          <  1:6990725/150‑1 (MQ=255)
                                                             gttgttgttGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCCGGGTAAACCCCACCAACAACTAGTGCACCCCCAGAGAGGGCGGTGGTGAACCAGGAGGAGCgggg                                                                                         >  1:6391992/1‑149 (MQ=255)
                                                                  gcgttctatgagtccggtttttcgctgtaccacgcggggcagtgcctgagggttgtgttgtaagtactGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTAcc                                                                                    >  2:5110671‑M1/69‑150 (MQ=255)
                                                                    ttgCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTTAAACCGCGCCTAAAAGTAGTGCACCGCCGGAGAGGGCTGTTGTGTACCGGTAGGGGCGGGttttaccg                                                                                  >  1:3530113/1‑145 (MQ=255)
                                                                        aGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCAGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGaa                                                                              >  2:2883630/1‑150 (MQ=255)
                                                                        aGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCGCCTACAAGTAGTGCACCGCCGGAGAGGGCTGTGGTGAACCGGTAGGCGCGGGTTTACCcgcgca                                                                              >  2:8405192/1‑148 (MQ=39)
                                                                        aGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGGTGGTTGCACAAGTGAGCATGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGTAGTGCCCCGCCGGAGAGGGCGGGGGGGAACCGGAAGGAGCGGGTTTACCCCCGaa                                                                              >  2:4256230/1‑150 (MQ=38)
                                                                          aaGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGGTAAACCCCGCCGACCACTAGTGCACCCCCAGAGAGGGCGGTGGTGAAACAGGGGGGGCGGGGTTTAccccccagg                                                                             >  2:11691/1‑144 (MQ=255)
                                                                           aGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGGAAACCCCCTCCGACAAGTAGTGCCCCCCCAGAGAGGGCGGTGGTGAACCGGTAGGAGCGGGTTTTCCCCCCAagag                                                                           >  1:7329322/1‑150 (MQ=255)
                                                                             ccGGCTTCAATCACGAAGGGTGGCTTCCTGTGTTGGTTGCAAAAGTGAGCCTGTGCCGGCCTCTTCGGGGGTAAACCCTCTCCTACAAGTAGTGCACCGCCTGAGTGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                         <  2:8002977/150‑1 (MQ=38)
                                                                               ggCTTCAATCAAGAAGGGTGGCTTCCTGCGTTGGTTGCACACGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGATTACCCGCGAAGAGGCCa                                                                       <  2:583969/150‑1 (MQ=255)
                                                                                      atcaaGACGGGTGGCTTCCTGTGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGTACTTa                                                                <  2:157498/150‑1 (MQ=255)
                                                                                         gctgtaccacgcggggcagtgcctgagggttgtgttgtaagtaccGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCgg                                                                                           <  1:1080899‑M1/75‑1 (MQ=255)
                                                                                         gctgtaccacgcggggcagtgcctgagggttgtgttgtaagtaccGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCgg                                                                                           >  2:1080899‑M1/46‑120 (MQ=255)
                                                                                         aagaagGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACAAGTAGTGCACCGCCTGGGAGGGCTGTGGTGAACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACATacaa                                                             >  2:3520210/1‑150 (MQ=255)
                                                                                           gaagGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCACCGACAAGTAGTGCACCCCCAGAGAGGGCGGTGGTGAACCGGTAGGCGCGGGTTTACCCGCGAAGAGGCCCCTACATacaaca                                                           >  1:5763234/1‑150 (MQ=37)
                                                                                           gaagGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCGACAAGTAGTGCACCGCCGGGGGGGGCTGTGGGGAACCGGTGGGAGCGGGGTTAACCCCGAAGAGGGCCCTACTTacaaca                                                           >  2:5467753/1‑150 (MQ=255)
                                                                                                                   gCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACAGCCCcgcgg                                   >  2:6570219/1‑149 (MQ=255)
                                                                                                                     acaAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCGACAAGTAGTGCACCCCCCGAGAGGGCGGTGGTGTACCGGGAGGGGCGGGGTTAACCGCGAACAGGCCCCGAATTACAACACAACCCCCAGGCACGGcccccccggg                                 >  2:4686891/1‑144 (MQ=255)
                                                                                                                        aGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCAAAAAGTAGTGCACCCCCCGAGGGGGCTGTGGTGAAACGGTGGGGGCGGGTTTAACCCCCAAGAGGCCAGTTCATAAAAAAAAACCCCCAGCCCCAGCCCCCCGGGGTAc                              >  2:5250941/1‑150 (MQ=255)
                                                                                                                          tGAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCAACAACAAATAGTTCACCCCCGGAGAGGGCGGTGGTGTACCATTTGGAGCGGGTTTACCCCCCAAGAGGCCCCTAAATAAAACACAAACCACACGCCCAGCCCCCCGGGGTACac                            >  2:4117570/1‑149 (MQ=255)
                                                                                                                          ggcgCCTGGGCCGGCCTCTGCGGGGGCAACCCCGCCCCACCACGGAGGGCCCCGCTTGAGCGGCCTGTGGTCACCCTGTTGTCGGGGTTTTCCCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAg                            <  1:7630735/147‑1 (MQ=255)
                                                                                                                           gAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGc                           <  2:8539305/150‑1 (MQ=255)
                                                                                                                               cTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACGTAGTGCCCCGCCTGCGAGGGCTGTGGTGAACCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGCGaaa                       <  2:7478487/150‑1 (MQ=255)
                                                                                                                               cTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCCCCGGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGTAGAGGCCCGGACTTACAACACACCCCTCAGGCCCTGCCCCGCGGGGTGCAGCGaaa                       >  2:5038906/1‑150 (MQ=37)
                                                                                                                                   gCCGGCCTCTTCGCGGGTAACCCCGCTCCTACACGTCGTGCACCGCCTGCGAGGGGTGTGGGGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGCGAAAAAcc                   <  2:3800027/150‑1 (MQ=39)
                                                                                                                                      GGCCTCTTCGGGGGTAACCCGGCTCCTCCACGTCGGGCACCGCGTGTGCGGGCTGTGGTGAACCTGTTGGCGGGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGCGAAAAACCGGa                >  1:7804738‑M1/1‑150 (MQ=255)
                                                                                                                                          TCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGCGAAAAACCGGACTCa            >  2:4782338‑M1/1‑150 (MQ=255)
                                                                                                                                           CTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGCGAAAAACCGGACTCAt           >  2:4167994‑M1/1‑150 (MQ=255)
                                                                                                                                                 CGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGCGAAAAACCGGACTCATAGAAc      <  1:4961192‑M1/150‑1 (MQ=255)
                                                                                                                                                     GTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGCGAAAAACCGGACTCATAGAACGCt   >  2:2266472‑M1/1‑150 (MQ=255)
                                                                                                                                                     GTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGCGAAAAACCGGACTCATAGAACGCt   <  1:877857‑M1/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GTGGCGAAGTGCAGACGGCGTTGAGCAAATTGCAGCGGTTCCGGCGGCAGAGCCGGGAGTCGTTGTTGTTGCAGAAGCCGGCTTCAATCAAGAAGGGTGGCTTCCTGAGTTGGTTGCACAAGTGAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTAGTGCACCGCCTGAGAGGGCTGTGGTGAACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTTACAACACAACCCTCAGGCACTGCCCCGCGTGGTACAGCGAAAAACCGGACTCATAGAACGCTG  >  NC_002947/6005774‑6006071

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: