Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,252,814 T→A 8.2% intergenic (+70/‑49) leuD → / → PP_1987 3‑isopropylmalate dehydratase small subunit/UbiE/COQ5 family methlytransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,252,8140TA 8.2% 202.1 / 4.4 70intergenic (+70/‑49)leuD/PP_19873‑isopropylmalate dehydratase small subunit/UbiE/COQ5 family methlytransferase
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base A (3/3);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCAGTGAAT  >  NC_002947/2252668‑2252959
                                                                                                                                                  |                                                                                                                                                   
tGCAGGACAGCGACGCGATCAAGGCCTACGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGAAGCTGACACTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTAc                                                                                                                                                  >  4:6238/1‑149 (MQ=255)
 gCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTAcc                                                                                                                                                 >  4:194017/1‑149 (MQ=255)
   aGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGCGTATACCcg                                                                                                                                               >  4:251928/1‑149 (MQ=255)
       cAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGaa                                                                                                                                          <  1:284400/149‑1 (MQ=255)
       cAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTAAACCcgcgc                                                                                                                                           <  8:99734/148‑2 (MQ=255)
         gCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGaa                                                                                                                                          >  6:145089/1‑147 (MQ=255)
          cGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAagag                                                                                                                                       <  2:261730/149‑1 (MQ=255)
           gACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAgg                                                                                                                                      <  8:177296/149‑1 (MQ=255)
              gcgATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                                   >  3:101786/1‑149 (MQ=255)
                gATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                                  >  2:290073/1‑148 (MQ=255)
                gATCAAGGCCTGCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGACGCCTCGCTCCTACAGAAAGGTACGCGGTCGCGGTAGGCGCGGTTTTACCCGCGAAGAGGCCGGc                                                                                                                                 >  2:63702/1‑149 (MQ=255)
                          ttCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGcaac                                                                                                                        >  3:125069/1‑148 (MQ=255)
                             gAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACaaa                                                                                                                     >  8:112565/1‑148 (MQ=255)
                              aaGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACaaaaa                                                                                                                   >  2:25097/1‑149 (MQ=255)
                                gCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTAAACCcgagccta                                                                                                                                         >  1:281288/1‑120 (MQ=255)
                                gCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTAAACCcgagccta                                                                                                                                         <  2:281288/126‑7 (MQ=255)
                                   aaGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTaa                                                                                                              <  4:128785/149‑1 (MQ=255)
                                    aGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTaaa                                                                                                             <  4:191206/149‑1 (MQ=255)
                                    aGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTaaa                                                                                                             <  4:101786/149‑1 (MQ=255)
                                      cACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACaaaaa                                                                                                                   >  1:218075/1‑141 (MQ=255)
                                      cACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACaaaaa                                                                                                                   <  2:218075/141‑1 (MQ=255)
                                         cgcgCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAAc                                                                                                        >  3:190091/1‑149 (MQ=255)
                                                  cagccTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCAGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGa                                                                                               >  7:234886/1‑149 (MQ=255)
                                                   agccTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGAc                                                                                              <  1:52676/149‑1 (MQ=255)
                                                      cTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGaa                                                                                                         >  6:356542/1‑135 (MQ=255)
                                                      cTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGaa                                                                                                         <  5:356542/135‑1 (MQ=255)
                                                         ggCTGGTCCGTGATGCCTGAAGTTCGTCGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGAccagca                                                                                         <  8:105542/148‑1 (MQ=255)
                                                          gCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGAccagcacc                                                                                       <  1:25097/149‑1 (MQ=255)
                                                                        ccTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCCGATGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGAt                                                                           <  1:290073/149‑1 (MQ=255)
                                                                         cTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGAtgtg                                                                        <  7:32242/149‑1 (MQ=255)
                                                                           gAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGt                                                                      >  7:324478/1‑149 (MQ=255)
                                                                                      gTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTgg                                                                       <  2:347060/137‑1 (MQ=255)
                                                                                      gTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTgg                                                                       >  1:347060/1‑137 (MQ=255)
                                                                                         aCAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCg                                                        >  7:239395/1‑149 (MQ=255)
                                                                                         aCAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCg                                                        >  5:170754/1‑149 (MQ=255)
                                                                                           aGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCTGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGc                                                      <  3:251928/149‑1 (MQ=255)
                                                                                           aGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGc                                                      <  8:324478/149‑1 (MQ=255)
                                                                                               cGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAg                                                   >  7:22952/1‑148 (MQ=255)
                                                                                                 ggCCTCCTCGCTCCTACAGAAAGGTACCCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAg                                                <  3:6238/149‑1 (MQ=255)
                                                                                                 ggCATCCTCGATCCTACAGAAATTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAg                                                <  8:239395/149‑1 (MQ=255)
                                                                                                   cctcctCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGc                                              <  5:145089/149‑1 (MQ=255)
                                                                                                    ctcctcGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAAc                                                                 >  5:264673/1‑129 (MQ=255)
                                                                                                    ctcctcGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAAc                                                                 <  6:264673/129‑1 (MQ=255)
                                                                                                       ctcGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAgc                                          >  8:283969/1‑149 (MQ=255)
                                                                                                       ctcGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAgc                                          <  7:283969/149‑1 (MQ=255)
                                                                                                              cTACAGACAGTTCCGCGGGCGCTGTAGGCGCGGGTTTACCCGCTAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCTGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTAcc                                   <  8:234886/149‑1 (MQ=39)
                                                                                                              cTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTAc                                    <  4:125069/148‑1 (MQ=255)
                                                                                                               tACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCt                                  >  5:43063/1‑149 (MQ=255)
                                                                                                                  aGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGt                                                                      >  8:248520/1‑110 (MQ=255)
                                                                                                                  aGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGt                                                                      <  7:248520/110‑1 (MQ=255)
                                                                                                                   gAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCt                                  >  3:326096/1‑145 (MQ=255)
                                                                                                                   gAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCt                                  <  4:326096/145‑1 (MQ=255)
                                                                                                                   gAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTcagc                              <  7:112565/149‑1 (MQ=255)
                                                                                                                   gAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTcag                               >  4:275517/1‑148 (MQ=255)
                                                                                                                    aaaGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTcag                               >  1:163954/1‑147 (MQ=255)
                                                                                                                        ttACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAgcgc                         >  3:285245/1‑149 (MQ=255)
                                                                                                                            gcgGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTg                     <  4:190091/149‑1 (MQ=255)
                                                                                                                               gTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGc                                                      >  6:99091/1‑113 (MQ=255)
                                                                                                                               gTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGc                                                      <  5:99091/113‑1 (MQ=255)
                                                                                                                                   cTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGcc               >  3:333064/1‑148 (MQ=255)
                                                                                                                                      tcttcgcgGGTAAACCCGCGCAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGCCCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAg           >  6:224041/5‑149 (MQ=255)
                                                                                                                                        ttcgcgGGTAAACCCGCGAAGAGGCCGGCAGAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAgg          >  1:33091/3‑148 (MQ=255)
                                                                                                                                         gcgcGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGc                                                               <  2:85345/94‑1 (MQ=255)
                                                                                                                                         gcgcGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGc                                                               >  1:85345/1‑94 (MQ=255)
                                                                                                                                         gcgcGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCa        <  6:43063/149‑1 (MQ=255)
                                                                                                                                         gcgcGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCa        <  4:333064/149‑1 (MQ=255)
                                                                                                                                           gcgGGTAAACCCGCGCAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCAGt      <  5:224041/149‑1 (MQ=255)
                                                                                                                                            cgGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCAGTg     <  4:285245/149‑1 (MQ=255)
                                                                                                                                             gggTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCAGTg     >  5:53386/1‑148 (MQ=255)
                                                                                                                                               gTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCAAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCAGTGAAt  <  1:63702/149‑1 (MQ=255)
                                                                                                                                                 ttACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGAtgt                                                                         <  5:187893/76‑1 (MQ=255)
                                                                                                                                                 ttACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGAtgt                                                                         >  6:187893/1‑76 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGCAGGACAGCGACGCGATCAAGGCCTTCGAAGCCAAGCACCGCGCCAGCCAGCCTTGGCTGTTCCGTGATGCCTGAAGTTCGTAGGTGACAGTGCGGGCCTCCTCGCTCCTACAGAAAGTTACGCGGTCGCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACAAAAATCTAAAGGAACGAGACATGACCAGCACCCAGCACATCGATGTGGTCCAACGCCAGTTCGGCGAGCAGGCCAGCGCCTACCTCAGCAGCGCCGTGCACGCCCAAGGCAGTGAAT  >  NC_002947/2252668‑2252959

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: