Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,189,041 A→G 6.6% intergenic (+25/+17) PP_1038 → / ← tatC‑I hypothetical protein/Sec‑independent protein translocase protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,189,0410AG 6.6% 262.6 / 6.8 91intergenic (+25/+17)PP_1038/tatC‑Ihypothetical protein/Sec‑independent protein translocase protein
Reads supporting (aligned to +/- strand):  ref base A (43/42);  new base G (4/2);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

GAAGAGTGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCCGGCG  >  NC_002947/1188894‑1189182
                                                                                                                                                   |                                                                                                                                             
gAAGAGTGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCAc                                                                                                                                              >  2:1357919/1‑149 (MQ=255)
gAAGAGTGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCAc                                                                                                                                              <  1:1357919/149‑1 (MQ=255)
   gagTGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGccc                                                                                                                                          >  2:435297/1‑150 (MQ=255)
        gAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaa                                                                                                                                     <  2:2385338/150‑1 (MQ=255)
        gAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGaa                                                                                                                                     >  1:1036641/1‑150 (MQ=255)
             tGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGc                                                                                                                                <  1:465691/150‑1 (MQ=255)
             tGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGc                                                                                                                                <  1:906996/150‑1 (MQ=255)
                 gCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTa                                                                                                                            <  2:614527/150‑1 (MQ=255)
                   tggtggtggCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGcc                                                                                                                               <  2:359754/145‑1 (MQ=255)
                   tggtggtggCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGcc                                                                                                                               >  1:359754/1‑145 (MQ=255)
                   tggtggtggCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACa                                                                                                                          >  1:732537/1‑150 (MQ=255)
                        gtggCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTCTGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCTGGCACGCCCGCGAAAGGGCCCTACAATCCa                                                                                                                     <  1:1999585/150‑1 (MQ=255)
                        gtggCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCa                                                                                                                     <  2:223838/150‑1 (MQ=255)
                           gCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGACCGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGcc                                                                                                                  <  1:1502382/150‑1 (MQ=255)
                             tGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCcac                                                                                                                >  2:919350/1‑150 (MQ=255)
                                     tGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCt                                                                                                        >  1:1474858/1‑150 (MQ=255)
                                      gATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCtt                                                                                                       >  2:2269070/1‑150 (MQ=255)
                                        ttGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGc                                                                                                                                <  1:1403216/123‑1 (MQ=255)
                                        ttGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGc                                                                                                                                >  2:1403216/1‑123 (MQ=255)
                                          gtgCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCt                                                                                                   <  2:1474858/150‑1 (MQ=255)
                                            gCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGcc                                                                                                                               >  1:138823/1‑120 (MQ=255)
                                            gCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGcc                                                                                                                               <  2:138823/120‑1 (MQ=255)
                                            gCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTgg                                                                                                 <  1:919350/150‑1 (MQ=255)
                                             ctggTCGTCGCTGCGCGAAAGCAGGGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCGCGCCGAGTTCCAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGc                                                                                                <  2:1497462/148‑1 (MQ=255)
                                                gtcgtcGCTGCGCTCAAGCAGAGCCACCCTTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTcc                                                                                             >  1:2060632/1‑150 (MQ=255)
                                                  cgtcgCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATc                                                                                                                       >  1:1658221/1‑122 (MQ=255)
                                                  cgtcgCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATc                                                                                                                       <  2:1658221/122‑1 (MQ=255)
                                                      gCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGc                                                                                       <  2:2382395/150‑1 (MQ=255)
                                                        tGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCtt                                                                                     >  1:1711783/1‑150 (MQ=255)
                                                         gcgcTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCCTc                                                                                    <  1:2402326/150‑1 (MQ=255)
                                                           gcTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCtt                                                                                                       >  1:177003/1‑129 (MQ=255)
                                                           gcTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCtt                                                                                                       <  2:177003/129‑1 (MQ=255)
                                                             tCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACa                                                                                >  2:1756554/1‑150 (MQ=255)
                                                               aaGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGaaa                                                                                                                                    >  2:1055183/1‑96 (MQ=255)
                                                               aaGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGaaa                                                                                                                                    <  1:1055183/96‑1 (MQ=255)
                                                               aaGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACAtt                                                                              >  2:650163/1‑150 (MQ=255)
                                                                  cAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCca                                                                           <  2:1564459/150‑1 (MQ=255)
                                                                     agCCACCCGTACGAGACGCCAGCCTATGACGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCCTGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACaa                                                                        <  1:770367/150‑1 (MQ=255)
                                                                     agCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCcaaag                                                                        <  1:2269070/150‑4 (MQ=255)
                                                                        cACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGcc                                                                     >  1:1757784/1‑150 (MQ=255)
                                                                         aCCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCa                                                                    >  2:26312/1‑150 (MQ=255)
                                                                           ccGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGc                                                                                                                                <  1:2226128/88‑1 (MQ=255)
                                                                           ccGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGc                                                                                                                                >  2:2226128/1‑88 (MQ=255)
                                                                           ccGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCaa                                                                                  >  1:2140924/1‑134 (MQ=255)
                                                                           ccGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCaa                                                                                  <  2:2140924/134‑1 (MQ=255)
                                                                             gTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACa                                                                                >  1:1244995/1‑134 (MQ=255)
                                                                             gTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACa                                                                                <  2:1244995/134‑1 (MQ=255)
                                                                                      gCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCAtctc                                                       >  2:1810041/1‑150 (MQ=255)
                                                                                         aGCCTATGAATTCTGACGGCTCGCCGACTTCTAGAACGCATCACGACTATTCGCGGGCACGCCCGCGAAAGGGCCCTACATTCCAGCAACAATCACCTTGCCTGGCTCCTGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa                                                    <  1:553972/150‑1 (MQ=255)
                                                                                         aGCCTATGAAGTCTGGCGGCTCGTTAAGTGCTAGAACGCATCACGTCTATGCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa                                                    >  1:1261321/1‑150 (MQ=255)
                                                                                         aGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGTGGGCGTGCCCGCGaa                                                                                                                                     >  2:2312914/1‑69 (MQ=255)
                                                                                         aGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGTGGGCGTGCCCGCGaa                                                                                                                                     <  1:2312914/69‑1 (MQ=255)
                                                                                             tATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGc                                                >  2:1228322/1‑150 (MQ=255)
                                                                                               tGAAGTCTGGCGGCTAGCCGAGTGCGAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAc                                              <  2:1757784/150‑1 (MQ=255)
                                                                                                 aaGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAccc                                            >  1:1506614/1‑150 (MQ=255)
                                                                                                  aGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACaa                                                                        >  1:1556533/1‑121 (MQ=255)
                                                                                                  aGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACaa                                                                        <  2:1556533/121‑1 (MQ=255)
                                                                                                   gTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGAGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCcac                                          >  1:1111121/1‑150 (MQ=255)
                                                                                                        gcggcTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg                                     >  1:1447132/1‑150 (MQ=255)
                                                                                                         cggcTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGc                                    <  1:650163/150‑1 (MQ=255)
                                                                                                          ggcTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTc                                                                             <  1:603935/108‑1 (MQ=255)
                                                                                                          ggcTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTc                                                                             >  2:603935/1‑108 (MQ=255)
                                                                                                          ggcTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCa                                   >  2:786824/1‑150 (MQ=255)
                                                                                                             tCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACa                                                  >  2:1188656/1‑132 (MQ=255)
                                                                                                             tCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACa                                                  <  1:1188656/132‑1 (MQ=255)
                                                                                                               gCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGc                               <  1:2057047/149‑1 (MQ=255)
                                                                                                               gCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGc                               >  2:2057047/1‑149 (MQ=255)
                                                                                                                ccGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCa                             <  2:2060632/150‑1 (MQ=255)
                                                                                                                  gAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGc                           <  2:1447132/150‑1 (MQ=255)
                                                                                                                     ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAAcc                                                              <  1:1697648/112‑1 (MQ=255)
                                                                                                                     ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAAca                                                              >  2:1697648/1‑111 (MQ=255)
                                                                                                                     ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGAGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATc                        <  2:1111121/150‑1 (MQ=255)
                                                                                                                        tAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCAAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTc                     >  1:1819050/1‑150 (MQ=255)
                                                                                                                           aaCGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATc                                       >  2:577849/1‑129 (MQ=255)
                                                                                                                           aaCGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATc                                       <  1:577849/129‑1 (MQ=255)
                                                                                                                            aCGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCg                 >  2:2294400/1‑150 (MQ=255)
                                                                                                                             cGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAg                                                                   <  2:958338/99‑1 (MQ=255)
                                                                                                                             cGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAg                                                                   >  1:958338/1‑99 (MQ=255)
                                                                                                                                aTCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGaaaa             <  1:435297/150‑1 (MQ=255)
                                                                                                                                aTCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGaaaa             >  1:1257998/1‑150 (MQ=255)
                                                                                                                                 tCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAGTCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAac            >  1:556357/1‑150 (MQ=255)
                                                                                                                                 tCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAac            <  2:1711783/150‑1 (MQ=255)
                                                                                                                                   aCGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCTAg                            <  1:2407545/132‑1 (MQ=255)
                                                                                                                                   aCGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAg                            >  2:2407545/1‑132 (MQ=255)
                                                                                                                                   aCGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAacac          >  2:1358144/1‑150 (MQ=255)
                                                                                                                                       cctttCGCGGGCGTGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATcc      >  2:317048/4‑150 (MQ=255)
                                                                                                                                       cTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCGCAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATATCGAACAGCACCCACATCGGCACTGCCAGCATCCTCTGCGAAAACACATcc      >  1:28665/1‑150 (MQ=255)
                                                                                                                                         tttCGCGGGCGTGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCcgg    >  2:2051980/2‑150 (MQ=255)
                                                                                                                                         tttCGCGGGCGTGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGATTAAATATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCcgg    >  1:1145029/2‑150 (MQ=255)
                                                                                                                                         aTTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCACCCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCcgg    <  2:1036641/150‑1 (MQ=255)
                                                                                                                                           tCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCcggcc  <  1:1756554/150‑2 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GAAGAGTGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAAAGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTCTGCGAAAACACATCCGGCG  >  NC_002947/1188894‑1189182

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