Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,173,519 T→C 6.4% N374D (AAT→GAT)  glmU ← bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,173,5190TC 6.4% 274.6 / 6.5 93N374D (AAT→GAT) glmUbifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base T (41/46);  new base C (3/3);  total (44/49)
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.26e-01

TCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAAC  >  NC_002947/6173372‑6173664
                                                                                                                                                   |                                                                                                                                                 
tCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTg                                                                                                                                                  <  1:19394/149‑1 (MQ=255)
tCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTg                                                                                                                                                  <  3:211205/149‑1 (MQ=255)
tCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTg                                                                                                                                                  >  6:232704/1‑149 (MQ=255)
tCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTg                                                                                                                                                  <  6:357410/149‑1 (MQ=255)
 cGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTgg                                                                                                                                                 >  5:8953/1‑149 (MQ=255)
  gAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGt                                                                                                                                                <  4:88122/149‑1 (MQ=255)
  gAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGt                                                                                                                                                <  7:298303/149‑1 (MQ=255)
  gAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGt                                                                                                                                                <  6:385590/149‑1 (MQ=255)
   aGCCGGCTGCGGTGGTCGCGCCAGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTg                                                                                                                                               >  7:125457/1‑149 (MQ=255)
                  tCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCg                                                                                                                                >  4:147591/1‑149 (MQ=255)
                  tCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAAAACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGAGCGCACCGATTTCg                                                                                                                                >  6:301526/1‑149 (MQ=255)
                  tCGCGCAGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCg                                                                                                                                >  6:234182/1‑149 (MQ=255)
                    gcgcCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGAAATCGGTGCGCGCACCAATATCGGc                                                                                                                              <  1:335359/149‑1 (MQ=255)
                     cgcCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCCCCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGAGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCa                                                                                                                             <  1:278087/149‑1 (MQ=255)
                     cgcCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGc                                                                                                                              >  8:321386/1‑148 (MQ=255)
                      gcCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCAt                                                                                                                            >  5:273217/1‑149 (MQ=255)
                        cGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGGCTTCACCCATCACGGTCTTGCACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATcc                                                                                                                          <  4:69712/149‑1 (MQ=255)
                          gCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATcccc                                                                                                                        <  6:8953/149‑1 (MQ=255)
                            cTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATcccc                                                                                                                        >  6:18413/1‑147 (MQ=255)
                             ttGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCag                                                                                                                     >  5:343058/1‑149 (MQ=255)
                               gATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCagg                                                                                                                    >  5:116776/1‑148 (MQ=255)
                                aTTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCaggta                                                                                                                  <  1:149077/149‑1 (MQ=255)
                                aTTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCaggt                                                                                                                   <  3:96472/148‑1 (MQ=255)
                                 tttCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCaggtag                                                                                                                 <  5:232704/149‑1 (MQ=255)
                                 tttCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCaggta                                                                                                                  <  7:275744/148‑1 (MQ=255)
                                  ttCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCaggtagg                                                                                                                >  5:301669/1‑149 (MQ=255)
                                  ttCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCaggtagg                                                                                                                >  7:15883/1‑149 (MQ=255)
                                     cacaGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTc                                                                                                              <  5:18413/148‑1 (MQ=255)
                                           cgcCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGAc                                                                                                       <  2:248127/149‑1 (MQ=255)
                                            gcCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGAcc                                                                                                      <  7:332464/149‑1 (MQ=255)
                                                 caGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCGATCGTAGTTGCAGGTGATGGTGCCGGCGCTGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCt                                                                                                 >  8:99190/1‑149 (MQ=255)
                                                 caGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAgg                                                                                                           <  1:15358/139‑1 (MQ=255)
                                                 caGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAgg                                                                                                           >  2:15358/1‑139 (MQ=255)
                                                 caGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCt                                                                                                 <  6:301669/149‑1 (MQ=255)
                                                  aGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGAc                                                                                                       >  1:145529/1‑142 (MQ=255)
                                                  aGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGAc                                                                                                       <  2:145529/142‑1 (MQ=255)
                                                        ttgttCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGcccc                                                                                          <  1:23244/149‑1 (MQ=255)
                                                          gttCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCtt                                                                                        >  1:231344/1‑149 (MQ=255)
                                                           ttCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATCGGTGCGCGCACCAATATc                                                                                                                                 >  6:24103/1‑107 (MQ=255)
                                                           ttCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATCGGTGCGCGCACCAATATc                                                                                                                                 <  5:24103/107‑1 (MQ=255)
                                                             cGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCGCGCACCAATATCGGCGCCGGCAccatcac                                                                                                                                  >  6:114367/1‑99 (MQ=25)
                                                             cGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCGCGCACCAATATCGGCGCCGGCAccatcac                                                                                                                                  <  5:114367/104‑6 (MQ=25)
                                                              gAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCAcc                                                                                    >  1:212339/1‑149 (MQ=255)
                                                               aaCCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCAccc                                                                                   <  5:673/149‑1 (MQ=255)
                                                                     tGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCTCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTgg                                                                             <  6:116776/149‑1 (MQ=255)
                                                                     tGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTgg                                                                             <  8:15883/149‑1 (MQ=255)
                                                                      gAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTggg                                                                            >  8:211585/1‑149 (MQ=255)
                                                                         cacGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGt                                                                         >  8:302699/1‑149 (MQ=255)
                                                                         cacGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGAGCCGATATTGGTGCGCGCCACGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGt                                                                         >  3:224768/1‑149 (MQ=255)
                                                                          acGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGtt                                                                        >  7:168261/1‑149 (MQ=255)
                                                                          acGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGtt                                                                        >  8:268991/1‑149 (MQ=255)
                                                                           cGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGttt                                                                       <  3:358339/149‑1 (MQ=255)
                                                                            gTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGtttt                                                                      >  5:246168/1‑149 (MQ=255)
                                                                            gTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGtttt                                                                      >  3:114258/1‑149 (MQ=255)
                                                                             tCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGttttt                                                                     <  5:234182/149‑1 (MQ=255)
                                                                             tCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGttttt                                                                     <  5:301526/149‑1 (MQ=255)
                                                                               ttCACCCATCACTGTCTTTAACTTGTAGGCGCCATAGTAGTTGCAGGTGATGGTGCCGGAGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCa                                                                   <  4:224768/149‑1 (MQ=255)
                                                                                  aCCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGc                                                                                             <  1:206169/120‑1 (MQ=255)
                                                                                  aCCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGc                                                                                             >  2:206169/1‑120 (MQ=255)
                                                                                     cATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCaa                                                             >  6:136761/1‑149 (MQ=255)
                                                                                         aCGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGaa                                                         <  7:268991/149‑1 (MQ=255)
                                                                                         aCGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGaa                                                         <  3:147591/149‑1 (MQ=255)
                                                                                           ggTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGt                                                       >  4:324231/1‑149 (MQ=255)
                                                                                            ctcatGAACTTGTTGCCGCCATCGTACTTGCAGTTGATGGTGCCGGCGCCGATATTGGTGCGCTCACCGATTTCGGCATCCCCCATGTAGGGCAGGTGACCAGCCTTGTCCCCTTCACCCAGGTGGTCGTTTTTCAGTTCAACGAAGtt                                                      <  1:273685/145‑1 (MQ=255)
                                                                                             tCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGttttt                                                                     >  5:354650/1‑133 (MQ=255)
                                                                                             tCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGttttt                                                                     <  6:354650/133‑1 (MQ=255)
                                                                                                tGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAg                                                                                 >  4:22528/1‑118 (MQ=255)
                                                                                                tGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAg                                                                                 <  3:22528/118‑1 (MQ=255)
                                                                                                tGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTAcc                                                   >  6:226379/1‑148 (MQ=255)
                                                                                                  aaCTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCca                                                 >  6:29542/1‑148 (MQ=255)
                                                                                                    cTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACAt                                              <  2:231344/149‑1 (MQ=255)
                                                                                                       gttgGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATgg                                            <  8:125457/148‑1 (MQ=255)
                                                                                                       gttgGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATgg                                            <  7:321386/148‑1 (MQ=255)
                                                                                                       gttgGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATgg                                            <  6:246168/148‑1 (MQ=255)
                                                                                                            cgcCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTg                                      <  5:136761/149‑1 (MQ=255)
                                                                                                            cgcCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTg                                      <  2:222604/149‑1 (MQ=255)
                                                                                                                 tCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTc                                 <  6:273217/149‑1 (MQ=255)
                                                                                                                        tGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACg                          >  7:90530/1‑149 (MQ=255)
                                                                                                                           aGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTg                       >  5:105681/1‑149 (MQ=255)
                                                                                                                           aGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTg                       >  4:33525/1‑149 (MQ=255)
                                                                                                                            ggTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTc                                 <  8:266352/138‑1 (MQ=255)
                                                                                                                            ggTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTc                                 >  7:266352/1‑138 (MQ=255)
                                                                                                                               gATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCa                                                                   >  1:370077/1‑101 (MQ=255)
                                                                                                                               gATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCa                                                                   <  2:370077/101‑1 (MQ=255)
                                                                                                                                  ggTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGgc                >  2:191134/1‑149 (MQ=255)
                                                                                                                                   gTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCTTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGgcg               >  5:8848/1‑149 (MQ=255)
                                                                                                                                   gTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGgc                <  5:29542/148‑1 (MQ=255)
                                                                                                                                   gTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGgc                <  4:114258/148‑1 (MQ=255)
                                                                                                                                         atgcCGAAATCGGTGCGCGCACCAATATCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCg         >  3:357060/3‑149 (MQ=255)
                                                                                                                                          cgcCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGa        <  6:343058/149‑1 (MQ=255)
                                                                                                                                          cgcCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGGTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGa        <  2:212339/149‑1 (MQ=255)
                                                                                                                                              gATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCa     >  5:42503/1‑148 (MQ=255)
                                                                                                                                                tatTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAAc  >  1:343568/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAAC  >  NC_002947/6173372‑6173664

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: