Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,093 T→C 21.1% intergenic (+63/‑225) PP_0367 → / → PP_0368 hypothetical protein/acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947445,0930TC21.1% 297.7 / 64.4 143intergenic (+63/‑225)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (28/80);  new base C (4/26);  total (36/107)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 146 total aligned reads displayed.

GTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTA  >  NC_002947/444944‑445242
                                                                                                                                                     |                                                                                                                                                     
gTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGt                                                                                                                                                       <  1:2462938/150‑1 (MQ=255)
gTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAAGGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGt                                                                                                                                                       <  2:780840/150‑1 (MQ=255)
 tGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGTGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTc                                                                                                                                                      <  1:3278018/150‑1 (MQ=255)
  gggCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTcc                                                                                                                                                     <  1:1572509/150‑1 (MQ=255)
  gggCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTcc                                                                                                                                                     <  2:45007/150‑1 (MQ=255)
  gggCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTcc                                                                                                                                                     <  2:3322057/150‑1 (MQ=255)
  gggCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTcc                                                                                                                                                     <  2:3322009/150‑1 (MQ=255)
  gggCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTcc                                                                                                                                                     <  2:3217374/150‑1 (MQ=255)
    gCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGa                                                                                                                                                   <  2:936718/150‑1 (MQ=255)
    gCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGa                                                                                                                                                   <  1:4785389/150‑1 (MQ=255)
     cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGGCCCCGGATGCGGCCGAt                                                                                                                                                  >  1:6799500/1‑150 (MQ=255)
     cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                                  <  1:1125610/150‑1 (MQ=255)
     cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                                  <  2:5073585/150‑1 (MQ=255)
     cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                                  <  2:2609560/150‑1 (MQ=255)
     cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                                  <  1:4569745/150‑1 (MQ=255)
     cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                                  <  1:2566319/150‑1 (MQ=255)
      aaGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                                 <  2:5349945/150‑1 (MQ=255)
      aaGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                                 <  1:5218090/150‑1 (MQ=255)
       aGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa                                                                                                                                                <  2:5558380/150‑1 (MQ=255)
         aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                              <  2:479057/150‑1 (MQ=255)
         aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                              <  2:596451/150‑1 (MQ=255)
         aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                              <  1:908697/150‑1 (MQ=255)
         aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                              <  1:6872123/150‑1 (MQ=255)
         aCCGAAGTCGAGCTCGGCGCAAAGGTCGACAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                              <  1:5695215/150‑1 (MQ=255)
             aaGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCGATCAGacaac                                                                                                                                          >  2:1122191/1‑150 (MQ=255)
                 cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                      <  1:6656955/150‑1 (MQ=255)
                 cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                      <  1:3339792/150‑1 (MQ=255)
                   aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                    <  2:2650154/150‑1 (MQ=255)
                   aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                    <  1:1802119/150‑1 (MQ=255)
                   aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                    <  2:640288/150‑1 (MQ=255)
                   aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                    <  1:220968/150‑1 (MQ=255)
                   aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                    <  2:2840041/150‑1 (MQ=255)
                   aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                    <  1:3973241/150‑1 (MQ=255)
                     cTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGACGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCCGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATtgttgt                                                                                                                                  <  2:1616365/150‑1 (MQ=25)
                      tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                 <  1:1716499/150‑1 (MQ=255)
                      tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                 <  1:2364201/150‑1 (MQ=255)
                       cGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                                <  2:3917913/150‑1 (MQ=255)
                           gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAAAATTGTTGTCTGATc                                                                                                                            <  1:4917800/150‑1 (MQ=38)
                           gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTGATc                                                                                                                            <  2:5826886/150‑1 (MQ=38)
                             aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                          <  2:2979074/150‑1 (MQ=255)
                               aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTGATCggac                                                                                                                        <  1:5819972/150‑3 (MQ=37)
                               aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTGATCggac                                                                                                                        <  2:6248864/150‑3 (MQ=37)
                                       aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgc                                                                                                                <  1:6328727/150‑1 (MQ=255)
                                        aGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                               <  1:5579661/150‑1 (MQ=255)
                                        aGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                               <  2:4315992/150‑1 (MQ=255)
                                         ggCATCCGACGACAGGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                              <  2:5728282/150‑1 (MQ=17)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGGGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                              <  1:3226214/150‑1 (MQ=18)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGGGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                              <  1:6353561/150‑1 (MQ=18)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                              <  2:2212471/150‑1 (MQ=255)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                              <  1:987461/150‑1 (MQ=255)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                              <  2:3344718/150‑1 (MQ=255)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                              <  1:6089021/150‑1 (MQ=255)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                              <  1:2231658/150‑1 (MQ=255)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                              <  2:4642513/150‑1 (MQ=255)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                              <  2:5280060/150‑1 (MQ=18)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                              <  1:3887943/150‑1 (MQ=18)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                              <  1:2265969/150‑1 (MQ=18)
                                         ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                              <  2:6789821/150‑1 (MQ=18)
                                         ggCATCCGACGAAATACTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                              <  2:4665078/150‑1 (MQ=18)
                                             tCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCaa                                                                                                          <  1:19058/150‑1 (MQ=255)
                                               cgacgaCATGCTGGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATGTCCGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                        <  2:302151/150‑1 (MQ=255)
                                               cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGATCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                        >  2:4314774/1‑150 (MQ=255)
                                               cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGc                                                                                                        <  1:6597807/150‑1 (MQ=17)
                                                gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTTCCCGCGAATTCGTGCGATAAGACAAAAAATTTTTGTTGTGGGGGCTCTTAGCGGGGAAGcc                                                                                                       >  2:4646372/1‑150 (MQ=255)
                                                gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGcc                                                                                                       <  1:1272922/150‑1 (MQ=17)
                                                gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGCAGGAGCGGGCGTGCCCGCGGATGCGTCCGATGAGACAACAATTTTGTATGGTGGGGGCGCGTGGCGGGCAAGcc                                                                                                       >  2:3342643/1‑150 (MQ=255)
                                                 acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTCCCCCCACAGCGGCCCATCAAACAACAAATTTGTGTGGTGGGGCCGCGTGGCGGGCCAcccc                                                                                                      >  2:4128693/1‑150 (MQ=255)
                                                 acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCGCAAGCGTCCGATCAGACAACAAATTTGTCTGGTCGGGGCGCGTCGCGGGCAAGccc                                                                                                      >  2:1562800/1‑150 (MQ=255)
                                                   gaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTTCCCCCGCAAGCGGGCCGTCCGACAACCAATTTGTGTGGGGGGGGCGCGTAGCCGGCCAGccccc                                                                                                    >  2:2416704/1‑148 (MQ=255)
                                                   gaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGc                                                                                                    <  1:3342489/150‑1 (MQ=255)
                                                     aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTTCCCCCGAATGCGTCCGATAAGAAAACAATTTTGTGTGGTAGGGCCGCTTAGCGGGCAAGcccccgc                                                                                                  >  2:5062924/1‑146 (MQ=255)
                                                     aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCGATAAGACAACAAATTTGTCTGGTCGGGCCTCGTCGCGGGCAAGCCCGcgc                                                                                                  >  2:5614083/1‑148 (MQ=255)
                                                          ttGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCCCGAAAGCGTCCGAACACACAACAATTTTGTATGGTAGGGCCTCGACGCGGGCAAGCCCCCTCCAACa                                                                                             >  2:6200567/1‑150 (MQ=255)
                                                                   tcatcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCCGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGAGGGCAAGCCCGCTCCTACAAGGGTTtgt                                                                                    <  1:6475090/150‑1 (MQ=255)
                                                                    aatcGACACCCACGCCTCAATACTTGCGACACCCCCTATTCCAGCGCCTGCCTTGTTGGCGCGGGCTTGCCCGCTAATGCGTCCGATCACAAAAAATTTTTGTCTGGTCGGGCCTCTGCGCGGGCACGCCCGCTCCTACAAGGGTTtgtg                                                                                   <  2:1007000/149‑1 (MQ=255)
                                                                             cTACGCCTTATTACTTGCTTCATCCCCTATTTCTGCGCCGGCTTTGTTGGTGGGGGATTGCCCGCGATTGCGTCCTATCATAAAACATTTTTGTCTGGGGGGGCTTCTTCGGGGGCATGCCCGCTCCTACAAGGGTTTGTggcatggcat                                                                          <  2:1095043/150‑1 (MQ=255)
                                                                             cTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCCCGAATTCGTGCGATAAGACAACAAATTTTTGTGGTCGGGGCTCTTCGCGGGCAAGCCCGCACCAACAAGGGTTTGTggcatggcat                                                                          >  2:6154622/1‑150 (MQ=255)
                                                                              tACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTTCCCGCGAAATCCTCCCCACAGACAACAAATTTTTGTGGGGGGGGCGCGTGGGGGGGAAGCCCCCGCCAAAAACGGGTTGTGggaatgcatg                                                                         >  2:1551938/1‑150 (MQ=255)
                                                                              cccGCCTTATTCCTGGCGTCAGCCCCTATTCCCGCGCCGGCCTTGTTGGTGGGGGGGTGCCCGCTAATGCGTCCCATCACACAACAATTTTGTGTGGTGGGGCCTCTTGGCGGGAACGCCCGCTCCTACAAGGGTTTGTGgcatggcatg                                                                         <  2:915010/148‑1 (MQ=255)
                                                                                      aaTACTTGCTTCTTCCCCTTTTTCCGCGCCTGCCTTGTAGGTGCGGGCTGGCCCGCGATTGCGTCCGATCAGACAACAATTTTGTCTGGTGGGGCGGCTTCGCGGGCAAGCCAGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                                 <  2:6548226/150‑1 (MQ=255)
                                                                                                     cccTCTTTCCGCGCCTGCCTTGTCGGCGCGGGCTTGCCCGCGAATGCGGCCGATCTGACACCAATTTTGTCTGGGCGGGCCTGTCCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                                  <  2:4450460/150‑1 (MQ=255)
                                                                                                     cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATATTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                        >  2:1901925/1‑96 (MQ=255)
                                                                                                     cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCCATCAGACAACAAGTTTTTCTGGTCGGGCCACTTAGCGGGCAAGCCCCCACCAACAAGGGTTTGTGGCATGGCAATACTGATTTTCCAATTTTCGTGGc                                                  >  2:6155046/1‑150 (MQ=255)
                                                                                                         aTTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                          >  1:1037009/1‑74 (MQ=255)
                                                                                                               gcgcCGGCCGTGTTGGCGGGGGCTTGCCCGCGAATGCGTCCGACCAGACACAAATTTTGTCGGGTGGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCa                                        <  2:3782753/150‑1 (MQ=38)
                                                                                                                      ctttGTCGGTGGGGGTTGGCCCGCTAATGCGTCCTAAAACAAAAAATTTGTGTGGGGGCGGGCCTCTTGGCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                                 <  1:5025796/148‑1 (MQ=255)
                                                                                                                           ttgttgggggTTTGCCCGCTATTTCGTCCTTTTTTACAAAATTTTTGTGTGGTGGGGCTTCTGGGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATc                            <  2:3039172/143‑1 (MQ=255)
                                                                                                                             ggAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc                                                                                                                            <  1:1673956/52‑1 (MQ=11)
                                                                                                                             ggAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc                                                                                                                            >  2:1673956/1‑52 (MQ=11)
                                                                                                                               aGCGGGCTTGCCCGCGACTGGGTCCGATCAGACAACAATTTGGTGTGGTCGGGCCTCTTCGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTg                        <  2:4414413/150‑1 (MQ=255)
                                                                                                                                      ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                 >  1:2404366/1‑150 (MQ=17)
                                                                                                                                        gCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc               >  2:4254679/1‑150 (MQ=255)
                                                                                                                                         cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGGCGGGCCTCTTCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt              <  2:103982/150‑1 (MQ=255)
                                                                                                                                         cccGCCATTGCGCCCTAACACACAACAATTTGGTCTGGGCGGCCCTCTGCGCGGCCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt              <  2:104155/150‑1 (MQ=255)
                                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg             >  1:4788429/1‑150 (MQ=255)
                                                                                                                                          ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg             >  1:2768266/1‑150 (MQ=255)
                                                                                                                                          acgcgAATGAGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg             >  1:1544142/2‑150 (MQ=255)
                                                                                                                                           cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa            >  1:914217/1‑150 (MQ=255)
                                                                                                                                           cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa            >  1:5741130/1‑150 (MQ=255)
                                                                                                                                                     tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa  >  1:1822345/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTA  >  NC_002947/444944‑445242

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