Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,121,458 C→A 6.3% 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 6.3% 304.5 / 6.9 95R298R (CGC→CGAPP_0982YjgP/YjgQ family permease
Reads supporting (aligned to +/- strand):  ref base C (43/46);  new base A (3/3);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

GTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCCTGCTGATCG  >  NC_002947/1121309‑1121606
                                                                                                                                                       |                                                                                                                                                     
gTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCgcgc                                                                                                                                                       >  1:2379659/1‑150 (MQ=255)
    cctaGCTCGAAGCGCGTCCCGCCTCCGAACTGATCTGCAAGGATGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAgg                                                                                                                                                  <  1:1102448/147‑1 (MQ=255)
    ccGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAgg                                                                                                                                                  <  2:2292305/150‑1 (MQ=255)
         cgcgAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTTGACCCGCGACAgcggcac                                                                                                                                              <  1:726460/150‑7 (MQ=255)
            gAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCCTG‑‑GCGCGGGTTGACCCGCGACAgcggcacc                                                                                                                                             <  2:583080/147‑8 (MQ=255)
            gAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCCTG‑‑GCGCGGGTTGACCCGCGACAgcggcacc                                                                                                                                             >  1:583080/1‑140 (MQ=255)
             aaGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑GCGCGGGTTGACCCGCGACAgc                                                                                                                                                   >  1:1196468/1‑140 (MQ=255)
             aaGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑GCGCGGGTTGACCCGCGACAgc                                                                                                                                                   <  2:1196468/141‑2 (MQ=255)
              aGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGAGGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCCGCCCGCGCCCGGGCCGCTTCCt                                                                                                                                        >  2:316491/1‑150 (MQ=255)
                 gcTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCaa                                                                                                                                     <  2:621340/150‑1 (MQ=255)
                 gcTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCaa                                                                                                                                     <  2:1408544/150‑1 (MQ=255)
                  cTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAg                                                                                                                                    >  2:1446325/1‑150 (MQ=255)
                  cTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAg                                                                                                                                    <  2:57860/150‑1 (MQ=255)
                   ttCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgc                                                                                                                                   >  2:290039/1‑150 (MQ=255)
                      ccACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgctgc                                                                                                                                <  1:1316668/150‑1 (MQ=255)
                      ccACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgctgc                                                                                                                                >  1:157558/1‑150 (MQ=255)
                       cACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCACGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgctgct                                                                                                                               <  2:1304133/150‑1 (MQ=255)
                          cTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAggg                                                                                                                                                 >  1:1794843/1‑129 (MQ=255)
                          cTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAggg                                                                                                                                                 <  2:1794843/129‑1 (MQ=255)
                            ccGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCgg                                                                                                                          <  2:1347100/150‑1 (MQ=255)
                               aaCTGTTCGGCAAGGAGGCCCTGCGTAAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGAGCCTGGTGTTCGTCGTCACCCTGCTGGGGGTACTGCTG‑‑TCGCGCGTCACCCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCaa                                                                                                                       <  2:497126/150‑1 (MQ=255)
                                    aTCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAAttctt                                                                                                                  <  1:412453/150‑1 (MQ=255)
                                      cGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCGGCCCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATtcttct                                                                                                                >  1:22641/1‑150 (MQ=255)
                                      cGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATtcttct                                                                                                                <  2:1249246/150‑1 (MQ=255)
                                      cGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATtcttct                                                                                                                <  2:919265/150‑1 (MQ=255)
                                        gCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCaa                                                                                                                                     <  2:416800/127‑1 (MQ=255)
                                        gCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCaa                                                                                                                                     >  1:416800/1‑127 (MQ=255)
                                                ggCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGTGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGcc                                                                                                      <  2:347089/150‑1 (MQ=255)
                                                   cTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGCTCCTGGTGTTCGTGGTCCTCCTGCTGGCGGTGCCGCTA‑‑TCGCGGGGCAACCCGCGCCAGGGCCGCCTCCTCATTCTGCTGCCGGCAATTCTTCTGTACATGGCCTAc                                                                                                   >  2:585616/1‑150 (MQ=255)
                                                   cTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTAc                                                                                                   >  2:1401445/1‑150 (MQ=255)
                                                    tGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACAt                                                                                                          <  2:1041895/142‑1 (MQ=255)
                                                    tGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACAt                                                                                                          >  1:1041895/1‑142 (MQ=255)
                                                    tGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTCTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTAcc                                                                                                  <  1:749592/150‑1 (MQ=255)
                                                           gcgcgcCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAAt                                                                                           <  1:841010/150‑1 (MQ=255)
                                                              cgcCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCt                                                                                        <  1:1617063/150‑1 (MQ=255)
                                                              cgcCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCt                                                                                        <  1:1401445/150‑1 (MQ=255)
                                                                ccGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGa                                                                                      <  2:157558/150‑1 (MQ=255)
                                                                 cGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGAt                                                                                     >  1:1899694/1‑150 (MQ=255)
                                                                   aGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGAttt                                                                                   <  2:2379659/150‑1 (MQ=255)
                                                                          aTGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCt                                                                                                 <  1:289280/129‑1 (MQ=255)
                                                                          aTGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCt                                                                                                 >  2:289280/1‑129 (MQ=255)
                                                                           tGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTAcgc                                                                           >  2:1545659/1‑150 (MQ=255)
                                                                           tGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCAACCTGATGGCGGCCCCGCTG‑‑CCGCGGGTCAACCCTGGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTAcgc                                                                           >  1:656708/1‑150 (MQ=255)
                                                                                gTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCt                                                                                        >  1:397732/1‑132 (MQ=255)
                                                                                gTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCt                                                                                        <  2:397732/132‑1 (MQ=255)
                                                                                 ttGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGcc                                                                     >  1:1107890/1‑150 (MQ=255)
                                                                                    tCGCTGCCGATCCTGGTGTTCGTCGTCACCCTCCTCAAGGTGCGGCTGCATCGCCGG‑‑AACCCGCGCCAGGGCCTTTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTc                                                                  >  1:1022810/1‑150 (MQ=255)
                                                                                      gCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCga                                                                <  2:1022810/150‑1 (MQ=255)
                                                                                        tGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCgaga                                                              >  1:622110/1‑150 (MQ=255)
                                                                                          ccGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAg                                                            <  2:656708/150‑1 (MQ=255)
                                                                                          ccGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAg                                                            >  1:111004/1‑150 (MQ=255)
                                                                                           cGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAgg                                                           <  1:2230402/150‑1 (MQ=255)
                                                                                           cGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAgg                                                           <  1:1336432/150‑1 (MQ=255)
                                                                                             aTCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGccc                                                                    >  2:2106515/1‑139 (MQ=255)
                                                                                             aTCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGccc                                                                    <  1:2106515/139‑1 (MQ=255)
                                                                                                 tGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGc                                                                                         >  1:439668/1‑114 (MQ=255)
                                                                                                 tGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGc                                                                                         <  2:439668/114‑1 (MQ=255)
                                                                                                 tGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAAc                                                     >  2:2203739/1‑150 (MQ=255)
                                                                                                 tGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAAc                                                     >  2:2356274/1‑150 (MQ=255)
                                                                                                         gtcgtcACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCaa                                                                                                                       <  2:2255396/76‑1 (MQ=255)
                                                                                                         gtcgtcACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCaa                                                                                                                       >  1:2255396/1‑76 (MQ=255)
                                                                                                         gtcgtcACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATc                                             >  2:586832/1‑150 (MQ=255)
                                                                                                            gtcACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTAcgcg                                                                          <  1:2300057/118‑1 (MQ=255)
                                                                                                            gtcACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTAcgcg                                                                          >  2:2300057/1‑118 (MQ=255)
                                                                                                             tcACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCCGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGa                                               >  1:72418/1‑144 (MQ=255)
                                                                                                             tcACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCCGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGa                                               <  2:72418/144‑1 (MQ=255)
                                                                                                                cccTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTgg                                      <  2:1473677/150‑1 (MQ=255)
                                                                                                                cccTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTCCCGATCGCCCTgg                                      <  2:1899694/150‑1 (MQ=255)
                                                                                                                 ccTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTggg                                     >  1:1089033/1‑150 (MQ=255)
                                                                                                                 ccTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTggg                                     >  1:2229576/1‑150 (MQ=255)
                                                                                                                  ctgctgGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGc                                    >  1:579341/1‑150 (MQ=255)
                                                                                                                   tgctgGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCa                                   >  1:1570054/1‑150 (MQ=255)
                                                                                                                    gctgGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCAt                                  <  2:2229576/150‑1 (MQ=255)
                                                                                                                       ggCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAACCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCAtgtg                               >  2:1294584/1‑150 (MQ=255)
                                                                                                                         cGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATgtggt                             <  1:2203739/150‑1 (MQ=255)
                                                                                                                         cGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATgtggt                             <  1:2356274/150‑1 (MQ=255)
                                                                                                                          ggTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGTACGGCAAACTGCCGATCGCCCTGGGCATgtggtg                            <  1:1536617/150‑1 (MQ=255)
                                                                                                                           gTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGtggtgg                           <  2:622110/150‑1 (MQ=255)
                                                                                                                           gTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGtggtgg                           <  2:1604172/150‑1 (MQ=255)
                                                                                                                              ccGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGgcgc                                                                      >  2:14452/1‑104 (MQ=255)
                                                                                                                              ccGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGgcgc                                                                      <  1:14452/104‑1 (MQ=255)
                                                                                                                                gCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCa                      >  1:1809234/1‑150 (MQ=255)
                                                                                                                                  cctggcgcgGGTTGACCCGCGACAGCGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTcc                       >  2:260349/6‑150 (MQ=255)
                                                                                                                                      gcgcgGGTTGACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAAGCCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGc                  >  1:556808/2‑150 (MQ=255)
                                                                                                                                        gcgGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGa                                                                                               <  2:1331054/71‑1 (MQ=255)
                                                                                                                                        gcgGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGa                                                                                               >  1:1331054/1‑71 (MQ=255)
                                                                                                                                        gcgGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCAGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCt                >  1:2202833/1‑150 (MQ=255)
                                                                                                                                          gggTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTg                                                                                       <  2:2200166/77‑1 (MQ=255)
                                                                                                                                          gggTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTg                                                                                       >  1:2200166/1‑77 (MQ=255)
                                                                                                                                           ggTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGtt             <  1:586832/150‑1 (MQ=255)
                                                                                                                                             tCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTcc           >  2:1874551/1‑150 (MQ=255)
                                                                                                                                               aaCCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCctg         >  2:2289746/1‑150 (MQ=255)
                                                                                                                                                aCCCGCGCCAGTGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCctgc        >  2:131072/1‑150 (MQ=255)
                                                                                                                                                  ccGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATgtggt                             <  2:2307584/127‑1 (MQ=255)
                                                                                                                                                  ccGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGCCATgtggt                             >  1:2307584/1‑127 (MQ=255)
                                                                                                                                                  ccGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGCGCAAGAGGTCGGTCCACGGCCTGATCgcgctg      >  2:992269/1‑150 (MQ=255)
                                                                                                                                                      gcCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCCTGCTGATCg  <  2:111004/150‑1 (MQ=255)
                                                                                                                                                      gcCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCCTGCTGATCg  <  2:1107890/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCCTGCTGATCG  >  NC_002947/1121309‑1121606

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