Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,219 T→A 6.0% Y283N (TAT→AAT)  PP_4664 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,293,2190TA 6.0% 605.3 / 18.3 181Y283N (TAT→AAT) PP_4664LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (81/89);  new base A (5/6);  total (86/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

GCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATGGTGAC  >  NC_002947/5293072‑5293363
                                                                                                                                                   |                                                                                                                                                
gcgggcggGGGCTGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTAt                                                                                                                                                <  2:4959929/150‑1 (MQ=255)
gcgggcggGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTAt                                                                                                                                                <  1:5413974/150‑1 (MQ=255)
   ggcgggGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGcc                                                                                                                                             <  1:1343913/150‑1 (MQ=255)
   ggcgggGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGcc                                                                                                                                             <  1:1758934/150‑1 (MQ=255)
     cgggGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCGc                                                                                                                                           <  2:2269011/150‑1 (MQ=255)
      ggggCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCa                                                                                                                                          <  1:723595/150‑1 (MQ=255)
        ggCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGa                                                                                                                                        >  1:619591/1‑150 (MQ=255)
        ggCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCGCaaa                                                                                                                                        <  1:344887/150‑3 (MQ=255)
         gCCGGAATCGGGGTGGTGCCACTGGCAGCAGTCCAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGaa                                                                                                                                       >  2:5031580/1‑150 (MQ=255)
         gCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGaa                                                                                                                                       <  2:89345/150‑1 (MQ=255)
          ccGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAAc                                                                                                                                      <  2:2275627/150‑1 (MQ=255)
             gAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCGCaaaatccc                                                                                                                                   <  1:3681056/150‑5 (MQ=255)
               aTCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTa                                                                                                                                 <  2:3578041/150‑1 (MQ=255)
                     gtggtgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCAcc                                                                                                                           <  2:4364209/150‑1 (MQ=255)
                         tgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTa                                                                                                                       <  2:3676011/150‑1 (MQ=255)
                         tgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTa                                                                                                                       >  2:4140886/1‑150 (MQ=255)
                         tgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTa                                                                                                                       >  2:4432402/1‑150 (MQ=255)
                         tgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTa                                                                                                                       >  2:4432407/1‑150 (MQ=255)
                           ccACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTAcc                                                                                                                     <  1:5033483/150‑1 (MQ=255)
                            cACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCg                                                                                                                    <  1:2650734/150‑1 (MQ=255)
                             aCTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCgg                                                                                                                   >  1:2155267/1‑150 (MQ=255)
                              cTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGc                                                                                                                  <  1:680295/150‑1 (MQ=255)
                               tGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCt                                                                                                                 >  1:755395/1‑150 (MQ=255)
                                 gcagcaGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGa                                                                                                               <  1:1959837/150‑1 (MQ=255)
                                 gcagcaGTACAGCGCTGGCAAGGCAAGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGa                                                                                                               <  2:2831874/150‑1 (MQ=255)
                                   agcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAcg                                                                                                             >  1:1158151/1‑150 (MQ=255)
                                    gcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAcgc                                                                                                            <  2:2511167/150‑1 (MQ=255)
                                    gcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAcgc                                                                                                            <  2:2669397/150‑1 (MQ=255)
                                         aCAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGACCTGGTAGCCACCTGTACCGGCTGTCGCGCCAg                                                                                                       >  2:4341677/1‑150 (MQ=255)
                                         aCAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAg                                                                                                       <  1:4140886/150‑1 (MQ=255)
                                                tGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGcac                                                                                                >  1:4589056/1‑150 (MQ=255)
                                                tGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGcac                                                                                                >  2:4811922/1‑150 (MQ=255)
                                                    aaGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCAcacca                                                                                            >  1:2243791/1‑150 (MQ=255)
                                                    aaGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCAcacca                                                                                            >  2:422887/1‑150 (MQ=255)
                                                     aGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCAcaccac                                                                                           >  2:4997203/1‑150 (MQ=255)
                                                       gCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTg                                                                                         >  2:3405563/1‑150 (MQ=255)
                                                        cATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGt                                                                                        <  1:4341677/150‑1 (MQ=255)
                                                           gctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                                     <  1:1606765/150‑1 (MQ=255)
                                                           gctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                                     <  1:2390183/150‑1 (MQ=255)
                                                           gctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                                     <  2:4324607/150‑1 (MQ=255)
                                                           gctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                                     <  1:5409031/150‑1 (MQ=255)
                                                           gctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                                     <  2:1190888/150‑1 (MQ=255)
                                                            ctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAg                                                                                    <  1:4165301/150‑1 (MQ=255)
                                                               ccGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCa                                                                                 >  2:5193571/1‑150 (MQ=255)
                                                               ccGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCa                                                                                 >  1:3783359/1‑150 (MQ=255)
                                                                        tGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTg                                                                        >  1:763234/1‑150 (MQ=255)
                                                                         gggTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGc                                                                       <  1:2230632/150‑1 (MQ=255)
                                                                          ggTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGt                                                                                                                                  >  2:689666/1‑90 (MQ=255)
                                                                           gTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCa                                                                     <  2:570267/150‑1 (MQ=255)
                                                                             ggCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATa                                                                   <  2:185498/150‑1 (MQ=255)
                                                                              gCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAgc                                                                 <  2:5338491/150‑1 (MQ=255)
                                                                              gCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGACATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAg                                                                  <  1:4524189/150‑1 (MQ=255)
                                                                               cACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAgc                                                                 >  1:1655289/1‑150 (MQ=255)
                                                                                aCTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAgcg                                                                >  2:3021087/1‑150 (MQ=255)
                                                                                aCTGCAGGAGCCAGGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGAGTTTGCGGCATTGCCGGGCTCTGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCACACCGCCTGCCATAgcg                                                                >  2:3854375/1‑150 (MQ=255)
                                                                                  tGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGcgcg                                                              <  1:3880873/150‑1 (MQ=255)
                                                                                   gCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCgg                                                             <  1:1623437/150‑1 (MQ=255)
                                                                                    cAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCggg                                                            <  1:4432402/150‑1 (MQ=255)
                                                                                      ggAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCa                                                          >  1:5331307/1‑150 (MQ=255)
                                                                                        aGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGa                                                        <  2:182235/150‑1 (MQ=255)
                                                                                        aGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGa                                                        <  1:3225970/150‑1 (MQ=255)
                                                                                          ccATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGaca                                                      <  2:5364732/150‑1 (MQ=255)
                                                                                           cATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGacac                                                     <  1:4806100/150‑1 (MQ=255)
                                                                                           cATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGacac                                                     >  1:2325630/1‑150 (MQ=255)
                                                                                              gggCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTg                                                  <  1:328661/150‑1 (MQ=255)
                                                                                               ggCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTgg                                                 <  2:2082110/150‑1 (MQ=255)
                                                                                               ggCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTgg                                                 >  2:5107027/1‑150 (MQ=255)
                                                                                               ggCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTgg                                                 <  1:3021087/150‑1 (MQ=255)
                                                                                                 ccAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCTGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGag                                               >  2:4647687/1‑149 (MQ=255)
                                                                                                 ccAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGAc                                               >  1:3233438/1‑150 (MQ=255)
                                                                                                  cAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCGCa                                                                                                                                          >  1:425129/1‑58 (MQ=25)
                                                                                                  cAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCGCa                                                                                                                                          <  2:425129/58‑1 (MQ=25)
                                                                                                   aaCCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGt                                             >  1:2417735/1‑150 (MQ=255)
                                                                                                    aCCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCggg                                                            >  2:3173545/1‑134 (MQ=255)
                                                                                                    aCCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCggg                                                            <  1:3173545/134‑1 (MQ=255)
                                                                                                       gccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTAc                                                                                                                      <  2:4079173/73‑1 (MQ=255)
                                                                                                       gccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTAc                                                                                                                      >  1:4079173/1‑73 (MQ=255)
                                                                                                       gccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAg                                         <  2:5434325/150‑1 (MQ=255)
                                                                                                           gTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCa                                     <  1:3651606/150‑1 (MQ=255)
                                                                                                            ttGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAg                                    <  2:537155/150‑1 (MQ=255)
                                                                                                            ttGCTGGTGTGCGCACGGGATTTTGCGGCAGTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAg                                    <  2:5238742/150‑1 (MQ=255)
                                                                                                                tGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCac                                >  2:4735490/1‑149 (MQ=255)
                                                                                                                  gtgtGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTGGCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAAcg                              >  2:2640951/1‑149 (MQ=255)
                                                                                                                     tgCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcg                           >  2:144947/1‑150 (MQ=255)
                                                                                                                       cgcACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcggc                         <  1:3136121/150‑1 (MQ=255)
                                                                                                                          aCGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcggcgcg                      >  2:3109278/1‑150 (MQ=255)
                                                                                                                          aCGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcggcgcg                      >  2:809531/1‑150 (MQ=255)
                                                                                                                          aCGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcggcgcg                      >  1:1426173/1‑150 (MQ=255)
                                                                                                                          aCGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcggcgcg                      <  2:3281718/150‑1 (MQ=255)
                                                                                                                           cGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGc                                  <  2:2197056/137‑1 (MQ=255)
                                                                                                                              gATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCg                  >  1:5114269/1‑150 (MQ=255)
                                                                                                                              agttctGCGGCATAGCCCGGCAATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCg                  >  1:291383/6‑150 (MQ=255)
                                                                                                                                 tttGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTg               <  1:1128304/150‑1 (MQ=255)
                                                                                                                                  ttGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTgg              <  2:685579/150‑1 (MQ=255)
                                                                                                                                  ttGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGCCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTgg              <  1:4582687/150‑1 (MQ=255)
                                                                                                                                    gCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGgcg                                                                                                   >  1:2625443/1‑63 (MQ=255)
                                                                                                                                    gCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGc            <  2:831459/150‑1 (MQ=255)
                                                                                                                                      ggCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCt          >  2:1354050/1‑150 (MQ=255)
                                                                                                                                      ggCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCt          >  1:855562/1‑150 (MQ=255)
                                                                                                                                           tGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATGGt     <  2:755395/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATGGTGAC  >  NC_002947/5293072‑5293363

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: