Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,633,969 C→T 6.9% intergenic (+21/‑143) lepB → / → rnc signal peptidase I/ribonuclease 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,633,9690CT 6.9% 213.1 / 3.8 73intergenic (+21/‑143)lepB/rncsignal peptidase I/ribonuclease 3
Reads supporting (aligned to +/- strand):  ref base C (31/36);  new base T (2/3);  total (33/40)
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

CCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGT  >  NC_002947/1633821‑1634117
                                                                                                                                                    |                                                                                                                                                     
ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGTATGAGCTGGGCAGAGCCTAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                      <  1:1590892/150‑1 (MQ=255)
ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                      >  2:1767057/1‑150 (MQ=255)
ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                      >  1:401597/1‑150 (MQ=255)
ccAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCcagca                                                                                                                                                      >  1:1366916/1‑150 (MQ=255)
 cAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCACGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTCTCcagcag                                                                                                                                                     >  1:1599495/1‑150 (MQ=255)
  aaGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGa                                                                                                                                                    <  2:1798264/150‑1 (MQ=255)
  aaGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGa                                                                                                                                                    >  1:1042424/1‑150 (MQ=255)
  aaGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGa                                                                                                                                                    >  2:102741/1‑150 (MQ=255)
    ggAACTGCACGGCATGGTTCCGGACCTGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACa                                                                                                                                                  <  1:819424/150‑1 (MQ=255)
    ggAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACa                                                                                                                                                  <  2:1366916/150‑1 (MQ=255)
         tGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGcc                                                                                                                                             >  1:1904548/1‑150 (MQ=255)
         cacccgatcgttgttcaacgcgttgttcaaatcaacatcctgaaagactgtttggctgtggtccgccgattgatgtgtggcgaatctcccgatccctgagaacacagcctatgtcagaaatgcatTCGGCGCTGTCCAGCAGACAGCGcc                                                                                                                                             <  2:105262‑M1/25‑1 (MQ=255)
                gaTGGTGCCGGACCGGAACATTGACGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCa                                                                                                                                      <  2:159502/149‑1 (MQ=255)
                cATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCACCAGGTGACCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCa                                                                                                                                      <  2:553375/150‑1 (MQ=255)
                 aTGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCAt                                                                                                                                     <  1:1473946/150‑1 (MQ=255)
                   ggTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGCCAGGCCCGAATGCAttt                                                                                                                                   <  2:1050510/150‑1 (MQ=255)
                    gTTCCGGTCCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGcc                                                                                                                                             <  2:1081518/139‑1 (MQ=255)
                    gTTCCGGTCCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGcc                                                                                                                                             >  1:1081518/1‑139 (MQ=255)
                       ccggaccggaACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGa                                                                                                                               <  2:1042424/150‑1 (MQ=255)
                                 aCATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtgt                                                                                                                     <  2:401597/150‑1 (MQ=255)
                                  cATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCtgtg                                                                                                                    <  2:1904548/150‑1 (MQ=255)
                                        cGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCa                                                                                                              <  2:96130/150‑1 (MQ=255)
                                         ggCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAg                                                                                                             >  2:598732/1‑150 (MQ=255)
                                               gCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCg                                                                                                       <  1:1798689/150‑1 (MQ=255)
                                                 cTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTAAGCCACTTGCCCAATCCACCGCGGGTTGGCCTGATCCATTTATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCggg                                                                                                     <  2:707297/150‑1 (MQ=255)
                                                         tgtgGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGGCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGc                                                                                             <  2:1682002/150‑1 (MQ=255)
                                                           tgGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTNCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCca                                                                                          <  1:1767057/150‑1 (MQ=255)
                                                             gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAgcgc                                                                                                                                              <  2:1224237/97‑1 (MQ=255)
                                                             gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAgcgc                                                                                                                                              >  1:1224237/1‑97 (MQ=255)
                                                             gATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCcaca                                                                                        >  1:1386092/1‑150 (MQ=255)
                                                                   tggtccgccgattgatgtgtggcgaatctcccgatccctgagaacacagcctatgtcagaaatgcatTCGGCGCTGTCcagc                                                                                                                                                       >  1:230124‑M1/68‑82 (MQ=255)
                                                                   tggtccgccgattgatgtgtggcgaatctcccgatccctgagaacacagcctatgtcagaaatgcatTCGGCGCTGTCcagc                                                                                                                                                       <  2:230124‑M1/15‑1 (MQ=255)
                                                                         agagCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATcggcg                                                                            <  1:1898189/150‑1 (MQ=255)
                                                                             cccAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTAGCGGGGATCGGGAGATTCGCCACACATCAATCGGCGGTcc                                                                        >  1:1373603/1‑150 (MQ=255)
                                                                              ccAAACTCAGCCACTTGCCCAACCTGTCTCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACca                                                                       >  1:346214/1‑150 (MQ=255)
                                                                               cAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGc                                                                                                                                       <  2:1608109/86‑1 (MQ=255)
                                                                               cAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGc                                                                                                                                       >  1:1608109/1‑86 (MQ=255)
                                                                               cAAACTCAGCCACTTCGACAACCTGTCGCTGGATGGACATATTAAATAATACCCCTCGGCGATATCCAGAATACAGCGTAAAATGCATTTTTGACATAAGCTGTGTGCTCAGGGATCGGGAGATTCGCCACATAAAAATCGGCGGACcac                                                                      <  2:121748/150‑1 (MQ=255)
                                                                                aaaCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACcaca                                                                     >  2:280922/1‑150 (MQ=255)
                                                                                  aCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATCCCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGc                                                                   <  2:297993/150‑1 (MQ=255)
                                                                                  aCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAg                                                                                                             <  1:1489264/109‑1 (MQ=255)
                                                                                  aCTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAg                                                                                                             >  2:1489264/1‑109 (MQ=255)
                                                                                    tCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACAt                                                                                     <  2:165031/130‑1 (MQ=255)
                                                                                    tCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACAt                                                                                     >  1:165031/1‑130 (MQ=255)
                                                                                        ccACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAg                                                                                                             >  1:1820644/1‑103 (MQ=255)
                                                                                        ccACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAg                                                                                                             <  2:1820644/103‑1 (MQ=255)
                                                                                        ccACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGc                                                                                             <  1:286905/119‑1 (MQ=255)
                                                                                        ccACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGc                                                                                             >  2:286905/1‑119 (MQ=255)
                                                                                        ccACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACa                                                             >  2:1621021/1‑150 (MQ=255)
                                                                                              gCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCa                                                                                                                                      <  1:1880408/72‑1 (MQ=255)
                                                                                              gCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCa                                                                                                                                      >  2:1880408/1‑72 (MQ=255)
                                                                                               cccAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTc                                                      <  1:102741/150‑1 (MQ=255)
                                                                                                       gTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGtt                                              <  1:598732/150‑1 (MQ=255)
                                                                                                        tCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCaaa                                                               >  1:1013815/1‑132 (MQ=255)
                                                                                                        tCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCaaa                                                               <  2:1013815/132‑1 (MQ=255)
                                                                                                        ctgagaacacagcctatgtcagaaatgcatTCGGCGCTGTCTGCTGGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTg                                             <  1:105262‑M1/120‑1 (MQ=255)
                                                                                                         tgagaacacagcctatgtcagaaatgcatTCGGCGCTGTCTGCTGGACAGCGCCGATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGAt                                           <  2:378181‑M1/121‑1 (MQ=255)
                                                                                                          gcgGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGAt                                           >  2:180019/1‑150 (MQ=255)
                                                                                                          gagaacacagcctatgtcagaaatgcatTCGGCGCTGTCTGCTGGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACa                                                             >  1:1906673‑M1/29‑132 (MQ=255)
                                                                                                          gagaacacagcctatgtcagaaatgcatTCGGCGCTGTCTGCTGGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACa                                                             <  2:1906673‑M1/104‑1 (MQ=255)
                                                                                                          gagaacacagcctatgtcagaaatgcatTCGGCGCTGTCTGCTGGACAGCGCCGATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGAtt                                          >  1:378181‑M1/29‑150 (MQ=255)
                                                                                                              gTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGa                                       >  2:1185240/1‑150 (MQ=255)
                                                                                                                  gCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCTGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaacaa                                   >  1:1666937/1‑150 (MQ=255)
                                                                                                                  gCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaacaa                                   >  1:1225204/1‑150 (MQ=255)
                                                                                                                  gCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaacaa                                   >  1:945574/1‑150 (MQ=255)
                                                                                                                  gCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGAACGGGTGCCCGCCTAGACATCAATCGGCTGACCACAGCCAAACAGTCTTTCAGGATGATGATTTGaacaa                                   >  1:1684035/1‑150 (MQ=255)
                                                                                                                   ccTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATcggc                                                                             >  1:713369/1‑107 (MQ=255)
                                                                                                                   ccTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATcggc                                                                             <  2:713369/107‑1 (MQ=255)
                                                                                                                   ccTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCa                                                     >  2:1730710/1‑131 (MQ=255)
                                                                                                                   ccTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCa                                                     <  1:1730710/131‑1 (MQ=255)
                                                                                                                         ccATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTg                            >  2:1642305/1‑150 (MQ=255)
                                                                                                                           aTTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCgg                                                                                                      <  2:1642437/75‑1 (MQ=255)
                                                                                                                           aTTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCgg                                                                                                      >  1:1642437/1‑75 (MQ=255)
                                                                                                                           aTTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGaa                          <  1:180019/150‑1 (MQ=255)
                                                                                                                                  cccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGAt                   <  1:280922/150‑1 (MQ=255)
                                                                                                                                   ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaacaac                                  >  1:1419086/1‑134 (MQ=255)
                                                                                                                                   ccATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGaacaac                                  <  2:1419086/134‑1 (MQ=255)
                                                                                                                                    cATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCg                 >  1:936553/1‑150 (MQ=255)
                                                                                                                                        GGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGAtt                                          >  1:1095268‑M1/1‑121 (MQ=255)
                                                                                                                                        GGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGAtt                                          <  2:1095268‑M1/121‑1 (MQ=255)
                                                                                                                                                   gcagACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGt  <  1:1621021/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCAAGGAACTGCACGGCATGGTTCCGGACCGGAACATTGTCGGCAAGGCCTTCGCGGTGTGGATGAGCTGGCCAGAGCCCAAACTCAGCCACTTGCCCAACCTGTCGCGGGTTGGCCTGATCCATTGATACCCATCGGCGCTGTCCAGCAGACAGCGCCGAATGCATTTCTGACATAGGCTGTGTTCTCAGGGATCGGGAGATTCGCCACACATCAATCGGCGGACCACAGCCAAACAGTCTTTCAGGATGTTGATTTGAACAACGCGTTGAACAACGATCGGGTGGCCGCATGAGT  >  NC_002947/1633821‑1634117

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: