Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,470 T→A 5.6% Q160L (CAG→CTG)  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,4700TA 5.6% 215.9 / 2.2 71Q160L (CAG→CTG) galUUTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (37/30);  new base A (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

AGGATGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTGAGAA  >  NC_002947/4349323‑4349614
                                                                                                                                                   |                                                                                                                                                
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                                                                                                                                    >  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                                                                                                                   <  4:129642/136‑1 (MQ=255)
                                           ccGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGttt                                                                                                                   >  3:129642/1‑136 (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                                                                                                  >  1:30582/1‑120 (MQ=255)
                                                                            ggAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCt                                                                                                  <  2:30582/120‑1 (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                                                          >  3:411403/1‑149 (MQ=255)
                                                                                       tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGgtc                                                          >  1:304267/1‑149 (MQ=255)
                                                                                       tcgCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGgtc                                                          <  4:411403/149‑1 (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                                                   <  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)
                                                                                                                                                   |                                                                                                                                                
AGGATGTAACGGCCGATGATCGCCAGGTTCGAAGGCGCGTCTTCCGGGGCCGGTTTTTCGACCATGTTGGTCACACGGAAGATGTCGTCGCGGATCATCTCACCGGCGATCACACCATACTTGTTGGTTTCCTGCGGGTCGACTTCCTGGATGGCGACGATCGAGCAGCGGTACTGTTTGTACAGCTTGACCATCTGAGTCAGCACGCCGTCACCTTCCGGGTTGACGCACAGGTCGTCCGCCAGCACCACCGCGAACGGCTCGTCACCGATCAGTGGGCGGCCAGTGAGAA  >  NC_002947/4349323‑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: ‑