Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,773,965 A→T 6.9% T209S (ACA→TCA)  PP_5743 → TonB‑dependent receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,773,9650AT 6.9% 317.8 / 11.4 102T209S (ACA→TCA) PP_5743TonB‑dependent receptor protein
Reads supporting (aligned to +/- strand):  ref base A (49/46);  new base T (3/4);  total (52/50)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 104 total aligned reads displayed.

GACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCTTTTCTAC  >  NC_002947/3773820‑3774114
                                                                                                                                                 |                                                                                                                                                     
gACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATg                                                                                                                                                   <  2:964419/150‑1 (MQ=255)
 aCGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATgg                                                                                                                                                  <  2:1416693/150‑1 (MQ=255)
    cccAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCg                                                                                                                                               >  1:1028345/1‑150 (MQ=255)
     ccAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATgg                                                                                                                                                  >  2:2372086/1‑146 (MQ=255)
     ccAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATgg                                                                                                                                                  <  1:2372086/146‑1 (MQ=255)
     ccAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGCGCCCATGTCATGGGcac                                                                                                                                              <  2:2305306/150‑3 (MQ=255)
       aaGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGTCATGGGcacgt                                                                                                                                            <  1:2087301/150‑5 (MQ=255)
       aaGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACAt                                                                                                                                                    >  1:1089916/1‑142 (MQ=255)
       aaGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACAt                                                                                                                                                    <  2:1089916/142‑1 (MQ=255)
       aaGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCt                                                                                                                                            <  1:443671/150‑1 (MQ=255)
             gCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATAc                                                                                                                                      <  1:2225987/150‑1 (MQ=255)
             gCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATAc                                                                                                                                      <  1:1868086/150‑1 (MQ=255)
                    aGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGCGCCCATGTCATGGGcacgt                                                                                                                                            >  2:839844/1‑133 (MQ=255)
                    aGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGCGCCCATGTCATGGGcacgt                                                                                                                                            <  1:839844/137‑5 (MQ=255)
                      ggCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGt                                                                                                                             >  2:1783126/1‑150 (MQ=255)
                          aCAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCt                                                                                                                         >  2:583679/1‑150 (MQ=255)
                            aaCCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTgg                                                                                                                       >  2:859332/1‑150 (MQ=255)
                             accaccATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgt                                                                                                                      <  2:2051405/150‑1 (MQ=255)
                               caccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATgg                                                                                                                                                  >  1:569399/1‑120 (MQ=255)
                               caccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATgg                                                                                                                                                  <  2:569399/120‑1 (MQ=255)
                               caccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgtgt                                                                                                                    >  1:116555/1‑150 (MQ=255)
                               caccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgtgt                                                                                                                    <  2:307482/150‑1 (MQ=255)
                               caccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTATCTGgtgt                                                                                                                    >  2:1310722/1‑150 (MQ=255)
                                accaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCt                                                                                                                         <  2:1202063/144‑1 (MQ=255)
                                accaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCt                                                                                                                         >  1:1202063/1‑144 (MQ=255)
                                 ccaTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTaa                                                                                                                  <  2:2123666/150‑1 (MQ=255)
                                  caTGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAAc                                                                                                                 >  2:1521148/1‑150 (MQ=255)
                                     ggggAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGt                                                                                                              >  1:410912/1‑150 (MQ=255)
                                            tGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGc                                                                                                       >  1:1313322/1‑150 (MQ=255)
                                              cAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAg                                                                                                     >  2:1864998/1‑150 (MQ=255)
                                                aCTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTg                                                                                                                        >  1:2425726/1‑129 (MQ=255)
                                                aCTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTg                                                                                                                        <  2:2425726/129‑1 (MQ=255)
                                                  tCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTg                                                                                                 >  1:1699006/1‑150 (MQ=255)
                                                  tCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTg                                                                                                 >  1:2142226/1‑150 (MQ=255)
                                                     ccGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCAAGGTTTCTGGTGTAACGGTCTGTCGCAGCCTgcgc                                                                                              <  2:2404697/150‑1 (MQ=255)
                                                         gggcAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCgcggcg                                                                                          <  1:2240730/149‑1 (MQ=255)
                                                           gcAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTa                                                                                        >  2:1208948/1‑150 (MQ=255)
                                                             aGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAAt                                                                                      <  1:2306940/150‑1 (MQ=255)
                                                                ttCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACCACCACACGTCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGc                                                                                   <  2:71923/150‑1 (MQ=255)
                                                                    gCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTAAAGAATCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAg                                                                               >  1:1636372/1‑150 (MQ=255)
                                                                        ggCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCgcgg                                                                                            >  1:648354/1‑133 (MQ=255)
                                                                        ggCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCgcgg                                                                                            <  2:648354/133‑1 (MQ=255)
                                                                        ggCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGt                                                                           <  1:859332/150‑1 (MQ=255)
                                                                          ccACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTgg                                                                                                                       >  1:1141811/1‑104 (MQ=255)
                                                                          ccACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTgg                                                                                                                       <  2:1141811/104‑1 (MQ=255)
                                                                           cACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgt                                                                                                                      >  1:827764/1‑104 (MQ=255)
                                                                           cACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGgt                                                                                                                      <  2:827764/104‑1 (MQ=255)
                                                                            atcGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAg                                                                       >  2:1073249/3‑150 (MQ=255)
                                                                            aCCGAACAAGCCAGCCTGATTGGCCGCGAATGCGATTTCGACTACCACACCCCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAg                                                                       <  1:815867/150‑1 (MQ=255)
                                                                            aCCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAg                                                                       >  2:2288859/1‑150 (MQ=255)
                                                                               gAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGttt                                                                                                                           <  1:1591727/95‑1 (MQ=255)
                                                                               gAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGttt                                                                                                                           >  2:1591727/1‑95 (MQ=255)
                                                                                    agccagccTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGa                                                               <  1:1853957/150‑1 (MQ=255)
                                                                                     gccagccTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGAc                                                              >  1:1503097/1‑150 (MQ=255)
                                                                                      ccagccTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACt                                                             >  1:741339/1‑150 (MQ=255)
                                                                                         gccTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTAc                                                          <  1:491765/150‑1 (MQ=255)
                                                                                          ccTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTaa                                                                                       <  2:1146681/120‑1 (MQ=255)
                                                                                          ccTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTaa                                                                                       >  1:1146681/1‑120 (MQ=255)
                                                                                             gATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAAc                                                                                                                 >  2:2150185/1‑91 (MQ=255)
                                                                                                 ggCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCCCGCCATCTGGTGGAACGGGCGGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTcgc                                                  >  1:2121510/1‑150 (MQ=255)
                                                                                                  gCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTcgcg                                                 <  1:844946/150‑1 (MQ=255)
                                                                                                   ccGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCgg                                                <  2:116555/150‑1 (MQ=255)
                                                                                                   ccGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCgg                                                >  1:561633/1‑150 (MQ=255)
                                                                                                    cgcgAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCggg                                               >  1:2213363/1‑150 (MQ=255)
                                                                                                    cgcgAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTCGTGTAACGGTCTGTCCGTTCCTGTGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCggg                                               <  1:888080/150‑1 (MQ=255)
                                                                                                      cgAGGGCGATGTCGACGACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGGAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCa                                             <  2:428295/150‑1 (MQ=255)
                                                                                                      cgAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCa                                             <  1:1377817/150‑1 (MQ=255)
                                                                                                        aaGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGAt                                                                     <  1:751262/124‑1 (MQ=255)
                                                                                                        aaGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGAt                                                                     >  2:751262/1‑124 (MQ=255)
                                                                                                          ggCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTc                                         <  2:805461/150‑1 (MQ=255)
                                                                                                           gCGATTTCGACTACCACACACCCAGGTACCTCCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTcc                                        >  1:62750/1‑150 (MQ=255)
                                                                                                            cGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCa                                       <  2:1699006/150‑1 (MQ=255)
                                                                                                             gATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGc                                                                                                    <  2:915124/88‑1 (MQ=255)
                                                                                                              aTTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTTTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGc                                     >  2:325874/1‑150 (MQ=255)
                                                                                                               tttCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGcc                                    >  2:753981/1‑149 (MQ=255)
                                                                                                                 tCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAAc                                  <  1:1521148/150‑1 (MQ=255)
                                                                                                                 tCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAAc                                  <  2:741339/150‑1 (MQ=255)
                                                                                                                  cGACTACCACACACCCAGGCACGTGCCCATGACATGTGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCTGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACt                                 <  1:110711/150‑1 (MQ=255)
                                                                                                                   gACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTg                                >  2:1860420/1‑150 (MQ=255)
                                                                                                                   gACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTg                                >  2:1815080/1‑150 (MQ=255)
                                                                                                                      tACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTtg                             >  1:2418738/1‑150 (MQ=255)
                                                                                                                       aCCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTtgc                            <  2:2213363/150‑1 (MQ=255)
                                                                                                                        ccACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTtgct                           <  2:561633/150‑1 (MQ=255)
                                                                                                                         cacacaCCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTtgctg                          >  1:1803783/1‑150 (MQ=255)
                                                                                                                            acacCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGGGTAACGGGCTGTCGCAGACTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCAc                       >  1:484806/1‑150 (MQ=255)
                                                                                                                              acCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACtt                     <  2:1908032/150‑1 (MQ=255)
                                                                                                                               cccAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGATTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTc                    <  2:1313322/150‑1 (MQ=255)
                                                                                                                                ccAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCg                   >  2:2184960/1‑150 (MQ=255)
                                                                                                                                ccAGGCACGTGCCCATCACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCg                   >  2:147432/1‑150 (MQ=255)
                                                                                                                                 cAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTAGCAGCCTGCCAAATTCAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGc                  >  1:1484291/1‑150 (MQ=255)
                                                                                                                                 cAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGGGTAACGGTCTGTCGCAGCCGTCGCGGCGTCATGGGTTAGGAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGc                  >  1:767395/1‑150 (MQ=255)
                                                                                                                                      agacgCCCATGTCATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCg             >  2:137669/5‑150 (MQ=255)
                                                                                                                                        acgCCCATGTCATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGc           >  1:1837703/3‑150 (MQ=255)
                                                                                                                                           cccATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCtt        <  2:2121510/150‑1 (MQ=255)
                                                                                                                                            ccATGTCATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCttt       <  2:1837703/150‑1 (MQ=255)
                                                                                                                                            ccATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCttt       >  1:1148070/1‑150 (MQ=255)
                                                                                                                                             cATGACATGGGCGTCTGCATACTGCGACGGTTTGTGGTGTAAGGGGGTGTTGCAGACTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCtttt      >  1:2263381/1‑150 (MQ=255)
                                                                                                                                             cATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTCCACTTCGCTCTCGGCCtttt      <  1:1208948/150‑1 (MQ=255)
                                                                                                                                                 aCATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATTGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCTTTTCTAc  >  2:1399797/1‑150 (MQ=255)
                                                                                                                                                 aCATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTCCCGATCGTGGACTTACCCGGTAGCGTGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCTTTTCTAc  <  1:1864998/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GACGCCCAAGGTTGCCGAGCAGGGCGACAACCACCATGGGGAGCTGCAACTCACCGGCGGCAGCTTCCGCACGGCCACCGAACAAGCCAGCCTGATTGGCCGCGAAGGCGATTTCGACTACCACACACCCAGGCACGTGCCCATGACATGGGCGTCTGCATACTGCCACGGTTTCTGGTGTAACGGTCTGTCGCAGCCTGCGCGGCGTAATGGCTTAGCAGTGCAGATCGTGGACTTACCCGGTCGCGGGCAGGTCCAGCAACTGTTGCTGCACTTCGCTCTCGGCCTTTTCTAC  >  NC_002947/3773820‑3774114

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: