Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,121,458 C→A 7.0% R298R (CGC→CGA PP_0982 → YjgP/YjgQ family permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,121,4580CA 7.0% 345.4 / 6.7 107R298R (CGC→CGAPP_0982YjgP/YjgQ family permease
Reads supporting (aligned to +/- strand):  ref base C (48/51);  new base A (4/4);  total (52/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCCTGCTGATCG  >  NC_002947/1121311‑1121606
                                                                                                                                                   |                                                                                                                                                    
gACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTTGACCCGCGACa                                                                                                                                                    <  2:229810/150‑1 (MQ=255)
   cGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCGAACTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAggg                                                                                                                                                 >  1:1776628/1‑150 (MQ=255)
   cGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCGGTCGCGGGTCAACCCGCCCCAggg                                                                                                                                                 >  1:1005622/1‑150 (MQ=255)
    gAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGc                                                                                                                                                >  1:2297954/1‑150 (MQ=255)
     aGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGcc                                                                                                                                               <  1:1029236/150‑1 (MQ=255)
      gcgcGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCg                                                                                                                                              <  2:865691/150‑1 (MQ=255)
         cgAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCtt                                                                                                                                           >  2:101011/1‑150 (MQ=255)
          gAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGCGCCGCTTc                                                                                                                                          <  2:2144009/150‑1 (MQ=255)
            aGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCt                                                                                                                                        <  2:1443842/150‑1 (MQ=255)
               gcTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCaa                                                                                                                                     >  1:1745880/1‑150 (MQ=255)
               gcTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCCGATCTTGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCaa                                                                                                                                     <  1:1779730/150‑1 (MQ=255)
                cTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAg                                                                                                                                    <  1:1267187/150‑1 (MQ=255)
                 ttCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgc                                                                                                                                   >  1:1465961/1‑150 (MQ=255)
                 ttCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgc                                                                                                                                   >  1:205190/1‑150 (MQ=255)
                 ttCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGACAgcg                                                                                                                                                 <  2:1242685/136‑3 (MQ=255)
                    ccACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGTGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgctgc                                                                                                                                <  1:1009235/150‑1 (MQ=255)
                     cACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgctgct                                                                                                                               >  2:1616669/1‑150 (MQ=255)
                        cTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGcc                                                                                                                            <  2:1484446/150‑1 (MQ=255)
                        cTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGcc                                                                                                                            >  1:1484446/1‑150 (MQ=255)
                         tCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCg                                                                                                                           >  2:2179295/1‑150 (MQ=255)
                          ccGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCgg                                                                                                                          <  2:948206/150‑1 (MQ=255)
                             aaCTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCaa                                                                                                                       <  2:806783/150‑1 (MQ=255)
                                tGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAAttc                                                                                                                    <  1:2115327/150‑1 (MQ=255)
                                    cGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATtcttct                                                                                                                <  2:1872408/150‑1 (MQ=255)
                                    cGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATtcttct                                                                                                                <  1:2205334/150‑1 (MQ=255)
                                        aaGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTAGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTAAAGCTGCTGCCGGCAATTCTTCTGTAc                                                                                                            >  1:2231822/1‑150 (MQ=255)
                                                ccTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTa                                                                                                    >  1:1047499/1‑150 (MQ=255)
                                                  tGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACGTGGCCTAcc                                                                                                  >  1:198180/1‑150 (MQ=255)
                                                  tGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTAcc                                                                                                  <  1:851703/150‑1 (MQ=255)
                                                     gTGAGCGCCCCGAGCTGCAAGGGCGGTTGTCGCTGCCGTTCCTGGTGTTCGTCGCCCCCCTGCTGGCGGTGCCGCTGTCGCGGGTCACCCCGCGCCAGGGCCCCTCCCCAACGCTCCAGCCGGAACTCCTTCTCTACTTGGCATACCTGa                                                                                               <  2:1786011/150‑1 (MQ=255)
                                                      tGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGAc                                                                                              <  2:1047499/150‑1 (MQ=255)
                                                         gcgcgcCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAAt                                                                                           <  1:386934/150‑1 (MQ=255)
                                                            cgcCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACGTGGCCTACCTGACAATGCt                                                                                        <  2:198180/150‑1 (MQ=255)
                                                             gcCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTg                                                                                       <  1:618544/150‑1 (MQ=255)
                                                                   cTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCgg                                                                                                                          <  2:281299/109‑1 (MQ=255)
                                                                   cTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTcc                                                                                 <  1:2179295/150‑1 (MQ=255)
                                                                   cTGCAAGGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTCATGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTcc                                                                                 <  1:1218428/150‑1 (MQ=255)
                                                                      cAATGGCGGTTGTCGCTGCCGATCCTGNTGTTCGTCGTCANCCTGCTNGCGGTGCCGCTGTAAGGGATCAAAACAGCACACGTCTGAACTCCAGTGCTGCTTCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTa                                                                              >  1:1694606/1‑150 (MQ=255)
                                                                      cAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGc                                                                                                       >  1:46788/1‑125 (MQ=255)
                                                                      cAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGc                                                                                                       <  2:46788/125‑1 (MQ=255)
                                                                          ggCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTAcgcg                                                                          <  1:2291056/150‑1 (MQ=255)
                                                                                tGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATtcttct                                                                                                                >  1:1319975/1‑106 (MQ=255)
                                                                                tGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATtcttct                                                                                                                <  2:1319975/106‑1 (MQ=255)
                                                                                  tCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGATAAAACGGATTTCCGTACGCGGCGCCCTc                                                                  >  1:1219/1‑150 (MQ=255)
                                                                                  tCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTc                                                                  >  1:1123155/1‑150 (MQ=255)
                                                                                    gCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCt                                                                                        <  1:416327/126‑1 (MQ=255)
                                                                                    gCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCga                                                                <  2:205190/150‑1 (MQ=255)
                                                                                    gCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCga                                                                >  1:796572/1‑150 (MQ=255)
                                                                                    gCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCACAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCga                                                                <  2:1917792/150‑1 (MQ=255)
                                                                                       gCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGaa                                                             <  2:1745880/150‑1 (MQ=255)
                                                                                        ccGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAg                                                            <  1:648245/150‑1 (MQ=255)
                                                                                         cGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCACCCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTCTTCATCGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAgg                                                           <  2:87019/150‑1 (MQ=255)
                                                                                            tCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAAttctt                                                                                                                  <  1:2093319/92‑1 (MQ=255)
                                                                                                ggTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGc                                                                                                       <  1:457110/99‑1 (MQ=255)
                                                                                                ggTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGc                                                                                                       >  2:457110/1‑99 (MQ=255)
                                                                                                ggTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACt                                                    >  2:327164/1‑150 (MQ=255)
                                                                                                    ttCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCCAGGACACTTTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCg                                                >  1:145272/1‑150 (MQ=255)
                                                                                                    ttCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAAACCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCg                                                >  1:516856/1‑150 (MQ=255)
                                                                                                     tcgtcgTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGa                                               <  2:160073/150‑1 (MQ=255)
                                                                                                       gtcgtcACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATc                                             >  1:106721/1‑150 (MQ=255)
                                                                                                        tcgtcACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCg                                            >  2:1578039/1‑150 (MQ=255)
                                                                                                          gtcACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGcc                                          >  2:642413/1‑150 (MQ=255)
                                                                                                           tcACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGccc                                         >  2:1789383/1‑150 (MQ=255)
                                                                                                               ccTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTAc                                                                                                   >  2:393955/1‑88 (MQ=255)
                                                                                                               ccTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTAc                                                                                                   <  1:393955/88‑1 (MQ=255)
                                                                                                               ccTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTgtg                                     >  1:2146716/1‑148 (MQ=255)
                                                                                                               ccTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTggg                                     >  1:250579/1‑150 (MQ=255)
                                                                                                                ctgctgGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAGACTGCCGATCGCCCTGGGc                                    >  2:1533787/1‑150 (MQ=255)
                                                                                                                ctgctgGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGc                                    >  1:1922478/1‑150 (MQ=255)
                                                                                                                    tgGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCATGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCAtgt                                >  2:795211/1‑150 (MQ=255)
                                                                                                                    tgGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCAtgt                                <  2:2146716/150‑1 (MQ=255)
                                                                                                                      gCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATgtgg                              <  1:101011/150‑1 (MQ=255)
                                                                                                                       cGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATgtggt                             <  2:2231822/150‑1 (MQ=255)
                                                                                                                         gTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATgtggttg                           <  1:1616669/150‑3 (MQ=255)
                                                                                                                         gTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGtggtgg                           <  1:642413/150‑1 (MQ=255)
                                                                                                                            ccGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCATGAGCATGTGGTCGGta                        <  2:1788469/150‑2 (MQ=255)
                                                                                                                              cctggcgcgGGTTGACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAgg                                                           <  1:75410/108‑1 (MQ=255)
                                                                                                                              ccTGGCGCGGGTTGACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTcc                                                                                 >  1:1405143/2‑91 (MQ=255)
                                                                                                                              ccTGGCGCGGGTTGACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAgg                                                           >  2:75410/2‑113 (MQ=255)
                                                                                                                               ctggcgcgGGTTGACCCGCGACAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCaaa                                                      <  1:2056926/113‑1 (MQ=255)
                                                                                                                               cTGGCGCGGGTTGACCCGCGACAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCaaa                                                      >  2:2056552/1‑117 (MQ=255)
                                                                                                                               cTGGCGCGGGTTGACCCGCGACAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCaaa                                                      >  2:2056926/1‑117 (MQ=255)
                                                                                                                                tGTCGCGGGTCAACCCGCGCCATGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACg                    <  1:795211/150‑1 (MQ=255)
                                                                                                                                tGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTg                                       <  1:1968809/131‑1 (MQ=255)
                                                                                                                                tGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTg                                       >  2:1968809/1‑131 (MQ=255)
                                                                                                                                  tCGCGGGTCAACCCGCTCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGc                  <  1:1789383/150‑1 (MQ=255)
                                                                                                                                  tCGCGGGTCAACCCGCGCCAGGGAGGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCTACGGc                  <  2:1123155/150‑1 (MQ=255)
                                                                                                                                    gcgGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCt                >  2:2202090/1‑150 (MQ=255)
                                                                                                                                        gTCAACCCGCGCCAGGGCCGCTTCCTCACGCTGCTGCCGGCAATTCTTCTGTACCCCAGCGACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTc            <  1:1194415/150‑1 (MQ=255)
                                                                                                                                        gTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCGGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTc            >  1:1891473/1‑150 (MQ=255)
                                                                                                                                            aCCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCctgc        <  1:983186/150‑1 (MQ=255)
                                                                                                                                             cccGCGCCCAGGACGCCTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCctgct       >  2:1264614/1‑150 (MQ=255)
                                                                                                                                             cccGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTAc                                                                                                   <  1:1536759/58‑1 (MQ=255)
                                                                                                                                             cccGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTAc                                                                                                   >  2:1536759/1‑58 (MQ=255)
                                                                                                                                             cccGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCctgct       >  2:54640/1‑150 (MQ=255)
                                                                                                                                             cccGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCctgct       >  1:1542246/1‑150 (MQ=255)
                                                                                                                                              ccGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCctgctg      >  2:346453/1‑150 (MQ=255)
                                                                                                                                              acgcaccAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGcc                 <  1:1475952/134‑1 (MQ=255)
                                                                                                                                              acgcaCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGcc                 >  2:1475952/2‑139 (MQ=255)
                                                                                                                                                  gcCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCCTGCTGATCg  >  2:1210959/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTGTCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCCTGCTGATCG  >  NC_002947/1121311‑1121606

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