Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,466 T→A 5.6% E161D (GAA→GAT)  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,4660TA 5.6% 220.1 / 2.6 71E161D (GAA→GAT) galUUTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base A (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 = 7.19e-01

GTCAGGATGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTGAGAA  >  NC_002947/4349320‑4349614
                                                                                                                                                  |                                                                                                                                                    
gTCAGGATGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGAGTTc                                                                                                                                                    <  7:384407/149‑1 (MQ=255)
gTCAGGATGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTc                                                                                                                                                    <  4:55777/149‑1 (MQ=255)
   aGGATGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCt                                                                                                                                                  <  5:419291/148‑1 (MQ=255)
     gATGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGa                                                                                                                                               >  3:403800/1‑149 (MQ=255)
       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                                                                                                                                    <  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                                                                    >  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                                                                         >  3:65816/1‑134 (MQ=255)
                                                                                          tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCgg                                                                         <  4:65816/134‑1 (MQ=255)
                                                                                          tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGgtc                                                          <  4:411403/149‑1 (MQ=255)
                                                                                          tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGgtc                                                          >  3:411403/1‑149 (MQ=255)
                                                                                          tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGgtc                                                          >  1:304267/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                                  <  5:83953/149‑1 (MQ=255)
                                                                                                                  acacCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACgg                                  <  4:77131/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                                                   >  7:389949/1‑109 (MQ=255)
                                                                                                                                         cGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGcc                                                   <  8:389949/109‑1 (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)
                                                                                                                                                  ttCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCAGCCAGTGAGaa  >  2:425872/1‑149 (MQ=255)
                                                                                                                                                  ttCCTGGATGGCGACGATCGAGCAGCGGTACTGATTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTgaga   >  6:351592/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTCAGGATGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTGAGAA  >  NC_002947/4349320‑4349614

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