Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,717,954 G→T 6.8% E173D (GAG→GAT)  fabB → 3‑oxoacyl‑ACP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,717,9540GT 6.8% 234.4 / 6.5 73E173D (GAG→GAT) fabB3‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base T (3/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-01

GGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCGCGCCTAC  >  NC_002947/4717817‑4718095
                                                                                                                                         |                                                                                                                                             
ggCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGt                                                                                                                                     <  3:206954/148‑1 (MQ=255)
   ccATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGATCTGCTCCAGt                                                                                                                                     >  5:196840/1‑145 (MQ=255)
   ccATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGATCTGCTCCAGt                                                                                                                                     <  6:196840/145‑1 (MQ=255)
   ccATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAg                                                                                                                                            >  5:380474/1‑138 (MQ=255)
   ccATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAg                                                                                                                                            <  6:380474/138‑1 (MQ=255)
       aCCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAg                                                                                                                             >  1:293974/1‑149 (MQ=255)
           cgTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGc                                                                                                                                              <  6:328249/128‑1 (MQ=255)
           cgTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGc                                                                                                                                              >  5:328249/1‑128 (MQ=255)
                 ccGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCg                                                                                                                   >  3:18718/1‑149 (MQ=255)
                   gCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTc                                                                                                                 <  4:18718/149‑1 (MQ=255)
                       cATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCg                                                                                                             >  3:302160/1‑149 (MQ=255)
                       cATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGACTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCg                                                                                                             >  5:207866/1‑149 (MQ=255)
                        aTGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGc                                                                                                            <  4:302160/149‑1 (MQ=255)
                         tGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGcc                                                                                                           >  7:167293/1‑149 (MQ=255)
                              agcaCCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCg                                                                                                          <  5:26305/145‑1 (MQ=255)
                              agcaCCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCg                                                                                                          >  6:26305/1‑145 (MQ=255)
                                 aCCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCgg                                                                                                   <  2:147539/149‑1 (MQ=255)
                                       tCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTgaaga                                                                                             <  3:386437/149‑1 (MQ=255)
                                         ggCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTgaag                                                                                              >  8:300398/1‑146 (MQ=255)
                                              gCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCAc                                                                                      >  2:2603/1‑149 (MQ=255)
                                                cTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAagaaga                                                                                          <  8:349/143‑1 (MQ=255)
                                                cTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAagaaga                                                                                          >  7:349/1‑143 (MQ=255)
                                                        cccGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTggg                                                                                                                                  >  3:384017/1‑95 (MQ=255)
                                                        cccGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTggg                                                                                                                                  <  4:384017/95‑1 (MQ=255)
                                                        cccGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTgg                                                                                   >  4:196231/1‑142 (MQ=255)
                                                        cccGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTgg                                                                                   <  3:196231/142‑1 (MQ=255)
                                                        cccGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGt                                                                            >  5:317689/1‑149 (MQ=255)
                                                        cccGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGt                                                                            >  8:35911/1‑149 (MQ=255)
                                                               aaGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGGGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCt                                                                     >  5:129033/1‑149 (MQ=255)
                                                               aaGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCt                                                                     >  6:274861/1‑149 (MQ=255)
                                                                 gATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCa                                                                                       >  6:67918/1‑129 (MQ=255)
                                                                 gATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCa                                                                                       <  5:67918/129‑1 (MQ=255)
                                                                   tCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTgaagaa                                                                                            <  5:15125/122‑1 (MQ=255)
                                                                   tCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTgaagaa                                                                                            >  6:15125/1‑122 (MQ=255)
                                                                     aaGGGCATCAACTACTCGATCTCGTCGGCCAGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTgaagaa                                                                                            >  1:32326/1‑120 (MQ=255)
                                                                     aaGGGCATCAACTACTCGATCTCGTCGGCCAGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTgaagaa                                                                                            <  2:32326/120‑1 (MQ=255)
                                                                      aGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGAt                                                              <  4:187134/149‑1 (MQ=255)
                                                                         gCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGGGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGc                                                            <  6:129033/148‑1 (MQ=255)
                                                                               actactCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGgc                                                     >  7:262687/1‑149 (MQ=255)
                                                                                 tactCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGTACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGgcgc                                                   >  7:214892/1‑149 (MQ=255)
                                                                                   ctCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGa                                                                                  <  4:316794/116‑1 (MQ=255)
                                                                                   ctCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGa                                                                                  >  3:316794/1‑116 (MQ=255)
                                                                                   ctCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGccc                                                 >  7:373230/1‑149 (MQ=255)
                                                                                      gATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCt                                                <  1:2603/147‑1 (MQ=255)
                                                                                      gATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGt                                              >  8:1768/1‑149 (MQ=255)
                                                                                           cgtcgGCCTGCGCCACCTCGGCTCACTGGATCGGTACCGCACTGGAGCAGATCCTGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAg                                                                                 <  3:352905/109‑1 (MQ=255)
                                                                                           cgtcgGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGAc                                          >  5:304510/1‑148 (MQ=255)
                                                                                            gtcgGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCg                                            <  2:293974/145‑1 (MQ=255)
                                                                                            gtcgGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCa                                        <  7:300398/149‑1 (MQ=255)
                                                                                            gtcgGCCTGCGCCACCTCGACTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTc                                             <  6:207866/144‑1 (MQ=255)
                                                                                             tcgGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCaa                                       <  8:167293/149‑1 (MQ=255)
                                                                                             tcgGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCaa                                       <  7:1768/149‑1 (MQ=255)
                                                                                               ggCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAgc                                     <  7:38230/149‑1 (MQ=255)
                                                                                                  cTGCGCCACCTCGGCTCACTGCATCGGTACCGTACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCa                                  <  8:214892/149‑1 (MQ=255)
                                                                                                     cgcCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCaa                                 <  8:262687/147‑1 (MQ=255)
                                                                                                       ccACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGaa                             <  5:274861/149‑1 (MQ=255)
                                                                                                        cACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGaaa                            <  6:317689/149‑1 (MQ=255)
                                                                                                        cACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGaaa                            >  7:80822/1‑149 (MQ=255)
                                                                                                        cACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGaaa                            <  5:206947/149‑1 (MQ=255)
                                                                                                              ggCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCTGTTAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGTGCGCCCTGTCGACCAAGCGCAACGAAACCCCgg                      <  5:65678/149‑1 (MQ=255)
                                                                                                               gCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCgg                      >  7:293735/1‑148 (MQ=255)
                                                                                                                cTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCaa                                       >  4:268341/1‑130 (MQ=255)
                                                                                                                cTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCaa                                       <  3:268341/130‑1 (MQ=255)
                                                                                                                   aCTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTg                 >  6:142884/1‑149 (MQ=255)
                                                                                                                   aCTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTg                 <  8:293735/149‑1 (MQ=255)
                                                                                                                       cATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGcc              <  6:304510/148‑1 (MQ=255)
                                                                                                                                cccACTGGATCTGCTCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTcgcgc        >  6:191965/3‑145 (MQ=255)
                                                                                                                                cccACTGGATCTGCTCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTcgcgc        <  5:191965/143‑1 (MQ=255)
                                                                                                                                 gCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCGCGCCTa   >  7:212963/1‑149 (MQ=255)
                                                                                                                                  cACTGGATCTGCTCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCGCGCCTac  >  3:32757/1‑149 (MQ=255)
                                                                                                                                  cACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCGCGCCTac  >  3:233138/1‑149 (MQ=255)
                                                                                                                                     tGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTcgc          <  5:357547/138‑1 (MQ=255)
                                                                                                                                     tGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTcgc          >  6:357547/1‑138 (MQ=255)
                                                                                                                                         |                                                                                                                                             
GGCCCATACCGCGTTACCCGCACCATGGGCAGCACCGTGTCGGCATGCCTGGCCACCCCGTTCAAGATCAAGGGCATCAACTACTCGATCTCGTCGGCCTGCGCCACCTCGGCTCACTGCATCGGTACCGCACTGGAGCAGATCCAGTGGGGCAAGCAGGACATCGTCTTCGCCGGTGGCGGTGAAGAAGAGCACTGGAGCCAGTCGTTCCTGTTCGATGCCATGGGCGCCCTGTCGACCAAGCGCAACGAAACCCCGGAACTGGCCTCGCGCGCCTAC  >  NC_002947/4717817‑4718095

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