Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,090 G→A 14.3% intergenic (+60/‑228) 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,0900GA14.3% 362.9 / 45.7 149intergenic (+60/‑228)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (27/94);  new base A (4/17);  total (37/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

Only 100 of 149 total aligned reads displayed.

CGGTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGA  >  NC_002947/444942‑445232
                                                                                                                                                    |                                                                                                                                              
cGGTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAACTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGc                                                                                                                                               <  1:2194623/150‑1 (MQ=255)
cGGTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGc                                                                                                                                               <  2:1851584/150‑1 (MQ=255)
  gTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAAGGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGt                                                                                                                                             <  2:780840/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:3217374/150‑1 (MQ=255)
      gCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGa                                                                                                                                         <  2:936718/150‑1 (MQ=255)
      gCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGa                                                                                                                                         <  1:566598/150‑1 (MQ=255)
      gCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGa                                                                                                                                         <  1:3115066/150‑1 (MQ=255)
       cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGGCCCCGGATGCGGCCGAt                                                                                                                                        >  1:6799500/1‑150 (MQ=255)
       cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCGCGCATGCgggcgct                                                                                                                                        >  2:6478185/1‑144 (MQ=255)
       cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                        <  1:1125610/150‑1 (MQ=255)
       cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                        <  2:2609560/150‑1 (MQ=255)
       cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                        <  1:4569745/150‑1 (MQ=255)
       cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                        <  1:4210284/150‑1 (MQ=255)
       cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                        <  1:2566319/150‑1 (MQ=255)
        aaGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGGCCGCCCCTGCGTCCGCTc                                                                                                                                       >  2:4305729/1‑150 (MQ=255)
         aGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa                                                                                                                                      <  2:5558380/150‑1 (MQ=255)
           aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                    <  1:908697/150‑1 (MQ=255)
           aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                    <  2:4387337/150‑1 (MQ=255)
           aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                    <  2:596451/150‑1 (MQ=255)
           aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                    <  2:5975970/150‑1 (MQ=255)
           aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                    <  1:6872123/150‑1 (MQ=255)
               aaGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCGATCAGacaac                                                                                                                                >  2:1122191/1‑150 (MQ=255)
                aGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGacaaca                                                                                                                               <  2:4003263/150‑1 (MQ=255)
                   cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                            <  1:3339792/150‑1 (MQ=255)
                   cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                            <  1:6656955/150‑1 (MQ=255)
                     aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                          <  1:220968/150‑1 (MQ=255)
                     aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                          <  2:2840041/150‑1 (MQ=255)
                     aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                          <  1:1802119/150‑1 (MQ=255)
                     aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                          <  2:4888161/150‑1 (MQ=255)
                       cTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGACGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCCGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATtgttgt                                                                                                                        <  2:1616365/150‑1 (MQ=25)
                        tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                       <  1:1716499/150‑1 (MQ=255)
                         cGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCt                                                                                                                      <  2:3917913/150‑1 (MQ=255)
                             gcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTc                                                                                                                  <  1:5745383/150‑1 (MQ=255)
                                 aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                              <  2:3597597/150‑1 (MQ=255)
                                 aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTGATCggac                                                                                                              <  2:6248864/150‑3 (MQ=37)
                                    tcgtcAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCctc                                                                                                           <  1:3335383/150‑1 (MQ=255)
                                          aGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                     <  2:4315992/150‑1 (MQ=255)
                                          aGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgcg                                                                                                     <  1:5579661/150‑1 (MQ=255)
                                           ggCATCCGACGACAGGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  2:5728282/150‑1 (MQ=17)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCGTCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  2:6197076/150‑1 (MQ=18)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGGGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  1:6353561/150‑1 (MQ=18)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGGGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  1:3226214/150‑1 (MQ=18)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                    <  2:2212471/150‑1 (MQ=255)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                    <  1:2231658/150‑1 (MQ=255)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                    <  1:3433796/150‑1 (MQ=255)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                    <  2:4642513/150‑1 (MQ=255)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                    <  1:987461/150‑1 (MQ=255)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                    <  1:6089021/150‑1 (MQ=255)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                    <  2:3344718/150‑1 (MQ=255)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  2:6308824/150‑1 (MQ=21)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  1:3887943/150‑1 (MQ=18)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  2:3327848/150‑1 (MQ=18)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  2:5280060/150‑1 (MQ=18)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  1:6434391/150‑1 (MQ=18)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  1:2265969/150‑1 (MQ=18)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  1:4899146/150‑1 (MQ=18)
                                           ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  2:6789821/150‑1 (MQ=18)
                                           ggCATCCGACGAAATCCTTGGCCGGCGCATCGACACCTCCGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  2:4664967/150‑1 (MQ=17)
                                           ggCATCCGACGAAATACTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                    <  2:4665078/150‑1 (MQ=18)
                                                 cgacgaCATGCTGGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATGTCCGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                              <  2:302151/150‑1 (MQ=255)
                                                 cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGc                                                                                              <  1:6597807/150‑1 (MQ=17)
                                                  gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGcc                                                                                             <  1:1272922/150‑1 (MQ=17)
                                                   acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCGCAAGCGTCCGATCAGACAACAAATTTGTCTGGTCGGGGCGCGTCGCGGGCAAGccc                                                                                            >  2:1562800/1‑150 (MQ=255)
                                                       aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTTCCCCCGAATGCGTCCGATAAGAAAACAATTTTGTGTGGTAGGGCCGCTTAGCGGGCAAGcccccgc                                                                                        >  2:5062924/1‑146 (MQ=255)
                                                                      aatcGACACCCACGCCTCAATACTTGCGACACCCCCTATTCCAGCGCCTGCCTTGTTGGCGCGGGCTTGCCCGCTAATGCGTCCGATCACAAAAAATTTTTGTCTGGTCGGGCCTCTGCGCGGGCACGCCCGCTCCTACAAGGGTTtgtg                                                                         <  2:1007000/149‑1 (MQ=255)
                                                                               cTACGCCTTATTACTTGCTTCATCCCCTATTTCTGCGCCGGCTTTGTTGGTGGGGGATTGCCCGCGATTGCGTCCTATCATAAAACATTTTTGTCTGGGGGGGCTTCTTCGGGGGCATGCCCGCTCCTACAAGGGTTTGTggcatggcat                                                                <  2:1095043/150‑1 (MQ=255)
                                                                                tACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTTCCCGCGAAATCCTCCCCACAGACAACAAATTTTTGTGGGGGGGGCGCGTGGGGGGGAAGCCCCCGCCAAAAACGGGTTGTGggaatgcatg                                                               >  2:1551938/1‑150 (MQ=255)
                                                                                cccGCCTTATTCCTGGCGTCAGCCCCTATTCCCGCGCCGGCCTTGTTGGTGGGGGGGTGCCCGCTAATGCGTCCCATCACACAACAATTTTGTGTGGTGGGGCCTCTTGGCGGGAACGCCCGCTCCTACAAGGGTTTGTGgcatggcatg                                                               <  2:915010/148‑1 (MQ=255)
                                                                                       taaTCCTTGCTCCACCCCCTTTTCCCGCGCGTGCTTTGTTGGCGGGGGCTCGCCGTCTAATGCGCCCTAAAAAAAAAAATTTTTGTGTGGGCGGCCCTCTTCGCGGGAACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAtt                                                        <  2:1971625/149‑1 (MQ=255)
                                                                                        aaTACTTGCTTCTTCCCCTTTTTCCGCGCCTGCCTTGTAGGTGCGGGCTGGCCCGCGATTGCGTCCGATCAGACAACAATTTTGTCTGGTGGGGCGGCTTCGCGGGCAAGCCAGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                       <  2:6548226/150‑1 (MQ=255)
                                                                                        aaTACTTGCGTCACCCCCAATTCCCGCGCCGGCCTTGTTGGCGGGGGCGGGCCCGCGAATGCGCCCGAACACAAAACATTTTTGGGGGGGGGGCCCTCTCCGGGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAttt                                                       <  2:2495574/150‑1 (MQ=255)
                                                                                                       ccttttttCCGCGCCTGCCTGTTTGGAGGGGGTTGGCCGGCTACGGCTCCCGATAACCAAACATTTTTGTGGGGTGGGGCATCTTCGGGGGAAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                        <  2:2980893/145‑1 (MQ=255)
                                                                                                       cccTCTTTCCGCGCCTGCCTTGTCGGCGCGGGCTTGCCCGCGAATGCGGCCGATCTGACACCAATTTTGTCTGGGCGGGCCTGTCCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGc                                        <  2:4450460/150‑1 (MQ=255)
                                                                                                       cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATCTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                              <  1:1901925/96‑1 (MQ=255)
                                                                                                       cccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCCATCAGACAACAAGTTTTTCTGGTCGGGCCACTTAGCGGGCAAGCCCCCACCAACAAGGGTTTGTGGCATGGCAATACTGATTTTCCAATTTTCGTGGc                                        >  2:6155046/1‑150 (MQ=255)
                                                                                                        ccTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTGCCCCGCATTCGTCCGATCAGACAACAAATTTTTGAGGTCGGGCCACATAGCGGGCCAGCCCCCTCCAACAAGGCTTTGTTGGAAGGCCTGACTGATTTTACAATTTTAGTgggg                                       >  2:5521312/1‑148 (MQ=255)
                                                                                                         cTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCCCGGATGCGTGCCGTAACACAACAAGTTTGTCTGGGGCGGCCTGATCGCGGGCAAGCCCGCGCCAACAAGGGTTTGTGGCATGGCAAGACTGATTTTCCAATTTTAGGGGCGt                                      >  2:2663887/1‑150 (MQ=255)
                                                                                                          tATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTCCCCCCAATTCGTCCGATAACACAACAATGTTTTCTGGGCGGGCCCCATCGCCGGCAAGCCCCCGCCGACAAGGGGTTGTGGCATGGCCAGACTGATTTGCCAAATTTAGTGGCGTa                                     >  2:5125112/1‑150 (MQ=255)
                                                                                                           aTTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                >  1:1037009/1‑74 (MQ=255)
                                                                                                            tttCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAAAATTTGTCTGGTCGGGCCTCTTAGCGGGCAAGCCCGCTCCTACAAGGGTTTTTGGCCTGGCATGACTGATTTGCCATTTTTGGGGGCGTtaa                                   >  2:5720526/1‑150 (MQ=37)
                                                                                                                        ctttGTCGGTGGGGGTTGGCCCGCTAATGCGTCCTAAAACAAAAAATTTGTGTGGGGGCGGGCCTCTTGGCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                       <  1:5025796/148‑1 (MQ=255)
                                                                                                                               ggAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc                                                                                                                  <  1:1673956/52‑1 (MQ=11)
                                                                                                                               ggAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAATTGTTGTCTGATc                                                                                                                  >  2:1673956/1‑52 (MQ=11)
                                                                                                                                        ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa       >  1:2485132/1‑150 (MQ=255)
                                                                                                                                        ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa       >  1:5089658/1‑150 (MQ=17)
                                                                                                                                        ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa       >  1:2404310/1‑150 (MQ=17)
                                                                                                                                        ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa       >  1:2404366/1‑150 (MQ=17)
                                                                                                                                          gCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGc     >  2:4254679/1‑150 (MQ=255)
                                                                                                                                           cccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGGCGGGCCTCTTCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt    <  2:103982/150‑1 (MQ=255)
                                                                                                                                           cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt    >  2:3285808/1‑150 (MQ=18)
                                                                                                                                           cccGCCATTGCGCCCTAACACACAACAATTTGGTCTGGGCGGCCCTCTGCGCGGCCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt    <  2:104155/150‑1 (MQ=255)
                                                                                                                                            ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg   >  1:6219391/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg   >  1:4788429/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATGCGTCCGATCAGACAACAATTATGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg   >  2:458198/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)
                                                                                                                                                    |                                                                                                                                              
CGGTGGGCAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGA  >  NC_002947/444942‑445232

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: