Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,173,519 T→C 5.6% 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 5.6% 214.0 / 2.1 71N374D (AAT→GAT) glmUbifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base C (2/2);  total (36/35)
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.43e-01

GAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGC  >  NC_002947/6173374‑6173656
                                                                                                                                                  |                                                                                                                                         
gAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGt                                                                                                                                        <  6:204257/149‑1 (MQ=255)
          cGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCg                                                                                                                              <  6:100871/149‑1 (MQ=255)
            gtggtCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATCAACACGTCTTCACCCATCACGCTCTTGAACTTCTTGGCGCCATCGTAGTTGCAGCTGATGGTGCCGGCTCCGATATTGGTGCGCGCACCGAt                                                                                                                            >  6:41717/1‑149 (MQ=255)
                tCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCg                                                                                                                        <  3:107614/149‑1 (MQ=255)
                 cgcgCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGAAATCGGTGCGCGCACCAATATCgg                                                                                                                       <  7:165306/149‑1 (MQ=255)
                 cgcgCAGGACTTGTTTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAAGACGGCTTCACCCATCACGGTCTTTATCTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGGCGATATTGGTGCGCGCACCGTTTTCgg                                                                                                                       <  7:283186/149‑1 (MQ=255)
                  gcgcCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGc                                                                                                                      <  6:157364/149‑1 (MQ=255)
                     ccGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATc                                                                                                                   <  8:35581/149‑1 (MQ=255)
                      cGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATcc                                                                                                                  >  8:107834/1‑149 (MQ=255)
                            tGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCTCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCagg                                                                                                            >  5:47425/1‑149 (MQ=255)
                             gATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCagg                                                                                                            >  6:270053/1‑148 (MQ=255)
                              aTTTCCACAGGTGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCTATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCagg                                                                                                            >  7:180566/1‑147 (MQ=255)
                                 tCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCaggtaggt                                                                                                       >  6:232851/1‑149 (MQ=255)
                                   cacaGGCGGCACCAGCGAGTTGTTCGAACCGATGATCACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCATGTAGGTc                                                                                                      <  5:341622/148‑1 (MQ=255)
                                    acaGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAg                                                                                                    >  6:319734/1‑149 (MQ=255)
                                      aGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGt                                                                                                  <  7:174777/149‑1 (MQ=255)
                                      aGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGt                                                                                                  <  7:24039/149‑1 (MQ=255)
                                        gcgcCACCAGCGAGTTGTTCGAACCGATGAACATCGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCag                                                                                                             >  1:218414/1‑137 (MQ=255)
                                        gcgcCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGa                                                                                                >  6:196337/1‑149 (MQ=255)
                                               caGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCt                                                                                         >  4:210693/1‑149 (MQ=255)
                                                 gCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTg                                                                                       <  4:5224/149‑1 (MQ=255)
                                                 gCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGAAATCGGTGCGCGCACCAATATCGGc                                                                                                                      <  5:315387/118‑1 (MQ=39)
                                                 gCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGAAATCGGTGCGCGCACCAATATCGGc                                                                                                                      >  6:315387/1‑118 (MQ=39)
                                                   gtgttgttCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGTTGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGc                                                                                     <  5:41717/147‑1 (MQ=255)
                                                         ttCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCtt                                                                                >  4:281151/1‑148 (MQ=255)
                                                           cGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGc                                                                                     >  2:64447/1‑141 (MQ=255)
                                                           cGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGc                                                                                     <  1:64447/141‑1 (MQ=255)
                                                           cGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCAc                                                                             >  1:209067/1‑149 (MQ=255)
                                                             aaCCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCTCCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCAcc                                                                            >  4:185167/1‑148 (MQ=255)
                                                             aaCCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGAc                                                                                               <  2:231880/129‑1 (MQ=255)
                                                             aaCCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGAc                                                                                               >  1:231880/1‑129 (MQ=255)
                                                                cGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTc                                                                                                      <  7:75813/119‑1 (MQ=255)
                                                                cGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTc                                                                                                      >  8:75813/1‑119 (MQ=255)
                                                                      acacGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCg                                                                  >  5:298618/1‑149 (MQ=255)
                                                                        acGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGtt                                                                >  8:121386/1‑149 (MQ=255)
                                                                           gTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGAGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCACAGGTAGGTCAGGTGCCCAGCCTTGGCCCCTTCACACAGGTGGGCGtttt                                                              >  6:30267/1‑149 (MQ=255)
                                                                                cACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCa                                                           <  5:282337/147‑1 (MQ=255)
                                                                                 aCCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGCGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGtt                                                        >  2:18646/1‑149 (MQ=255)
                                                                                  cccATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGtt                                                        >  4:134047/1‑148 (MQ=255)
                                                                                   ccATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTc                                                       >  7:184769/1‑148 (MQ=255)
                                                                                    cATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCaa                                                     <  6:298618/149‑1 (MQ=255)
                                                                                     aTCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAAc                                                    <  1:24847/149‑1 (MQ=255)
                                                                                        aCGGTCTTGAACTTGTTGGCGGCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCACCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGaa                                                 <  5:30267/149‑1 (MQ=255)
                                                                                        aCGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGTTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCACAAGGTAGGTCAGGTGACCAGCATTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCaaag                                                   >  7:305232/1‑145 (MQ=255)
                                                                                        aCGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACg                                                   >  5:128650/1‑147 (MQ=255)
                                                                                        aCGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACg                                                   >  4:155556/1‑147 (MQ=255)
                                                                                          ggTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGt                                               >  2:112170/1‑149 (MQ=255)
                                                                                           gTCTTGAACTTGTTGGCACCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGc                                                                                     <  7:27702/110‑1 (MQ=255)
                                                                                           gTCTTGAACTTGTTGGCACCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGCCATCCCCCAGGTAGGTCAGGTCACCAGCATAGGc                                                                                     >  8:27702/1‑110 (MQ=255)
                                                                                             cTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTAc                                            >  1:218891/1‑149 (MQ=255)
                                                                                                gAACTTGTTGGCACCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTAccc                                          >  6:128127/1‑148 (MQ=255)
                                                                                                      gttgGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATgg                                    <  3:14757/148‑1 (MQ=255)
                                                                                                       ttgGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGc                                  >  5:231629/1‑149 (MQ=255)
                                                                                                            gcCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTg                              >  3:160792/1‑148 (MQ=255)
                                                                                                            gcCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTg                              <  4:160792/148‑1 (MQ=255)
                                                                                                             ccATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCCCCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGc                            >  4:258118/1‑149 (MQ=255)
                                                                                                             ccATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGc                            <  2:218891/149‑1 (MQ=255)
                                                                                                             ccATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGc                            <  5:270053/149‑1 (MQ=255)
                                                                                                                tCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTc                         <  3:210693/149‑1 (MQ=255)
                                                                                                                tCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTc                         <  5:196337/149‑1 (MQ=255)
                                                                                                                tCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTc                         <  8:180566/149‑1 (MQ=255)
                                                                                                                  gTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCa                       <  3:281151/149‑1 (MQ=255)
                                                                                                                              gATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCg           <  7:121386/149‑1 (MQ=255)
                                                                                                                              gATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCg           >  2:190040/1‑149 (MQ=255)
                                                                                                                               aTGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCgg          <  5:319734/149‑1 (MQ=255)
                                                                                                                                tGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCgg          <  5:128127/148‑1 (MQ=255)
                                                                                                                                  gTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGgc        <  1:112170/148‑1 (MQ=255)
                                                                                                                                   tGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGgcgc      <  6:47425/149‑1 (MQ=255)
                                                                                                                                     ccGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAg    <  2:209067/149‑1 (MQ=255)
                                                                                                                                       ggCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGc  >  3:286625/1‑149 (MQ=255)
                                                                                                                                       gatgcCGAAATCGGTGCGCGCACCAATATCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGc  >  8:131525/4‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                         
GAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGC  >  NC_002947/6173374‑6173656

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: