Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,792,128 A→G 13.2% intergenic (+29/+39) ynfM → / ← PP_5073 efflux transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,792,1280AG13.2% 422.1 / 43.3 159intergenic (+29/+39)ynfM/PP_5073efflux transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (68/70);  new base G (8/13);  total (76/83)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 161 total aligned reads displayed.

GGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCA  >  NC_002947/5791983‑5792276
                                                                                                                                                    |                                                                                                                                                    
   gtgtTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                                               <  1:4547983/150‑1 (MQ=255)
   gtgtTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                                               <  1:5170956/150‑1 (MQ=255)
   gtgtTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                                               <  1:5170984/150‑1 (MQ=255)
      ttCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAg                                                                                                                                            <  2:1685116/150‑1 (MQ=255)
       tCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAGGTTGTGGGTGc                                                                                                                                           <  2:2209076/150‑1 (MQ=255)
         tGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCTCCCACAGGTTGTGGGTGCgg                                                                                                                                         <  1:2670427/150‑1 (MQ=255)
          ggCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCggg                                                                                                                                        <  1:1504951/150‑1 (MQ=255)
           gCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCCTTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgt                                                                                                                                       <  2:4843674/150‑1 (MQ=255)
           gCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCTCCCACAGGTTGTGGGTGCgggc                                                                                                                                       <  2:3227868/150‑2 (MQ=255)
           gCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgt                                                                                                                                       <  2:4843506/150‑1 (MQ=255)
             aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGTCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGGCCTCTTCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCggg                                                                                                                                        <  1:3777451/150‑1 (MQ=25)
             aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                     <  1:3279561/150‑1 (MQ=255)
             aCCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                     <  1:554264/150‑1 (MQ=255)
                aGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGcgtgcc                                                                                                                                  <  2:1625902/150‑3 (MQ=37)
                 gTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGAcc                                                                                                                                 <  1:2861041/150‑1 (MQ=255)
                       cTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAag                                                                                                                           <  1:5073283/150‑1 (MQ=255)
                       cTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGaat                                                                                                                           <  2:2568118/150‑2 (MQ=37)
                        tGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGTGCACGCCCGCACCCACAACCtgtg                                                                                                                                                >  2:46350/1‑128 (MQ=255)
                        tGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAaga                                                                                                                          <  1:2119870/150‑1 (MQ=255)
                         ggAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAagag                                                                                                                         <  1:1974456/150‑1 (MQ=255)
                         ggAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAagag                                                                                                                         <  2:723391/150‑1 (MQ=255)
                           aaCGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGaatgcat                                                                                                                       <  2:2691114/150‑6 (MQ=25)
                             cGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                     >  1:87289/1‑134 (MQ=255)
                             cGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                     <  2:87289/134‑1 (MQ=255)
                                gATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCGCCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGc                                                                                                                                           <  1:3705060/125‑1 (MQ=255)
                                gATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGTGCgg                                                                                                                                         >  1:1282627/1‑127 (MQ=255)
                                gATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatc                                                                                                                   <  1:142064/149‑1 (MQ=255)
                                gATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatc                                                                                                                   >  2:142064/1‑149 (MQ=255)
                                 aTCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcag                                                                                                                 >  2:1491391/1‑150 (MQ=255)
                                     ggCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTTAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagatca                                                                                                             <  1:1102423/150‑1 (MQ=255)
                                        tGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGtt                                                                                                          <  1:3932582/150‑1 (MQ=255)
                                        tGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGtt                                                                                                          >  2:5050827/1‑150 (MQ=255)
                                          ttCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtga                                                                                                        <  1:4915558/150‑1 (MQ=255)
                                             aTTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTtgatga                                                                                                     <  1:1794165/150‑1 (MQ=255)
                                              ttGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTgatgat                                                                                                    >  1:983767/1‑150 (MQ=255)
                                                ggCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtct                                                                                                  >  2:3183447/1‑150 (MQ=255)
                                                  cAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCTCCCACAGGTTGTGGGTGCgggc                                                                                                                                       <  1:1478987/111‑2 (MQ=38)
                                                  cAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCg                                                                                                >  1:332206/1‑150 (MQ=255)
                                                    gCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGAc                                                                                              <  2:150220/150‑1 (MQ=255)
                                                          tGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAAcc                                                                                                                                                    >  1:3508499/1‑90 (MQ=255)
                                                          tGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                     <  1:3864750/105‑1 (MQ=255)
                                                          tGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                                     >  2:3864750/1‑105 (MQ=255)
                                                          tGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatc                                                                                        <  2:4917008/150‑1 (MQ=255)
                                                          tGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatc                                                                                        <  1:176876/150‑1 (MQ=255)
                                                           ggCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatcc                                                                                       <  1:3101118/150‑1 (MQ=255)
                                                               gTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccacat                                                                                   >  2:649400/1‑150 (MQ=255)
                                                                    cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtc                                                                                                   >  2:1728911/1‑129 (MQ=255)
                                                                    cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtc                                                                                                   <  1:1728911/129‑1 (MQ=255)
                                                                    cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACt                                                                              >  2:1066903/1‑150 (MQ=255)
                                                                    cTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACt                                                                              >  1:1335989/1‑150 (MQ=255)
                                                                       gATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTcc                                                                           <  1:1763317/150‑1 (MQ=255)
                                                                         tggtggCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCg                                                                         >  2:2553426/1‑150 (MQ=255)
                                                                            tggCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCt                                                                      >  2:3653494/1‑150 (MQ=255)
                                                                            tggCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCt                                                                      >  1:1539776/1‑150 (MQ=255)
                                                                              gCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtca                                                                    >  2:1820004/1‑150 (MQ=255)
                                                                                gTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                                  >  2:360540/1‑150 (MQ=255)
                                                                                gTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                                  >  2:5256582/1‑150 (MQ=255)
                                                                                 ttGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatca                                                                 >  2:226219/1‑150 (MQ=255)
                                                                                     aCCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTg                                                             >  1:2126068/1‑150 (MQ=255)
                                                                                        tGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCtt                                                          >  2:1925457/1‑150 (MQ=255)
                                                                                        tGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCtt                                                          >  2:1291363/1‑150 (MQ=255)
                                                                                         gAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCttt                                                         <  2:5169560/150‑1 (MQ=255)
                                                                                            cAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCa                                                      >  2:5066208/1‑150 (MQ=255)
                                                                                              aGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAAc                                                    <  1:5066118/150‑1 (MQ=255)
                                                                                              aGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGACCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAAc                                                    >  2:3026412/1‑150 (MQ=255)
                                                                                                   ccACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCAc                                               <  1:3043856/150‑1 (MQ=255)
                                                                                                    cACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGAGCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCAcc                                              >  2:2318837/1‑150 (MQ=255)
                                                                                                       ccAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtct                                                                                                  <  1:3176083/95‑1 (MQ=255)
                                                                                                       ccAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtct                                                                                                  >  2:3176083/1‑95 (MQ=255)
                                                                                                       ccAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACaa                                                                                            >  1:1682992/1‑101 (MQ=255)
                                                                                                       ccAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGt                                           >  1:4270076/1‑150 (MQ=255)
                                                                                                       ccAAAACGGAGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCTACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGt                                           >  2:794339/1‑150 (MQ=255)
                                                                                                        cAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTc                                          <  2:2423096/150‑1 (MQ=255)
                                                                                                        cAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCAATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACATCACCGGTc                                          >  2:4420523/1‑150 (MQ=255)
                                                                                                                    aTG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTg                              >  1:5260766/1‑150 (MQ=255)
                                                                                                                    aTG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTg                              <  2:856972/150‑1 (MQ=255)
                                                                                                                      g‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTggg                            >  1:4543396/1‑149 (MQ=255)
                                                                                                                         cATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTc                           >  2:4217227/1‑150 (MQ=255)
                                                                                                                           aTTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCt                          >  2:4126673/1‑150 (MQ=255)
                                                                                                                            ttCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCtt                         <  1:2469604/150‑1 (MQ=255)
                                                                                                                                tGGTCACACCCGCTCCCACAGGTTGTGGGTGCGGGCGTGCCCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCg                     >  2:5434178/1‑150 (MQ=37)
                                                                                                                                tGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACaa                                                                                            <  2:5188373/79‑1 (MQ=255)
                                                                                                                                tGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCg                     >  1:4278141/1‑150 (MQ=255)
                                                                                                                                   gCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTg                  <  1:3653150/150‑1 (MQ=255)
                                                                                                                                    cACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGa                 <  2:946342/150‑1 (MQ=255)
                                                                                                                                     acacccGCTCCCACAGGTTGTGGGTGCGGGCGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTg                        >  2:5281959/4‑142 (MQ=38)
                                                                                                                                     acacccGCTCCCACAGGTTGTGGGTGCGGGCGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTg                        <  1:5281959/139‑1 (MQ=38)
                                                                                                                                     aCGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGaa                >  2:5087122/1‑150 (MQ=255)
                                                                                                                                       acccGCTCCCACAGGTTGTGCGTGCGGGCGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACg              >  2:381905/2‑150 (MQ=255)
                                                                                                                                         ccGCTCCCACAGGTTGTGGGTGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGgt            >  2:1058443/1‑150 (MQ=255)
                                                                                                                                         ccGCTCCCACAGGTTGTGGGTGCGGGTGTGACCACGAAGAGGCCATAAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGgt            >  2:4964351/1‑150 (MQ=255)
                                                                                                                                         ccGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGgt            >  1:2546188/1‑150 (MQ=255)
                                                                                                                                            cACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTc         <  2:3361903/150‑1 (MQ=255)
                                                                                                                                               ccACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCCGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTg                              <  1:1634738/126‑1 (MQ=255)
                                                                                                                                               ccACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTc                                          <  1:3240446/114‑1 (MQ=255)
                                                                                                                                               ccACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTg                              >  2:1635057/1‑126 (MQ=255)
                                                                                                                                               ccACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTg                              >  2:1634738/1‑126 (MQ=255)
                                                                                                                                               ccACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTg                              <  1:1634699/126‑1 (MQ=255)
                                                                                                                                               ccACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTcggc      >  2:1253110/1‑150 (MQ=255)
                                                                                                                                                   aaCCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCa  >  1:4792694/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATG‑‑CATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAACGGTGTCGGCGGCA  >  NC_002947/5791983‑5792276

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