Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,476 G→T 5.6% A158D (GCC→GAC)  galU ← UTP‑glucose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,349,4760GT 5.6% 231.9 / 2.9 71A158D (GCC→GAC) galUUTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (37/30);  new base T (2/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTGAGAATGGCGTG  >  NC_002947/4349329‑4349621
                                                                                                                                                   |                                                                                                                                                 
tAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTCGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGc                                                                                                                                                  >  6:210123/1‑149 (MQ=255)
  aCGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGTGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGcga                                                                                                                                                <  7:6655/149‑1 (MQ=255)
  aCGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGCCATCCAGGAAGTcga                                                                                                                                                >  5:362728/1‑149 (MQ=255)
    ggCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGCCATCCAGGAAGTcgacc                                                                                                                                              >  6:404567/1‑148 (MQ=255)
        gatgatCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGCCATCCAGGAAGTcgacccgca                                                                                                                                          <  5:404567/149‑6 (MQ=255)
        gatgatCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATc                                                                                                                                           >  6:83953/1‑148 (MQ=255)
         atgatCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGCCATCCAGGAAGTcgacccgcag                                                                                                                                         <  6:362728/149‑7 (MQ=255)
         atgatCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGa                                                                                                                                         <  7:329399/149‑1 (MQ=255)
               gCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATc                                                                                                                                           <  6:338326/141‑1 (MQ=255)
               gCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATc                                                                                                                                           >  5:338326/1‑141 (MQ=255)
                         aaGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTg                                                                                                                         <  1:9085/149‑1 (MQ=255)
                           ggCGCGTCTTCCGGGGCCGGTTTTTGGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTa                                                                                                                       <  1:426662/149‑1 (MQ=255)
                            gcgcGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTAc                                                                                                                      <  2:358247/149‑1 (MQ=255)
                                    tCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGAc                                                                                                              >  5:379730/1‑149 (MQ=255)
                                     ccGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGttt                                                                                                                          >  3:129642/1‑136 (MQ=255)
                                     ccGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGttt                                                                                                                          <  4:129642/136‑1 (MQ=255)
                                     ccGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGAc                                                                                                              >  2:199170/1‑148 (MQ=255)
                                       ggggCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCAt                                                                                                           <  2:122639/149‑1 (MQ=255)
                                           ccGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGa                                                                                                       <  1:199170/149‑1 (MQ=255)
                                              gTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATc                                                                                                                                           <  1:98168/110‑1 (MQ=255)
                                              gTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGATGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATc                                                                                                                                           >  2:98168/1‑110 (MQ=255)
                                                     gACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGc                                                                                             >  5:181545/1‑149 (MQ=255)
                                                      aCCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGcc                                                                                            >  2:280302/1‑149 (MQ=255)
                                                                 cacaCGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCg                                                                                 >  3:425829/1‑149 (MQ=255)
                                                                  acacGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCgg                                                                                <  4:403800/149‑1 (MQ=255)
                                                                    acGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGa                                                                                                       >  6:175223/1‑124 (MQ=255)
                                                                    acGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGa                                                                                                       <  5:175223/124‑1 (MQ=255)
                                                                      ggAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCt                                                                                                         <  2:30582/120‑1 (MQ=255)
                                                                      ggAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCt                                                                                                         >  1:30582/1‑120 (MQ=255)
                                                                      ggAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGtt                                                                             >  6:106047/1‑148 (MQ=255)
                                                                       gAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGa                                                                           <  3:41057/149‑1 (MQ=255)
                                                                       gAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGa                                                                           >  5:114403/1‑149 (MQ=255)
                                                                        aaGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGAc                                                                          <  6:379730/149‑1 (MQ=255)
                                                                         aGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACg                                                                         <  6:114403/149‑1 (MQ=255)
                                                                           aTGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGca                                                                       >  8:175562/1‑149 (MQ=255)
                                                                             gtcgtcGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGca                                                                       >  2:305086/1‑147 (MQ=255)
                                                                                 tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCgg                                                                                <  4:65816/134‑1 (MQ=255)
                                                                                 tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCgg                                                                                >  3:65816/1‑134 (MQ=255)
                                                                                 tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGgtc                                                                 >  1:304267/1‑149 (MQ=255)
                                                                                 tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGgtc                                                                 <  4:411403/149‑1 (MQ=255)
                                                                                 tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGgtc                                                                 >  3:411403/1‑149 (MQ=255)
                                                                                     ggATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTcc                                                             <  1:280302/149‑1 (MQ=255)
                                                                                       atcatcTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGc                                                           >  5:419121/1‑149 (MQ=255)
                                                                                          atcTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAg                                                        <  2:304267/149‑1 (MQ=255)
                                                                                            ctcACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGca                                                      >  6:46347/1‑149 (MQ=255)
                                                                                              cACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGcac                                                     >  5:40940/1‑148 (MQ=255)
                                                                                               aCCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGcacca                                                   >  3:184307/1‑149 (MQ=255)
                                                                                               aCCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGcac                                                     <  5:106047/147‑1 (MQ=255)
                                                                                                 cGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCaccacc                                                 >  3:77131/1‑149 (MQ=255)
                                                                                                    cGATCACACCATACTTGTTGCTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACcgcg                                              <  2:371286/149‑1 (MQ=255)
                                                                                                         acacCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACgg                                         <  4:77131/149‑1 (MQ=255)
                                                                                                         acacCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACgg                                         <  5:83953/149‑1 (MQ=255)
                                                                                                            ccATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGcgac                                                                                                                                               <  3:414969/44‑1 (MQ=255)
                                                                                                            ccATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGcgac                                                                                                                                               >  4:414969/1‑44 (MQ=255)
                                                                                                            ccATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTc                                      <  5:210123/149‑1 (MQ=255)
                                                                                                                       ggTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGAt                             >  7:13473/1‑147 (MQ=255)
                                                                                                                        gTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAg                          >  4:144392/1‑149 (MQ=255)
                                                                                                                         tttCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGt                         >  6:381266/1‑149 (MQ=255)
                                                                                                                            ccTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTggg                      >  5:129566/1‑149 (MQ=255)
                                                                                                                             cTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGggc                     >  6:15259/1‑149 (MQ=255)
                                                                                                                                cGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGcc                                                          <  8:389949/109‑1 (MQ=255)
                                                                                                                                cGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGcc                                                          >  7:389949/1‑109 (MQ=255)
                                                                                                                                 gggTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGcc                 <  7:175562/149‑1 (MQ=255)
                                                                                                                                  ggTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGggcggc                  >  2:15760/1‑147 (MQ=255)
                                                                                                                                      gACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTga            <  5:46347/149‑1 (MQ=255)
                                                                                                                                       aCTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGATTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGCCTCGTCACCGATCAGTGGGCGGCCAGt              >  3:318024/1‑146 (MQ=255)
                                                                                                                                         ttgcTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCAGCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACAGAGAAGGGCTGGTGACCGATAAGTGGGCGGCCAGGGCGaa         >  4:348490/4‑149 (MQ=255)
                                                                                                                                         ttCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCAGCCAGTGAGaa         >  2:425872/1‑149 (MQ=255)
                                                                                                                                         ttCCTGGATGGCGACGATCGAGCAGCGGTACTGATTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTgaga          >  6:351592/1‑148 (MQ=255)
                                                                                                                                              ggATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTGAGAATGGCg    >  4:409481/1‑149 (MQ=255)
                                                                                                                                                aTGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTGAGAATGGCGTg  <  1:15760/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTGAGAATGGCGTG  >  NC_002947/4349329‑4349621

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: ‑