Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,121,449 C→T 11.4% V295V (GTC→GTT 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,4490CT11.4% 245.9 / 19.0 97V295V (GTC→GTTPP_0982YjgP/YjgQ family permease
Reads supporting (aligned to +/- strand):  ref base C (39/46);  new base T (5/6);  total (45/52)
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.64e-01

GCCGAGGAAGTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCC  >  NC_002947/1121300‑1121597
                                                                                                                                                       |                                                                                                                                                     
gCCGAGGAAGTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTc                                                                                                                                                       >  1:621340/1‑150 (MQ=255)
  cGAGGAAGTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCaa                                                                                                                                                     <  2:220261/150‑1 (MQ=255)
   gAGGAAGTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCTCTG‑‑TCGCGGGTCAAc                                                                                                                                                    >  2:121590/1‑150 (MQ=255)
       aaGTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTTGACCcgc                                                                                                                                                <  2:860098/150‑1 (MQ=255)
       aaGTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTTGACCcgc                                                                                                                                                <  2:1551416/150‑1 (MQ=255)
         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                                                                                                                                    >  1:583080/1‑140 (MQ=255)
                     gAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCCTG‑‑GCGCGGGTTGACCCGCGACAgcggcacc                                                                                                                                    <  2:583080/147‑8 (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:1408544/150‑1 (MQ=255)
                          gcTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCaa                                                                                                                            <  2:621340/150‑1 (MQ=255)
                           cTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAg                                                                                                                           <  2:57860/150‑1 (MQ=255)
                           cTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAg                                                                                                                           >  2:1446325/1‑150 (MQ=255)
                            ttCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgc                                                                                                                          >  2:290039/1‑150 (MQ=255)
                               ccACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgctgc                                                                                                                       >  1:157558/1‑150 (MQ=255)
                               ccACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAgctgc                                                                                                                       <  1:1316668/150‑1 (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                                                                                                                            >  1:416800/1‑127 (MQ=255)
                                                 gCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCaa                                                                                                                            <  2:416800/127‑1 (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                                                                                                 >  1:1041895/1‑142 (MQ=255)
                                                             tGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACAt                                                                                                 <  2:1041895/142‑1 (MQ=255)
                                                             tGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTCTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTAcc                                                                                         <  1:749592/150‑1 (MQ=255)
                                                                    gcgcgcCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAAt                                                                                  <  1:841010/150‑1 (MQ=255)
                                                                       cgcCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCt                                                                               <  1:1401445/150‑1 (MQ=255)
                                                                       cgcCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCt                                                                               <  1:1617063/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)
                                                                              cTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTTGACCcgcga                                                                                                                                              >  1:786613/1‑80 (MQ=255)
                                                                              cTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTTGACCcgcga                                                                                                                                              <  2:786613/81‑2 (MQ=38)
                                                                                   aTGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCt                                                                                        >  2:289280/1‑129 (MQ=255)
                                                                                   aTGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCt                                                                                        <  1:289280/129‑1 (MQ=255)
                                                                                    tGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTAcgc                                                                  >  2:1545659/1‑150 (MQ=255)
                                                                                    tGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCAACCTGATGGCGGCCCCGCTG‑‑CCGCGGGTCAACCCTGGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTAcgc                                                                  >  1:656708/1‑150 (MQ=255)
                                                                                         gTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCt                                                                               <  2:397732/132‑1 (MQ=255)
                                                                                         gTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCt                                                                               >  1:397732/1‑132 (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:1336432/150‑1 (MQ=255)
                                                                                                    cGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAgg                                                  <  1:2230402/150‑1 (MQ=255)
                                                                                                      aTCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGccc                                                           <  1:2106515/139‑1 (MQ=255)
                                                                                                      aTCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGccc                                                           >  2:2106515/1‑139 (MQ=255)
                                                                                                          tGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGc                                                                                >  1:439668/1‑114 (MQ=255)
                                                                                                          tGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGc                                                                                <  2:439668/114‑1 (MQ=255)
                                                                                                          tGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAAc                                            >  2:2356274/1‑150 (MQ=255)
                                                                                                          tGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAAc                                            >  2:2203739/1‑150 (MQ=255)
                                                                                                                  gtcgtcACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCaa                                                                                                              >  1:2255396/1‑76 (MQ=255)
                                                                                                                  gtcgtcACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCaa                                                                                                              <  2:2255396/76‑1 (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                                      <  2:72418/144‑1 (MQ=255)
                                                                                                                      tcACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCCGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGa                                      >  1:72418/1‑144 (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:2356274/150‑1 (MQ=255)
                                                                                                                                  cGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATgtggt                    <  1:2203739/150‑1 (MQ=255)
                                                                                                                                   ggTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGTACGGCAAACTGCCGATCGCCCTGGGCATgtggtg                   <  1:1536617/150‑1 (MQ=255)
                                                                                                                                    gTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGtggtgg                  <  2:1604172/150‑1 (MQ=255)
                                                                                                                                    gTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGtggtgg                  <  2:622110/150‑1 (MQ=255)
                                                                                                                                       ccGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGgcgc                                                             <  1:14452/104‑1 (MQ=255)
                                                                                                                                       ccGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGgcgc                                                             >  2:14452/1‑104 (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                                                                                      >  1:1331054/1‑71 (MQ=255)
                                                                                                                                                 gcgGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGa                                                                                      <  2:1331054/71‑1 (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)
                                                                                                                                                       |                                                                                                                                                     
GCCGAGGAAGTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCC  >  NC_002947/1121300‑1121597

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: