Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,143,926 T→A 10.6% intergenic (+35/‑88) PP_3647 → / → PP_3648 oxidoreductase/carboxymuconolactone decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,143,9260TA10.6% 192.2 / 12.9 66intergenic (+35/‑88)PP_3647/PP_3648oxidoreductase/carboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base A (3/4);  total (30/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 = 7.68e-01

CACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTA  >  NC_002947/4143784‑4144075
                                                                                                                                              |                                                                                                                                                     
cACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCcgc                                                                                                                                                >  2:1263683/1‑150 (MQ=255)
  cgcgATCATTGTCGGCGCGGTGAAGCACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCATGCTCGGCGCCCTGCCCCAGGTCAGCAAAGACGTCTTGACTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGa                                                                                                                                              <  1:283691/150‑1 (MQ=255)
   gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                             <  1:393625/150‑1 (MQ=255)
     gATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCGTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGa                                                                                                                                           >  1:1648342/1‑150 (MQ=255)
     gATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGa                                                                                                                                           <  1:641512/150‑1 (MQ=255)
      aTCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGATGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGAt                                                                                                                                          >  1:803072/1‑150 (MQ=255)
           tGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGcccc                                                                                                                                     >  1:547383/1‑150 (MQ=255)
             tCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATAACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTg                                                                                                                                   <  2:134983/150‑1 (MQ=255)
               cgctcgGTGAAGAACCTGGTGTAGCCTCAGCAGCCGATGACCCCGATGGTGTAACCCAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCtttt                                                                                                                                 <  1:329977/146‑1 (MQ=255)
                   cgGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCa                                                                                                                             <  2:1153145/150‑1 (MQ=255)
                    ggTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAg                                                                                                                            <  2:1367905/150‑1 (MQ=255)
                        aagaaCCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGt                                                                                                                        >  2:647134/1‑150 (MQ=255)
                        aagaaCCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGa                                                                                                                                           <  2:31721/131‑1 (MQ=255)
                        aagaaCCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGa                                                                                                                                           >  1:31721/1‑131 (MQ=255)
                         agaaCCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTc                                                                                                                       >  2:587889/1‑150 (MQ=255)
                             ccTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCgcgc                                                                                                                   >  2:460185/1‑150 (MQ=255)
                             ccTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGAGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCgcgc                                                                                                                   <  1:1436613/150‑1 (MQ=255)
                                    ttGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGa                                                                                                                                           >  1:698756/1‑119 (MQ=255)
                                    ttGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGa                                                                                                                                           <  2:698756/119‑1 (MQ=255)
                                          gAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCAcc                                                                                                      >  1:503967/1‑150 (MQ=255)
                                             cagcCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTc                                                                                                   >  1:303797/1‑150 (MQ=255)
                                              agacGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGc                                                                                                                                        >  1:762455/4‑112 (MQ=255)
                                              agacGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGc                                                                                                                                        <  2:762455/109‑1 (MQ=255)
                                                       ccccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGcc                                                                                         >  1:34366/1‑150 (MQ=255)
                                                            aTGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttt                                                                                    <  2:1081385/150‑1 (MQ=255)
                                                                tgtACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTa                                                                                >  2:1310035/1‑150 (MQ=255)
                                                                  tACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATg                                                                              >  1:860035/1‑150 (MQ=255)
                                                                   aCCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGACATCTTCGCGGGCTTGCCCGCCAAGCTGCCCCTCTTTTCAGGGGTCGCGCCACTCGTCCCACCCCCCGCTGCGGCCCCTTTTTTATGc                                                                             >  1:356446/1‑150 (MQ=255)
                                                                     cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCg                                                                           <  2:860035/150‑1 (MQ=255)
                                                                     cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCg                                                                           <  2:547383/150‑1 (MQ=255)
                                                                       caAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                         <  2:1648342/150‑1 (MQ=255)
                                                                        aaGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCAcc                                                                                                      <  2:1580064/120‑1 (MQ=255)
                                                                        aaGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCAcc                                                                                                      >  1:1580064/1‑120 (MQ=255)
                                                                             cAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTg                                                                   <  1:460185/150‑1 (MQ=255)
                                                                              aaGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGc                                                                  <  2:803072/150‑1 (MQ=255)
                                                                                  tCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCACCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCTCCTTTTTTATGCCGGAGTACTGCACAt                                                              <  1:221574/150‑1 (MQ=255)
                                                                                       gcCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTAt                                                                               >  1:1280565/1‑128 (MQ=255)
                                                                                       gcCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTAt                                                                               <  2:1280565/128‑1 (MQ=255)
                                                                                       gcCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTa                                                                      <  2:814695/137‑1 (MQ=255)
                                                                                       gcCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTa                                                                      >  1:814695/1‑137 (MQ=255)
                                                                                        cccTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCaa                                                        <  1:647134/150‑1 (MQ=255)
                                                                                          cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGt                                                                                                                        >  2:1289786/1‑84 (MQ=255)
                                                                                          cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGt                                                                                                                        <  1:1289786/84‑1 (MQ=255)
                                                                                          cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATg                                                      <  2:34366/150‑1 (MQ=255)
                                                                                            gccgccACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGaa                                                    <  1:587889/150‑1 (MQ=255)
                                                                                                  aCTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTAc                                              <  1:1263683/150‑1 (MQ=255)
                                                                                                       gCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGaaaa                                         >  2:1147872/1‑150 (MQ=255)
                                                                                                         tGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAgg                                       >  2:138460/1‑150 (MQ=255)
                                                                                                           aTACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGc                                     <  2:1445092/150‑1 (MQ=255)
                                                                                                             aCATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGaaa                                                   >  1:114144/1‑134 (MQ=255)
                                                                                                             aCATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGaaa                                                   <  2:114144/134‑1 (MQ=255)
                                                                                                             aCATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCt                                   <  1:138460/150‑1 (MQ=255)
                                                                                                                 cTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGcc                                                                            <  2:622399/105‑1 (MQ=255)
                                                                                                                 cTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGcc                                                                            >  1:622399/1‑105 (MQ=255)
                                                                                                                         gATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCAc                       <  2:303797/150‑1 (MQ=255)
                                                                                                                                  aTCTTCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACt              >  1:1293237/1‑150 (MQ=255)
                                                                                                                                     ttCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg           >  2:1536889/1‑150 (MQ=255)
                                                                                                                                     ttCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATg                                                                              >  2:882839/1‑83 (MQ=255)
                                                                                                                                     ttCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATg                                                                              <  1:882839/83‑1 (MQ=255)
                                                                                                                                      tCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGcc                                                                            >  2:786721/1‑84 (MQ=255)
                                                                                                                                      tCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGcc                                                                            <  1:786721/84‑1 (MQ=255)
                                                                                                                                        gcgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGa        >  2:1044307/1‑150 (MQ=255)
                                                                                                                                         cgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAg       <  1:1310035/150‑1 (MQ=255)
                                                                                                                                          gggCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg           <  2:26763/145‑1 (MQ=255)
                                                                                                                                          gggCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg           >  1:26763/1‑145 (MQ=255)
                                                                                                                                            gCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGGCGCGCTAATCGGCTCAGCCTCGGTTGCGGCCCCTTTGTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGAc    <  2:356446/150‑1 (MQ=255)
                                                                                                                                              ttGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAAGCGTCTCAGCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTa  >  1:698150/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
CACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTA  >  NC_002947/4143784‑4144075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: