Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,121,450 A→G 11.5% N296D (AAC→GAC)  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,4500AG11.5% 245.4 / 18.3 96N296D (AAC→GAC) PP_0982YjgP/YjgQ family permease
Reads supporting (aligned to +/- strand):  ref base A (39/45);  new base G (5/6);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

CGAGGAAGTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCCTG  >  NC_002947/1121302‑1121599
                                                                                                                                                      |                                                                                                                                                      
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: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: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:2203739/1‑150 (MQ=255)
                                                                                                        tGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAAc                                              >  2:2356274/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)
                                                                                                                                                      aaCCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCctg  >  2:2289746/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CGAGGAAGTGACCGAGCGCGAAGCGCTTCCCACCTCCGAACTGATCGGCAAGGAGGGCCTGCGTGAGCGCGCCGAGCTGCAATGGCGGTTGTCGCTGCCGATCCTGGTGTTCGTCGTCACCCTGCTGGCGGTGCCGCTG‑‑TCGCGGGTCAACCCGCGCCAGGGCCGCTTCCTCAAGCTGCTGCCGGCAATTCTTCTGTACATGGCCTACCTGACAATGCTGATTTCCGTACGCGGCGCCCTCGAGAAGGGCAAACTGCCGATCGCCCTGGGCATGTGGTGGGTCCACGGCCTGTTCCTG  >  NC_002947/1121302‑1121599

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: