Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,474 T→A 5.5% I159F (ATC→TTC)  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,4740TA 5.5% 220.5 / 2.6 72I159F (ATC→TTC) galUUTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (37/31);  new base A (2/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

TGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTGAGAATGGCGTG  >  NC_002947/4349327‑4349621
                                                                                                                                                   |                                                                                                                                                   
tGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATg                                                                                                                                                    <  6:91834/149‑1 (MQ=255)
  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                                                                                                                                           >  5:338326/1‑141 (MQ=255)
                 gCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATc                                                                                                                                           <  6:338326/141‑1 (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                                                                                                       <  5:175223/124‑1 (MQ=255)
                                                                      acGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGa                                                                                                       >  6:175223/1‑124 (MQ=255)
                                                                        ggAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCt                                                                                                         <  2:30582/120‑1 (MQ=255)
                                                                        ggAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCt                                                                                                         >  1:30582/1‑120 (MQ=255)
                                                                        ggAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGtt                                                                             >  6:106047/1‑148 (MQ=255)
                                                                         gAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGa                                                                           >  5:114403/1‑149 (MQ=255)
                                                                         gAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGa                                                                           <  3:41057/149‑1 (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                                                                                                                                               >  4:414969/1‑44 (MQ=255)
                                                                                                              ccATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGcgac                                                                                                                                               <  3:414969/44‑1 (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)
                                                                                                                                                   |                                                                                                                                                   
TGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTGAGAATGGCGTG  >  NC_002947/4349327‑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: ‑