Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,876,231 C→G 5.1% L227V (CTG→GTG)  xcpS → type II secretion pathway protein XcpS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,876,2310CG 5.1% 382.7 / 8.7 117L227V (CTG→GTG) xcpStype II secretion pathway protein XcpS
Reads supporting (aligned to +/- strand):  ref base C (57/54);  new base G (3/3);  total (60/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 117 total aligned reads displayed.

GCTGGTGGGAGGCGGTGTGGTGCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTCTATCGCT  >  NC_002947/3876083‑3876378
                                                                                                                                                    |                                                                                                                                                   
gCTGGTGGGAGGCGGTGTGGTGCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACt                                                                                                                                                    <  2:762047/150‑1 (MQ=255)
gCTGGTGGGAGGCGGTGTGGTGCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACt                                                                                                                                                    <  1:2015938/150‑1 (MQ=255)
  tggtggGAGGCGGTGTGGTGATGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTgg                                                                                                                                                  >  2:2340224/1‑150 (MQ=255)
    gtggGAGGCGGTGTGGTGCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGgcg                                                                                                                                                >  1:2093418/1‑150 (MQ=255)
      gggAGGCGGTGTGGTGCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCt                                                                                                                                              <  2:1651549/150‑1 (MQ=255)
          ggCGGTGTGGTGCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCAGTGGCGCttgct                                                                                                                                          <  2:1821939/150‑5 (MQ=255)
             ggTGTGGTGCTTTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCAGTGGCGCttgctggg                                                                                                                                       <  1:370956/150‑8 (MQ=255)
              gtgtGGTGCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCg                                                                                                                                      >  1:812525/1‑150 (MQ=255)
                   gtgCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGt                                                                                                                                 >  1:58116/1‑150 (MQ=255)
                    tgCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTg                                                                                                                                <  1:1370911/150‑1 (MQ=255)
                     gCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGACCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTgg                                                                                                                               <  2:310628/150‑1 (MQ=255)
                      cTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTggg                                                                                                                              <  2:319841/150‑1 (MQ=255)
                      cTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTggg                                                                                                                              >  1:906898/1‑150 (MQ=255)
                      cTGTTCCTGCTCGGCGATGTGGTGCCGCGTTTAAACAAGATATTGGAGGCCATTTGCACCGAGCTACCCTGGTCATCTCGGGTGCTGATGCAGGTCGGCCTGTTTCTGAATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTggg                                                                                                                              <  1:1115998/150‑1 (MQ=255)
                        gTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGt                                                                                                                            >  2:1175826/1‑150 (MQ=255)
                         ttCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTc                                                                                                                           >  1:2052552/1‑150 (MQ=255)
                            cTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTcgcg                                                                                                                        >  1:1571413/1‑150 (MQ=255)
                                tCGTCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTAAACTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGgcgc                                                                                                                    >  2:1521281/1‑150 (MQ=255)
                                tCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGgcgc                                                                                                                    >  1:1577893/1‑150 (MQ=255)
                                    cTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTg                                                                                                                <  1:104947/150‑1 (MQ=255)
                                     tATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTgg                                                                                                               <  1:1175826/150‑1 (MQ=255)
                                         tggtgCCGCGTTTCAGCCAGGGGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCACGTCGGCCTGTTTCTGCATGCCCAGCACGCGCCACTGGCGCTGGGCACCGTCGGTGGGGGCGCGGCGCTTTGgctgc                                                                                                           >  1:783256/1‑150 (MQ=255)
                                               cgcgTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGcc                                                                                                     <  1:1907046/150‑1 (MQ=255)
                                               cgcgTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGcc                                                                                                     >  2:1847267/1‑150 (MQ=255)
                                                gcgTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCg                                                                                                    >  1:98505/1‑150 (MQ=255)
                                                      cAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCAccc                                                                                              <  2:58116/150‑1 (MQ=255)
                                                        gCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCcgc                                                                                            >  1:1589781/1‑150 (MQ=255)
                                                         ccAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCcgcg                                                                                           >  2:1524429/1‑150 (MQ=255)
                                                         ccAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCcgcg                                                                                           >  2:2170482/1‑150 (MQ=255)
                                                               gtTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACg                                                                                     <  1:920235/150‑1 (MQ=255)
                                                                ttCGAAGGCATGGGCACCGAGCTACCCTGGATGTCACGGAGGCTGATGCAGGTCGGCCTGTTTCTGCCGGCCGAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCAGACACCCGCGGGTACgg                                                                                    >  1:1355552/1‑150 (MQ=255)
                                                                  cGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTa                                                                                  <  1:1874931/150‑1 (MQ=255)
                                                                   gAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTAc                                                                                 >  2:898/1‑150 (MQ=255)
                                                                   gAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGAACCGGCACCCACGTGTGCGGTAc                                                                                 >  2:2274233/1‑150 (MQ=255)
                                                                    aaGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACt                                                                                <  1:193408/150‑1 (MQ=255)
                                                                      ggCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTgg                                                                              <  1:460290/150‑1 (MQ=255)
                                                                       gCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTggg                                                                             <  1:1471407/150‑1 (MQ=255)
                                                                          tGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCgt                                                                                          >  1:1500409/1‑134 (MQ=255)
                                                                          tGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGcct                                                                          <  2:139434/150‑1 (MQ=255)
                                                                           gggCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGcctc                                                                         >  2:1354644/1‑150 (MQ=255)
                                                                           gggCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGGCGGCCTGTGTCGGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGcctc                                                                         >  2:1942312/1‑150 (MQ=255)
                                                                              cACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGTCCAGCAAGCGCCACTGGCGCTgg                                                                                                                                            <  2:1184856/80‑1 (MQ=255)
                                                                              cACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTgg                                                                                                                                            >  1:1184856/1‑80 (MQ=255)
                                                                              cACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTg                                                                      <  1:1748425/150‑1 (MQ=255)
                                                                               cccGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGc                                                                     <  2:1954333/149‑1 (MQ=255)
                                                                                  gAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTTGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAg                                                                  >  2:2395441/1‑150 (MQ=255)
                                                                                   aGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGc                                                                 >  1:1385447/1‑150 (MQ=255)
                                                                                   aGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGc                                                                 <  1:2170482/150‑1 (MQ=255)
                                                                                   aGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGGCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGc                                                                 <  2:268474/150‑1 (MQ=255)
                                                                                    gCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCt                                                                <  2:2251033/150‑1 (MQ=255)
                                                                                      tACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTgc                                                              >  1:1868401/1‑150 (MQ=255)
                                                                                      tACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTgc                                                              >  1:1867516/1‑150 (MQ=255)
                                                                                         ccTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGc                                                                            >  2:2403377/1‑133 (MQ=255)
                                                                                           tGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGAGTGTACGGTACTGGGCCTCCTGCCAGCTGCGgcgc                                                         >  2:1690311/1‑150 (MQ=255)
                                                                                               tGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGc                                                     >  2:2173408/1‑150 (MQ=255)
                                                                                               tGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCCGCGGCGCCTGc                                                     <  1:576383/150‑1 (MQ=255)
                                                                                                gTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGcc                                                    <  1:2173408/150‑1 (MQ=255)
                                                                                                   aCGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGc                                                 >  1:2104677/1‑150 (MQ=255)
                                                                                                       gggCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGcctcct                                                                       <  2:1094864/122‑1 (MQ=255)
                                                                                                       gTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGcctcct                                                                       >  1:1094864/1‑124 (MQ=255)
                                                                                                       gTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTg                                             >  1:727540/1‑150 (MQ=255)
                                                                                                        tGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTgg                                                                              >  1:673914/1‑116 (MQ=255)
                                                                                                        tGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTgg                                                                              <  2:673914/116‑1 (MQ=255)
                                                                                                        tGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGc                                            <  1:2340224/150‑1 (MQ=255)
                                                                                                          cTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCAc                                          <  2:1355552/150‑1 (MQ=255)
                                                                                                           tGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGgc                                                                                                              <  1:610276/81‑1 (MQ=255)
                                                                                                           tGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGgc                                                                                                              >  2:610276/1‑81 (MQ=255)
                                                                                                             aTGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTa                                                                                       <  1:1890026/102‑1 (MQ=255)
                                                                                                               gCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAgcg                                     >  2:771647/1‑150 (MQ=255)
                                                                                                                 aGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGGCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGGGCGTGCCCGCACTGCACCAGCGcc                                   >  2:1208636/1‑150 (MQ=255)
                                                                                                                  ggTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCt                                  >  1:1555608/1‑150 (MQ=255)
                                                                                                                  ggTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCt                                  >  2:1250958/1‑150 (MQ=255)
                                                                                                                      ggCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTgg                                                                              <  2:1440270/102‑1 (MQ=255)
                                                                                                                      ggCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTgg                                                                              >  1:1440270/1‑102 (MQ=255)
                                                                                                                        ccTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTgatgat                            >  2:2181892/1‑150 (MQ=255)
                                                                                                                          tGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCAGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGt                          >  2:1900710/1‑150 (MQ=255)
                                                                                                                            tttCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTAc                        <  1:2430119/150‑1 (MQ=255)
                                                                                                                            tttCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTAc                        <  1:905224/150‑1 (MQ=255)
                                                                                                                                  cATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTg                  <  2:1225383/150‑1 (MQ=255)
                                                                                                                                  cATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTg                  <  1:771647/150‑1 (MQ=255)
                                                                                                                                   gtGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTgg                 >  1:905199/2‑150 (MQ=255)
                                                                                                                                   gtGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTgg                 >  2:1540614/2‑150 (MQ=255)
                                                                                                                                   aTGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTgg                 <  2:1385447/150‑1 (MQ=255)
                                                                                                                                    tGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGTACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGgc                <  2:1577893/150‑1 (MQ=255)
                                                                                                                                      cccAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGgcgc                                                         <  1:1715807/107‑1 (MQ=255)
                                                                                                                                      cccAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGgcgc                                                         >  2:1715807/1‑107 (MQ=255)
                                                                                                                                       ccAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGcgcg             >  1:1061649/1‑150 (MQ=255)
                                                                                                                                        cAGCAAGCGCCACTGGCGCTGGGCACCATCGGTGGGGTCGCGGCGCTTTGGCTTCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGt            >  1:545280/1‑150 (MQ=255)
                                                                                                                                          gCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGttt          <  2:812525/150‑1 (MQ=255)
                                                                                                                                           ccaGCGCCAGTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCa                                                                   >  1:1626544/3‑92 (MQ=255)
                                                                                                                                           ccaGCGCCAGTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCa                                                                   <  2:1626544/90‑1 (MQ=255)
                                                                                                                                           ccaGCGCCAGTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTc         >  2:155959/3‑150 (MQ=255)
                                                                                                                                           ccaGCGCCAGTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTc         >  1:2086395/3‑150 (MQ=255)
                                                                                                                                            aaGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTCt        <  2:2243391/150‑1 (MQ=255)
                                                                                                                                              gcgcCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTCTAt      <  2:1571413/150‑1 (MQ=255)
                                                                                                                                              gcgcCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTCTAt      <  1:1847267/150‑1 (MQ=255)
                                                                                                                                                gcCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCGTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTCTAtcg    >  1:1413330/1‑150 (MQ=255)
                                                                                                                                                  cACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGGGTGTACGGTACGGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTCTAtcgct  <  1:1942312/150‑1 (MQ=255)
                                                                                                                                                  cACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTCTAtcgct  >  2:1053706/1‑150 (MQ=255)
                                                                                                                                                  cACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTCTAtcgct  >  2:886222/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTGGTGGGAGGCGGTGTGGTGCTGTTCCTGCTCGGCTATGTGGTGCCGCGTTTCAGCCAGGTGTTCGAAGGCATGGGCACCGAGCTACCCTGGCTGTCACGGGTGCTGATGCAGGTCGGCCTGTTTCTGCATGCCCAGCAAGCGCCACTGGCGCTGGGCACCGTCGGTGGGGTCGCGGCGCTTTGGCTGCTGCGCCGGCACCCGCGTGTACGGTACTGGGCCTCCTGCCAGCTGCGGCGCCTGCCCGCACTGCACCAGCGCCTGATGATGTACGAGCTGGCGCGTTTCTATCGCT  >  NC_002947/3876083‑3876378

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