Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,173,535 T→A 5.2% E368D (GAA→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,5350TA 5.2% 240.9 / 2.6 77E368D (GAA→GATglmUbifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base A (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

TCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAACGGGCCAGCATC  >  NC_002947/6173390‑6173675
                                                                                                                                                  |                                                                                                                                            
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:24039/149‑1 (MQ=255)
                      aGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGt                                                                                                                     <  7:174777/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                                                                                                                                         >  6:315387/1‑118 (MQ=39)
                                 gCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGAAATCGGTGCGCGCACCAATATCGGc                                                                                                                                         <  5:315387/118‑1 (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                                                                                                                         >  8:75813/1‑119 (MQ=255)
                                                cGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTc                                                                                                                         <  7:75813/119‑1 (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                                                 <  4:160792/148‑1 (MQ=255)
                                                                                            gcCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTg                                                 >  3:160792/1‑148 (MQ=255)
                                                                                             ccATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCCCCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGc                                               >  4:258118/1‑149 (MQ=255)
                                                                                             ccATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGc                                               <  5:270053/149‑1 (MQ=255)
                                                                                             ccATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGc                                               <  2:218891/149‑1 (MQ=255)
                                                                                                tCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTc                                            <  8:180566/149‑1 (MQ=255)
                                                                                                tCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTc                                            <  5:196337/149‑1 (MQ=255)
                                                                                                tCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTc                                            <  3:210693/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)
                                                                                                                                   ggTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAACggg          <  1:130189/148‑1 (MQ=255)
                                                                                                                                    gTGCGCGCACAGATTTCGGCAGCCCCCAGGTATGTCAGGTGACCAGCCTTGTCCCCTTCACCCAGGTTGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAACGGTcc        <  8:305232/149‑1 (MQ=255)
                                                                                                                                     tGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAACGGGc         >  5:609/1‑147 (MQ=255)
                                                                                                                                      gcgcgcACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTg                                                 >  4:100212/1‑106 (MQ=255)
                                                                                                                                      gcgcgcACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTg                                                 <  3:100212/106‑1 (MQ=255)
                                                                                                                                          gcACCGATTTTGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAACGGGCCAGCATc  >  3:259962/1‑149 (MQ=255)
                                                                                                                                          gcACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAACGGGCCAGCATc  >  2:45743/1‑149 (MQ=255)
                                                                                                                                            aCCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAACgg           >  8:187708/1‑138 (MQ=255)
                                                                                                                                            aCCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAACgg           <  7:187708/138‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
TCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACCAGCCTTGGCCCCTTCACCCAGGTGGGCGTTTTTCAGTTCAACGAAGTTACCCACATGGGCCTTGGCGTCCAGTACGCTGCCCGGGCGCAGGCGAGCAAACGGGCCAGCATC  >  NC_002947/6173390‑6173675

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: