Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 134,397 A→C 5.4% D98A (GAT→GCT) ‡ PP_0129 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947134,3970AC 5.4% 292.7 / 6.4 93D98A (GAT→GCT) ‡PP_0129GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (46/42);  new base C (3/2);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

GGTGGTAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGCGGC  >  NC_002947/134251‑134542
                                                                                                                                                   |                                                                                                                                                 
ggtggtAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATg                                                                                                                                                 >  4:362974/1‑149 (MQ=255)
 gtggtAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATgc                                                                                                                                                >  1:167873/1‑149 (MQ=255)
 gtggtAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATgcg                                                                                                                                               >  4:367046/1‑149 (MQ=255)
    gtACACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCCACAACATGGACGCCGGCCTCGCTCGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCt                                                                                                                                             >  4:136322/1‑149 (MQ=255)
       gACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGa                                                                                                                                           <  3:369367/148‑1 (MQ=255)
       gACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGAc                                                                                                                                          >  3:343043/1‑149 (MQ=255)
            gTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCAc                                                                                                                                     >  8:85473/1‑149 (MQ=255)
                aTGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGt                                                                                                                                 >  8:27377/1‑149 (MQ=255)
                aTGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGt                                                                                                                                 >  7:295276/1‑149 (MQ=255)
                 tGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGtt                                                                                                                                >  4:332069/1‑149 (MQ=255)
                   aaGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTgg                                                                                                                                  >  1:128729/1‑145 (MQ=255)
                   aaGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTgg                                                                                                                                  <  2:128729/145‑1 (MQ=255)
                   aaGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTgg                                                                                                                              >  1:326025/1‑149 (MQ=255)
                   aaGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTg                                                                                                                               <  8:278430/148‑1 (MQ=255)
                    aGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGc                                                                                                                             <  8:202796/149‑1 (MQ=255)
                     ggAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGc                                                                                                                             >  3:104914/1‑148 (MQ=255)
                            cgcgATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGc                                                                                                                      <  1:26955/148‑1 (MQ=255)
                             gcgATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCaa                                                                                                                    <  2:167873/149‑1 (MQ=255)
                              cgATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAg                                                                                                                   <  7:281224/149‑1 (MQ=255)
                                  gTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcct                                                                                                               >  4:261544/1‑149 (MQ=255)
                                  gTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcct                                                                                                               >  3:129985/1‑149 (MQ=255)
                                    cATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGCATCGCTGAcctgg                                                                                                                                      <  1:3639/124‑4 (MQ=255)
                                    cATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGCATCGCTGAcctgg                                                                                                                                      >  2:3639/1‑121 (MQ=255)
                                     aTTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGc                                                                                                                      <  7:53737/139‑1 (MQ=255)
                                     aTTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgcc                                                                                                            >  8:86734/1‑149 (MQ=255)
                                     aTTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAAAAGCGTAGCGAAACCCGCATGAACCAGGTCAGCGATGCGATGACCGCACTGGTTGGCCAAATga                                                                                                                      >  8:53737/1‑138 (MQ=255)
                                      ttCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgccc                                                                                                           >  6:59110/1‑149 (MQ=255)
                                              gACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGc                                                                                                                 <  8:139488/135‑1 (MQ=255)
                                              gACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGc                                                                                                                 >  7:139488/1‑135 (MQ=255)
                                                   caTGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGata                                                                                              >  1:365753/1‑149 (MQ=255)
                                                    aTGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGatat                                                                                             >  8:75821/1‑149 (MQ=255)
                                                    aTGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGatat                                                                                             >  7:353587/1‑149 (MQ=255)
                                                         cGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTcgc                                                                                        >  2:108805/1‑149 (MQ=255)
                                                          gCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTcgcg                                                                                       <  2:365753/149‑1 (MQ=255)
                                                            cGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgc                                                                                                             >  7:249047/1‑125 (MQ=255)
                                                            cGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGcctgc                                                                                                             <  8:249047/125‑1 (MQ=255)
                                                             ggCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGcc                                                                                    <  2:326025/149‑1 (MQ=255)
                                                               ccTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGACAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCag                                                                                  >  4:293747/1‑148 (MQ=255)
                                                                 tCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCAtt                                                                                 >  4:83735/1‑148 (MQ=255)
                                                                 tCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCAtt                                                                                 >  7:255382/1‑148 (MQ=255)
                                                                      ggCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGaaaa                                                                            >  4:88280/1‑148 (MQ=255)
                                                                        cTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAAc                                                                         <  1:139197/149‑1 (MQ=255)
                                                                         ttGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACt                                                                        >  1:190436/1‑149 (MQ=255)
                                                                           gATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTg                                                                       <  3:332069/148‑1 (MQ=255)
                                                                            aTTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTgg                                                                      <  8:255382/148‑1 (MQ=255)
                                                                            aTTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGc                                                                     >  6:189609/1‑149 (MQ=255)
                                                                             ttCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGc                                                                     >  3:200473/1‑148 (MQ=255)
                                                                               ccGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCga                                                                   <  2:27198/148‑1 (MQ=255)
                                                                                cgcgCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCgaag                                                                 <  3:136322/149‑1 (MQ=255)
                                                                                 gcgcCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCgaag                                                                 >  3:116450/1‑148 (MQ=255)
                                                                                     cTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCt                                                            >  3:183027/1‑149 (MQ=255)
                                                                                         agaAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGa                                                         <  2:190436/148‑1 (MQ=255)
                                                                                          gaAAGCGGTACTGGATACCTAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAgg                                                       <  2:334220/149‑1 (MQ=255)
                                                                                                ggTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCgggg                                                 <  4:314785/149‑1 (MQ=255)
                                                                                                       gATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGccc                                           <  3:362974/148‑1 (MQ=255)
                                                                                                       gATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGccc                                           <  8:219958/148‑1 (MQ=255)
                                                                                                       gATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGcc                                            >  4:269621/1‑147 (MQ=255)
                                                                                                       gATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCa                                          <  8:353587/149‑1 (MQ=255)
                                                                                                       gATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCa                                          <  3:83735/149‑1 (MQ=255)
                                                                                                         ttCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGt                                        >  8:333214/1‑149 (MQ=255)
                                                                                                           cctaaCAGCGTCGCGAAACCCGTATTAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTc                                       <  7:296503/145‑1 (MQ=255)
                                                                                                           ccGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTc                                       <  4:343043/148‑1 (MQ=255)
                                                                                                           ccGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTc                                       <  4:99743/148‑1 (MQ=255)
                                                                                                            cGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAc                                     <  3:367046/149‑1 (MQ=255)
                                                                                                            cGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAc                                     <  4:104914/149‑1 (MQ=255)
                                                                                                             gAACAGCGTCGCGAAACCCGCATGAACCAGGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAc                                     >  1:33676/1‑149 (MQ=255)
                                                                                                             gAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAcg                                    >  1:368030/1‑149 (MQ=255)
                                                                                                             gAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAcg                                    >  3:3837/1‑149 (MQ=255)
                                                                                                             gAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAcg                                    >  1:16655/1‑149 (MQ=255)
                                                                                                              aaCAGGGTCGCGAAACCCTCATGAACCAGGTCAGCGATGCGCTGACCGAACTGGTTGGCAAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAATCTTGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAcgc                                   <  3:293747/149‑1 (MQ=255)
                                                                                                              aaCAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCAcgc                                   >  7:240829/1‑149 (MQ=255)
                                                                                                                cAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGa                                 <  2:33676/149‑1 (MQ=255)
                                                                                                                     tCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTgc                            >  2:112954/1‑149 (MQ=255)
                                                                                                                      cgcgAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTgcc                           >  2:140341/1‑149 (MQ=255)
                                                                                                                           aaCCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGt                                        <  1:162986/131‑1 (MQ=255)
                                                                                                                           aaCCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGt                                        >  2:162986/1‑131 (MQ=255)
                                                                                                                            aCCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTTAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCTCCGCt                     >  7:58972/1‑149 (MQ=255)
                                                                                                                              ccGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGt                   >  3:153467/1‑149 (MQ=255)
                                                                                                                               cGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGtt                  >  2:115582/1‑149 (MQ=255)
                                                                                                                               cGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGtt                  <  4:183027/149‑1 (MQ=255)
                                                                                                                               cGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGtt                  <  1:43226/149‑1 (MQ=255)
                                                                                                                                   tGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGc              <  3:269621/149‑1 (MQ=255)
                                                                                                                                     aaCCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGa            <  7:85473/149‑1 (MQ=255)
                                                                                                                                      agtgcggTCAGCGCATCGCTGACCTCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGa                                                         <  8:219494/99‑1 (MQ=255)
                                                                                                                                      agtgcggTCAGCGCATCGCTGACCTCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGa                                                         >  7:219494/6‑104 (MQ=255)
                                                                                                                                      aCCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAg           <  7:75821/149‑1 (MQ=255)
                                                                                                                                      aCCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAg           <  5:59110/149‑1 (MQ=255)
                                                                                                                                        cAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGtt         >  6:134898/1‑149 (MQ=255)
                                                                                                                                          cggTCAGCGCATCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGtt         >  2:158771/2‑147 (MQ=255)
                                                                                                                                            gTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAg      <  4:200473/149‑1 (MQ=255)
                                                                                                                                             tCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGTTGCCGCCGCTGTTGGGCGAGTTGAGc     <  7:27377/149‑1 (MQ=255)
                                                                                                                                               aGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGCg    <  4:116450/148‑1 (MQ=255)
                                                                                                                                                gCGATGCGCGGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGCGGc  <  8:295276/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGTGGTAGACGAGTGCATGAAGGAAATGCGCGATGTCATTCGCAGCGACAACATGGACGCCGGCCTCGCTGGCTTGATTCCGCGCCTTGAGAAAGCGGTACTGGATTCCGAACAGCGTCGCGAAACCCGCATGAACCAGGTCAGCGATGCGCTGACCGCACTGGTTGGCCAACTGCAAGGCCTGCCCCTGCCCAGCGATATCTCGCGGCCATTGAAAAAACTGGCGAAGAAGCTCGAGGGCGGGGTTGCCCAGTCACGCGAGCTGCCGCCGCTGTTGGGCGAGTTGAGCGGC  >  NC_002947/134251‑134542

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