Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,693,834 C→G 6.9% R87R (CGC→CGG) ‡ PP_1489 → CheW domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,693,8340CG 6.9% 270.0 / 10.2 87R87R (CGC→CGG) ‡PP_1489CheW domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (41/40);  new base G (3/3);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

AAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGATGC  >  NC_002947/1693686‑1693983
                                                                                                                                                      |                                                                                                                                                     
aaGTGTTGCCGCGTCAGCCGCTCAAGCCGATCACCCATGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGcc                                                                                                                                                      >  1:627826/1‑150 (MQ=255)
     ttGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCctggtg                                                                                                                                                 >  2:1333248/1‑150 (MQ=255)
        ccGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGTCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCAGGcgggtgctg                                                                                                                                              <  2:59224/150‑1 (MQ=255)
         cgcgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTTCAACGCACCA‑GCACCAGGcgggtgctgg                                                                                                                                             <  1:572015/150‑1 (MQ=255)
          gcgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCAGGcgggtgctggt                                                                                                                                            <  2:1515064/150‑1 (MQ=255)
            cgtgcacGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCtggtgctggtgc                                                                                                                                          <  2:372781/144‑1 (MQ=255)
                gCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTa                                                                                                                                      >  1:723924/1‑150 (MQ=255)
                   gCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCg                                                                                                                                   <  1:954886/150‑1 (MQ=255)
                    cTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGa                                                                                                                                  >  2:1217047/1‑150 (MQ=255)
                     tCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAg                                                                                                                                 >  2:803529/1‑150 (MQ=255)
                     tCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAg                                                                                                                                 <  1:94276/150‑1 (MQ=255)
                     tCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAg                                                                                                                                 >  2:803613/1‑150 (MQ=255)
                     tCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATAACGTTTGC‑CGCCTTGAACTTTGGCCAGCAAGTCGGCAAAGTCATCA‑GGATCGGCCTGGTGCTGGTGCACTACCGAg                                                                                                                                 >  1:424364/1‑150 (MQ=255)
                       aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTa                                                                                                                                      >  1:860116/1‑143 (MQ=255)
                       aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTa                                                                                                                                      <  2:860116/143‑1 (MQ=255)
                       aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAggc                                                                                                                               >  2:1143539/1‑150 (MQ=255)
                       aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAggc                                                                                                                               >  2:711642/1‑150 (MQ=255)
                       aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAggc                                                                                                                               >  1:734878/1‑150 (MQ=255)
                        aggcggtCGCCCAGGCACCTGCCTGGGTCGCGGGACTTCTCGCACATCGTTGCTAGCTGGTTCCGGTCCAC‑ACCAGGCGCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAggcg                                                                                                                              <  1:424945/144‑1 (MQ=255)
                        aGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAggcg                                                                                                                              <  1:830909/150‑1 (MQ=255)
                            gATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCt                                                                                                                          <  1:24676/150‑1 (MQ=255)
                            gATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCt                                                                                                                          >  2:1170074/1‑150 (MQ=255)
                                gCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAg                                                                                                                      >  1:434185/1‑150 (MQ=255)
                                gCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAg                                                                                                                      <  1:1324550/150‑1 (MQ=255)
                                  ctaGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCt                                                                                                                    >  1:141816/3‑150 (MQ=255)
                                  ccAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGcc                                                                                                                    >  2:981805/1‑149 (MQ=255)
                                  ccAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCt                                                                                                                    >  1:589214/1‑150 (MQ=255)
                                   cAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCtt                                                                                                                   <  1:315500/150‑1 (MQ=255)
                                    aGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTg                                                                                                                  <  1:1217047/150‑1 (MQ=255)
                                      gCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGc                                                                                                                <  1:803613/150‑1 (MQ=255)
                                      gCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGc                                                                                                                <  1:803529/150‑1 (MQ=255)
                                            gcctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGATAGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATc                                                                                                          <  1:628207/150‑1 (MQ=255)
                                                gggTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTgg                                                                                                      >  1:1392265/1‑150 (MQ=255)
                                                 ggTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGcc                                                                                                                                                      >  1:1331619/1‑101 (MQ=255)
                                                 ggTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGcc                                                                                                                                                      <  2:1331619/101‑1 (MQ=255)
                                                     gcgGGCATTCTCGCCCATCGTGGCCAGCTGGTTCTCCTCAGCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAg                                                                                                 <  2:627826/150‑1 (MQ=255)
                                                     gcgGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑AGCCTTGACATTTTGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAg                                                                                                 >  1:232765/1‑150 (MQ=255)
                                                       gggAATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGc                                                                                               <  1:968938/150‑1 (MQ=255)
                                                        ggCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGcc                                                                                              >  1:234977/1‑150 (MQ=255)
                                                           aTTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACATTGC‑AGAATTGAACTTTGGGCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCAcc                                                                                           >  1:515219/1‑150 (MQ=255)
                                                            ttCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTa                                                                                                                                      <  2:1327625/106‑1 (MQ=255)
                                                            ttCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTa                                                                                                                                      >  1:1327625/1‑106 (MQ=255)
                                                               tcGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTGGGCCTGATCCTGGCGAAGGCGACCGacc                                                                                       >  2:782052/1‑148 (MQ=255)
                                                                  cccATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAAccc                                                                                    >  1:1399146/1‑150 (MQ=255)
                                                                    cATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTACAGCTTGTCCTGATCCTGGAGCAGGCCACCGAAACCCTg                                                                                  >  2:1258764/1‑150 (MQ=255)
                                                                     aTCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGc                                                                                 <  2:1571484/150‑1 (MQ=255)
                                                                     aTCGTGGCCAGCTGGTTCCGGTCATCCACCATCC‑AGCCCAACGCCGGCGGCAGGGGGCGGTGCAACGCAAAG‑GCATTCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGc                                                                                 <  1:1867417/150‑1 (MQ=255)
                                                                       cGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGt                                                                               <  2:889909/150‑1 (MQ=255)
                                                                        gTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGtt                                                                              <  2:1659678/150‑1 (MQ=255)
                                                                         tGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTg                                                                             <  1:1143539/150‑1 (MQ=255)
                                                                                  gggtCCGGTCATCGACCGTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGa                                                                    <  2:722434/147‑1 (MQ=255)
                                                                                   gTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGAc                                                                   <  2:234977/150‑1 (MQ=255)
                                                                                    ttCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACg                                                                  >  1:1215295/1‑150 (MQ=255)
                                                                                      ccGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCACCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAACCTCCCATTATCCTGGATCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAg                                                                <  1:516084/150‑1 (MQ=255)
                                                                                        ggTCATTGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAA‑GCGTCATTGACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTTGCCTGATCCTGGAGCAGGCCACCAAAACCCTGCGATGCCAACCGGCCGAGtt                                                              >  1:987421/1‑150 (MQ=255)
                                                                                         gTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTc                                                             <  1:1238637/150‑1 (MQ=255)
                                                                                            aTCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAg                                                          >  2:1517683/1‑150 (MQ=255)
                                                                                             tCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGCCGAGTTCCAGc                                                         >  2:232836/1‑150 (MQ=255)
                                                                                              cgcccTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGcc                                                        <  1:283091/147‑1 (MQ=255)
                                                                                               gACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGccc                                                       <  2:232765/150‑1 (MQ=255)
                                                                                                    tGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATgg                                                  >  2:1546532/1‑150 (MQ=255)
                                                                                                         gccTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTg                                              >  1:413745/1‑150 (MQ=255)
                                                                                                           cTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGTACTACCGAGGCGGCCTGc                                                                                                                        >  2:1546518/1‑74 (MQ=255)
                                                                                                           cTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGTACTACCGAGGCGGCCTGc                                                                                                                        <  1:1546518/74‑1 (MQ=255)
                                                                                                           cTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGa                                            <  1:1418655/150‑1 (MQ=255)
                                                                                                             tGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGaca                                          <  2:1392265/150‑1 (MQ=255)
                                                                                                                gCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAg                                       >  2:1279028/1‑150 (MQ=255)
                                                                                                                    tGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGGGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACAACTTCCTGCCCTATGGGGTGGACACAGGCGa                                   <  1:519821/150‑1 (MQ=255)
                                                                                                                    tGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCtt                                                                                                                   <  1:1174095/70‑1 (MQ=255)
                                                                                                                    tGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCtt                                                                                                                   >  2:1174095/1‑70 (MQ=255)
                                                                                                                    tGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGa                                   <  2:589214/150‑1 (MQ=255)
                                                                                                                      gCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATg                                 >  2:93350/1‑150 (MQ=255)
                                                                                                                        cAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATgcg                               >  2:517956/1‑150 (MQ=255)
                                                                                                                         aGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATggg                                                 <  2:242480/131‑1 (MQ=255)
                                                                                                                         aGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATggg                                                 >  1:242480/1‑131 (MQ=255)
                                                                                                                         aGCCCGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTCTGGGGTGGACACAGGCGATgcgc                              <  2:947351/150‑1 (MQ=255)
                                                                                                                                  tGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGGTTAGTCTGAGGCTGTAGTATGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTgg                     <  1:1170074/150‑1 (MQ=255)
                                                                                                                                    cAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGc                   <  2:434185/150‑1 (MQ=255)
                                                                                                                                     agtgcacca‑gcaccaGGCGGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAATCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGcc                  >  1:428990/4‑150 (MQ=255)
                                                                                                                                      aCGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGccc                 >  1:1448728/1‑150 (MQ=255)
                                                                                                                                      aCGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGccc                 >  1:872433/1‑150 (MQ=255)
                                                                                                                                      aCGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGccc                 <  1:667113/150‑1 (MQ=255)
                                                                                                                                          aCCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTCAGGGGGCGACACAGGCGATGCGCCTTATCTGGGCCCGGTg             >  2:1589763/1‑150 (MQ=255)
                                                                                                                                          aCCA‑GCACCAGGCGGGTGCTGGTGCGCTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCCGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCACCATTGGGGGGACACAGGCGATGCGCCTTATCTGGGCCCGGTg             >  1:1644835/1‑150 (MQ=255)
                                                                                                                                           ccA‑GCACCAGGCGGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTgc            >  1:1710817/1‑150 (MQ=255)
                                                                                                                                            cA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTgcg           >  1:1172083/1‑150 (MQ=255)
                                                                                                                                                 aCCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAg       <  2:734878/150‑1 (MQ=255)
                                                                                                                                                   ccGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGa     >  1:179178/1‑150 (MQ=255)
                                                                                                                                                   ccGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGa     >  1:1606059/1‑150 (MQ=255)
                                                                                                                                                    cGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGAt    <  2:424364/150‑1 (MQ=255)
                                                                                                                                                      ccTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCTTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGATGc  <  1:1333248/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGC‑CGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCA‑GCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGATGC  >  NC_002947/1693686‑1693983

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: ‑