Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,524,047 G→C 6.5% A251P (GCT→CCT) ‡ murG → UDP‑N‑acetylglucosamine‑‑N‑acetylmuramyl‑ (pentapeptide) pyrophosphoryl‑undecaprenol N‑acetylglucosamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,524,0470GC 6.5% 240.5 / 6.8 77A251P (GCT→CCT) ‡murGUDP‑N‑acetylglucosamine‑‑N‑acetylmuramyl‑ (pentapeptide) pyrophosphoryl‑undecaprenol N‑acetylglucosamine transferase
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base C (2/3);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

CCTAAGGCCCTGTCCGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGCCCAATATCTGGCCC  >  NC_002947/1523900‑1524191
                                                                                                                                                   |                                                                                                                                                 
ccTAAGGCCCTGTCCGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGc                                                                                                                                                  >  1:218341/1‑149 (MQ=255)
 cTAAGGCCCTGTCCGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCt                                                                                                                                                 >  7:116298/1‑149 (MQ=255)
    aGGCCCTGTCCGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGTTCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCtt                                                                                                                                                <  4:164548/147‑1 (MQ=255)
      gCCCTGTCCGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATgg                                                                                                                                            >  5:155340/1‑149 (MQ=255)
        ccTGTCCGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCt                                                                                                                                          >  8:42986/1‑149 (MQ=255)
           gTCCGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGTCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTAGCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAATGAGATGGCCGAGGCTTATGGCTgg                                                                                                                                        >  7:108411/1‑148 (MQ=255)
           gTCCGAGGTGACCGCGGCGCTGCTTGCAGAGGGTTTCCACGAGGCGGGCAAGGAGCACGCACCGATCACCGCCGATCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTgg                                                                                                                                        <  8:108411/148‑1 (MQ=255)
            tCCGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTggg                                                                                                                                       <  7:22742/148‑1 (MQ=255)
            tCCGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGc                                                                                                                                     >  4:124775/1‑149 (MQ=255)
              cGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCGCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGCCATAAGCCTGGGcc                                                                                                                                    <  8:285489/148‑1 (MQ=255)
                aGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGAc                                                                                                                                 <  3:291003/149‑1 (MQ=255)
                 ggTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGCCTGGGCCGAc                                                                                                                                 >  8:106755/1‑148 (MQ=255)
                  gTGCCCGCTGCGCTGCGGCCAGAGGTGTTCCACCAGTCCGGCAAGCATCACGCACCGATCACCGCCGAGCTTTATCACGAGGCGGGCGTCGCGGCTCATGTAGAGCCTTTCATCAAATACATGGCCCAGCCATAAGCCTGGGccatgtc                                                                                                                               <  4:218165/149‑6 (MQ=255)
                     cccGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCtggt                                                                                                                            <  7:42986/149‑1 (MQ=255)
                     cccGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCtgg                                                                                                                             <  3:124775/148‑1 (MQ=255)
                      ccGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCtggtg                                                                                                                           <  1:137367/149‑1 (MQ=255)
                        gcggcgCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTggtggt                                                                                                                         >  1:38383/1‑149 (MQ=255)
                         cggcgCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGgtggtg                                                                                                                        >  7:26920/1‑149 (MQ=255)
                          ggcgCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGgtgt                                                                                                                       <  6:310723/149‑1 (MQ=255)
                          ggcgCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGgtgt                                                                                                                       <  2:218341/149‑1 (MQ=255)
                          ggcgCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGgtgt                                                                                                                       <  7:216305/149‑1 (MQ=255)
                               tGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTgccgg                                                                                                                  >  3:106159/1‑149 (MQ=255)
                                   gCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACACGGCCCAGCCATAAGCCTGGGccat                                                                                                                                  >  3:156618/1‑127 (MQ=255)
                                   gCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACACGGCCCAGCCATAAGCCTGGGccat                                                                                                                                  <  4:156618/129‑3 (MQ=255)
                                    ccAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTgccgggccgg                                                                                                             >  7:187697/1‑149 (MQ=255)
                                      agagGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTg                                                                                                           >  4:211475/1‑149 (MQ=255)
                                        agGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGcc                                                                                                         >  2:248581/1‑149 (MQ=255)
                                         ggTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCt                                                                                                        <  2:28888/149‑1 (MQ=255)
                                             ttCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTTGAGCCTTTCATCAATGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGAc                                                                                                    <  7:106755/149‑1 (MQ=255)
                                             ttCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGa                                                                                                     <  4:106159/148‑1 (MQ=255)
                                               ccaccaGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTg                                                                                                  >  7:301628/1‑149 (MQ=255)
                                               ccaccaGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTg                                                                                                  >  3:146110/1‑149 (MQ=255)
                                                  ccaGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTct                                                                                               >  1:293890/1‑148 (MQ=255)
                                                      gCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCAGGGCCGGTGCCTTGACTGTCAGCg                                                                                            >  2:17536/1‑148 (MQ=255)
                                                           cAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTgg                                                                                      <  4:146110/149‑1 (MQ=255)
                                                             agcagcACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGgca                                                                                    <  8:26920/149‑1 (MQ=255)
                                                                  cacgcacCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGgcagcag                                                                                <  1:248581/149‑1 (MQ=255)
                                                                      cacCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGgc                                                                           <  5:328/149‑1 (MQ=255)
                                                                      cacCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGgc                                                                           <  2:38383/149‑1 (MQ=255)
                                                                        ccGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCTGGTGCCTTGACTGTCAGCTAACTGGCAGCAGCAGgcct                                                                         <  3:229890/149‑1 (MQ=255)
                                                                         agATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGgcctg                                                                        >  2:27042/2‑149 (MQ=255)
                                                                              aCCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCtt                                                                    <  8:301628/148‑1 (MQ=255)
                                                                              aCCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTc                                                                   <  6:155340/149‑1 (MQ=255)
                                                                              aCCGCCGAGCGTTATCACGAGGCGGGCGTCGCGCCTCAGCTAGAGCCATTCATCAAAGACATGGCACAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCACCAGCAGGCCTGCCTTc                                                                   >  1:44550/1‑149 (MQ=255)
                                                                                cgccgAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGa                                                                 <  5:263556/149‑1 (MQ=255)
                                                                                cgccgAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGa                                                                 >  4:286245/1‑149 (MQ=255)
                                                                                  ccgAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATg                                                               >  1:123001/1‑149 (MQ=255)
                                                                                       cGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGt                                                          >  3:176056/1‑149 (MQ=255)
                                                                                                gAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGcc                                                 <  2:229708/149‑1 (MQ=255)
                                                                                                  ggcgggcgTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGccc                                                >  1:189066/1‑148 (MQ=255)
                                                                                                     gggcgTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCAcg                                            >  1:269807/1‑149 (MQ=255)
                                                                                                          tCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATc                                       <  2:189066/149‑1 (MQ=255)
                                                                                                          tCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATc                                       <  2:123001/149‑1 (MQ=255)
                                                                                                             cgGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCCGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATcgac                                    <  5:305757/149‑1 (MQ=255)
                                                                                                             cgGCTCAGCTAGATCCTTTCATCAAAGACATGGCCCCGGCTTATGGCTTGGCCGACCTGGGGTGGGGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATcgac                                    <  8:230624/149‑1 (MQ=255)
                                                                                                                cTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGGCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGAt                                 <  3:275376/149‑1 (MQ=255)
                                                                                                                 tCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGAt                                 >  3:203699/1‑148 (MQ=255)
                                                                                                                 tCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGAt                                 >  7:216146/1‑148 (MQ=255)
                                                                                                                   aGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATc                                <  6:9515/147‑1 (MQ=255)
                                                                                                                    ggTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGNTGCCTTTGCCCCACGCGATCGACGATCAcc                             <  1:27042/149‑1 (MQ=255)
                                                                                                                           ccTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGa                                         <  1:145072/130‑1 (MQ=255)
                                                                                                                           ccTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGa                                         >  2:145072/1‑130 (MQ=255)
                                                                                                                           ccTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATcga                                     >  2:169035/1‑134 (MQ=255)
                                                                                                                           ccTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATcga                                     <  1:169035/134‑1 (MQ=255)
                                                                                                                           ccTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTca                      <  8:216146/149‑1 (MQ=255)
                                                                                                                             tttCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTcac                     >  6:246274/1‑148 (MQ=255)
                                                                                                                                 atcaAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTAGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGc                <  6:240208/149‑1 (MQ=255)
                                                                                                                                    aaaGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGccc              >  6:58808/1‑148 (MQ=255)
                                                                                                                                       gACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGACCGGTGCCTTGACTGTCAGCGAACTGGCAGCGGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGCCCAata          >  4:78760/1‑149 (MQ=255)
                                                                                                                                         gtcggCCCAGCCATAAGCCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCATGCGATCGACGATCACCAGACTCACAACGCCCAATATc        >  8:194893/4‑149 (MQ=255)
                                                                                                                                          aTGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTTCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGCCCAATATCt       <  8:116298/149‑1 (MQ=255)
                                                                                                                                          aTGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGCCCAatat         >  4:225271/1‑147 (MQ=255)
                                                                                                                                          aTGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGCCCAATATCt       <  8:187697/149‑1 (MQ=255)
                                                                                                                                          aTGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGCCCAATATCt       >  5:250077/1‑149 (MQ=255)
                                                                                                                                          aTGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGACCGGTGCCTTGACTGTCAGCGAACTGGCAGCGGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGCCCAATATCt       <  3:78760/149‑1 (MQ=255)
                                                                                                                                             gCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGCCCAATATCTGGc    <  3:211475/149‑1 (MQ=255)
                                                                                                                                               ccAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGCCCAATATCTGGccc  >  1:291840/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTAAGGCCCTGTCCGAGGTGCCCGCGGCGCTGCGGCCAGAGGTGTTCCACCAGGCCGGCAAGCAGCACGCACCGATCACCGCCGAGCGTTATCACGAGGCGGGCGTCGCGGCTCAGGTAGAGCCTTTCATCAAAGACATGGCCCAGGCTTATGGCTGGGCCGACCTGGTGGTGTGCCGGGCCGGTGCCTTGACTGTCAGCGAACTGGCAGCAGCAGGCCTGCCTTCGATGCTGGTGCCTTTGCCCCACGCGATCGACGATCACCAGACTCACAACGCCCAATATCTGGCCC  >  NC_002947/1523900‑1524191

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