Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,098 T→C 8.7% intergenic (+68/‑220) 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,0980TC 8.7% 269.6 / 9.6 95intergenic (+68/‑220)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (33/51);  new base C (4/4);  total (38/57)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTT  >  NC_002947/444949‑445246
                                                                                                                                                      |                                                                                                                                                     
cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCATGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                                       <  2:2461800/150‑1 (MQ=255)
cAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGAt                                                                                                                                                       <  2:3618671/150‑1 (MQ=255)
 aaGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATc                                                                                                                                                      <  2:598326/150‑1 (MQ=255)
  aGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa                                                                                                                                                     <  1:4304611/150‑1 (MQ=255)
  aGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa                                                                                                                                                     <  1:4304607/150‑1 (MQ=255)
  aGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCa                                                                                                                                                     <  1:2929488/150‑1 (MQ=255)
    aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGGCCGCGAATGCGTCCGATCAGa                                                                                                                                                   <  2:3810431/150‑1 (MQ=255)
    aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                                   <  1:928917/150‑1 (MQ=255)
    aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                                   <  1:928267/150‑1 (MQ=255)
    aCCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGa                                                                                                                                                   <  1:2550864/150‑1 (MQ=255)
        aaGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGacaac                                                                                                                                               <  1:1211236/150‑1 (MQ=255)
           tCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAt                                                                                                                                           <  1:952007/150‑1 (MQ=255)
           tCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAt                                                                                                                                           <  1:951965/150‑1 (MQ=255)
            cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                          <  2:931152/150‑1 (MQ=255)
            cGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtt                                                                                                                                          <  2:152969/150‑1 (MQ=255)
              aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                        <  2:2591090/150‑1 (MQ=255)
              aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                        <  1:779342/150‑1 (MQ=255)
              aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                        <  1:4272677/150‑1 (MQ=255)
              aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                        <  1:4272587/150‑1 (MQ=255)
              aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                        <  1:2871040/150‑1 (MQ=255)
              aGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAAtttt                                                                                                                                        <  2:3949727/150‑1 (MQ=255)
                cTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTTAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGt                                                                                                                                      <  1:4239790/150‑1 (MQ=255)
                 tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCGCGCATGCGGGCCCTCAGCCAACACAtttttg                                                                                                                                      >  2:3966465/1‑148 (MQ=255)
                 tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                     <  1:3665901/150‑1 (MQ=255)
                 tCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTc                                                                                                                                     <  2:858133/150‑1 (MQ=255)
                    gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCCCGCAAGCGTCCCATCACACAACCAATTTGtgttg                                                                                                                                  >  2:3149771/1‑146 (MQ=255)
                    gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTgg                                                                                                                                  <  2:1216446/150‑1 (MQ=255)
                    gcgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTgg                                                                                                                                  <  2:1216374/150‑1 (MQ=255)
                     cgcAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGt                                                                                                                                 <  2:3106773/150‑1 (MQ=255)
                       cAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAATCGTCCGATCAGACAACAAATTTTTGTGGtgg                                                                                                                               >  2:1848665/1‑148 (MQ=255)
                        aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                              <  2:3269220/150‑1 (MQ=255)
                        aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCgg                                                                                                                              <  1:1119070/150‑1 (MQ=255)
                        aaaGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGACTTTCCCGCGAATGCGTCCGATAAGACAACAATTTTGTCTGGTCgg                                                                                                                              >  2:1334246/1‑150 (MQ=255)
                          aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTTCCCGCGAATGCGTCCGATCAGAAAAAAATTTTGTCTGGGAGGGc                                                                                                                            >  2:2167607/1‑150 (MQ=255)
                          aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                            <  1:1285967/150‑1 (MQ=255)
                          aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                            <  1:1286212/150‑1 (MQ=255)
                          aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGc                                                                                                                            <  1:1286283/150‑1 (MQ=255)
                          aGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCCAAAGCGTCCGATCAGAAAAAAAATTTGTGTGGTGGGGc                                                                                                                            >  2:2167661/1‑150 (MQ=255)
                              cgtcAAGGCATCCGACGAAATGCTTGGCCTGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCggacgcat                                                                                                                        <  1:1362749/150‑7 (MQ=21)
                                  aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgc                                                                                                                    <  1:2671737/150‑1 (MQ=255)
                                  aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgc                                                                                                                    <  1:415889/150‑1 (MQ=255)
                                   aGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTcgct                                                                                                                   <  1:4091379/150‑2 (MQ=255)
                                   aGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATGTTGTCTGGTCGGGCCTCTTcgcg                                                                                                                   <  2:3833341/150‑1 (MQ=255)
                                    ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACGTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTGTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                  <  2:4068803/150‑1 (MQ=255)
                                     gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCggg                                                                                                                 <  2:406792/150‑1 (MQ=18)
                                     gCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACGTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCggg                                                                                                                 <  2:1119492/150‑1 (MQ=255)
                                      cATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGc                                                                                                                <  1:235520/150‑1 (MQ=255)
                                        tCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTATGGTCGGGCCACTTCGCCGGCaa                                                                                                              >  2:2967771/1‑150 (MQ=255)
                                          cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGc                                                                                                            <  1:1532747/150‑1 (MQ=255)
                                          cgacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAAGCGTCCGATCAGACAACCATTTTGTATGGTGGGGCCTCGTCGCGGGCAAGc                                                                                                            >  2:2628377/1‑150 (MQ=255)
                                           gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGcc                                                                                                           <  1:1607538/150‑1 (MQ=255)
                                           gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGcc                                                                                                           <  1:1438010/150‑1 (MQ=255)
                                           gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGcc                                                                                                           <  2:837072/150‑1 (MQ=255)
                                           gacgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTGCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCGCGAATGCGTCCGATCAGACAACACTTGTGTCTGGTCGGGCCTCTTCGCGGGCAAGcc                                                                                                           <  1:2158132/150‑1 (MQ=255)
                                                      gggCCGCATCATCGACACCTACGCCTTAATACTTGCGACATCCCCTTTTTCAGCGCCTGCGTTGGCGGAGGGGGCTTGCCGGCTAATGCGTCCTATCTTAAAACATTTTTGTGTGGGGGGGCCTCTCCGGGGGCAAGCCCGCTCCTACaa                                                                                                <  2:1463008/149‑1 (MQ=255)
                                                         ccGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGATGCCCCCCAATTCCTACGATAAGACAACAATTTTGTCAGGTCGGGCCTCTTCCAGGGCCATCCCCCGCCGACAAcgg                                                                                             >  2:1019611/1‑150 (MQ=255)
                                                                 tcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCGGAAGCCGCCGGGCAGGCAACCATTTTGTGGGGGCGGGGCGCATCGCGGGCCAGCCCCCGCCGACAACGGGTTGtggg                                                                                     >  2:2531054/1‑149 (MQ=255)
                                                                 tcGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCCCATGCGGCCCCTCAGACAACCATTTTGTCGGGGGGGGCCGCATGGCGGGCCACCCCGCGCCGACAAGGGGTTGtggg                                                                                     >  2:2531046/1‑149 (MQ=255)
                                                                    acacCTCCGCCTCAATACTGGCGACCCCCCCTTTTGCCGCGCCTCCCTTGTGGGCGGGGGCTGGCCCGGTATGGGGCCCTATCACACACCATTTTTGTGGGGGGGGGCCGCTCCGGGGGCACGCCCGCTCCTACAAGGGTTTGtggcatg                                                                                  <  2:3895915/150‑1 (MQ=255)
                                                                      cccTACGCCTGACTCCGGGCGTCACCCCCTATTTCCGCGCCGGCCTTGTAGGAGCGGGCTTGCCCGCGACTGGGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatggc                                                                                <  2:2568861/149‑1 (MQ=25)
                                                                      acCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCTAATGCGCCCGATCAGACAACATTTTTGTCTGGTCGGGCCTCTTGGCGGGCAAGCCCGCTCCTACAAGGGTTTGtggcatggc                                                                                <  2:2568753/150‑1 (MQ=255)
                                                                       cctcccccTCAATCCTGGCGACACCCCCTATTCCAGCCCCGGCCTTGTAGGAGGGGGCTGGCCCGCCAAGGCGCCCGAACAGACACCAATTTTGGGGGGCGGGGCCTCCCCGGGGCAAACCCCGCTCCTACAAGGGTTTGtggcatggca                                                                               <  2:3440447/144‑1 (MQ=255)
                                                                              cTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTTCCCCCCCAAGGCGCCGCTCAGAAAAAAAATTTTTATGGTGGGGGCTCGTAGCCGGCCACCCCCCGACCAAAACGGGGTGTGGGATGGGATGaaata                                                                        >  2:1747790/1‑146 (MQ=255)
                                                                                           tCACCCCCTATTTCCGCGCCGGCCTGGTCGGCGGGGGCTGGCCCGCGAATGCGTCCGATCATCAAAAAATTTGGCCGGGTCGGGCCTCTTCGGGGGAACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTc                                                           <  2:3273040/150‑1 (MQ=255)
                                                                                            cATCCCCTATTTCAGCGCCTGCCTTGTCGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGCCAACAATTTGGTCTGGGCGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCg                                                          <  2:4107782/150‑1 (MQ=39)
                                                                                             ccccccTCTTCCCGCCCCGGCGTTGTCGGCGCGGGCTTGCCCGCGAATGGGCCCGAACACACAACAATTTTGGGGGGGCGGGCCCCTGCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGt                                                         <  2:3974916/148‑1 (MQ=255)
                                                                                             aTCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGACCATACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGt                                                         <  2:4013616/150‑1 (MQ=255)
                                                                                              tCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGTTTGCCGCCGCAGCGTACCCTACTACAACAAATTTTTGGGGTGGGGCGGCTTAGCGGGCAACCCCGCGCCTACCAGGGTTTGTTGGATGGGCTGACAGGTTTTACAattttattt                                                        >  2:4401789/1‑146 (MQ=255)
                                                                                                  cTATTCCCGCGCCTGCCTTGTCGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGt                                                    <  2:710547/150‑1 (MQ=255)
                                                                                                            ccctgcctTGTTGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTGTGTCTGGTCGGCCCTATTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCAcc                                          <  2:4428267/149‑1 (MQ=38)
                                                                                                                       attcgCGGGCAAGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATcc                               >  2:4095695/5‑150 (MQ=255)
                                                                                                                                  ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                     >  2:3317202/1‑150 (MQ=255)
                                                                                                                                  ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                     >  1:4325168/1‑150 (MQ=255)
                                                                                                                                  ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                     >  1:893715/1‑150 (MQ=17)
                                                                                                                                   tGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGAGTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAg                    >  1:3261976/1‑150 (MQ=255)
                                                                                                                                      ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGGCTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                 >  1:1762061/1‑150 (MQ=255)
                                                                                                                                      ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                 >  1:4502842/1‑150 (MQ=255)
                                                                                                                                      ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                 >  2:1427949/1‑150 (MQ=255)
                                                                                                                                      ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                 >  1:1996289/1‑150 (MQ=255)
                                                                                                                                      acgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                 >  1:3543912/2‑150 (MQ=255)
                                                                                                                                       cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCg                                                          >  2:3888550/1‑108 (MQ=255)
                                                                                                                                       cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCg                                                          <  1:3888550/108‑1 (MQ=255)
                                                                                                                                        gcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAg               >  2:3781449/1‑150 (MQ=255)
                                                                                                                                         ccacggcgcCCGACCACCCAACATGTTGGTCGGGCCGGCCCCCGCCGCGGCCACCCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGc              <  1:2842418/145‑1 (MQ=255)
                                                                                                                                            aTGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaa           >  2:3854862/1‑150 (MQ=255)
                                                                                                                                                gcccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt       >  1:1581323/3‑150 (MQ=38)
                                                                                                                                                gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt       >  1:1747552/1‑150 (MQ=255)
                                                                                                                                                 tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa      >  2:3311984/1‑150 (MQ=255)
                                                                                                                                                 tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa      >  1:2979977/1‑150 (MQ=255)
                                                                                                                                                  ccGACCAGACAAAATTGTTGTCTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa     >  1:2766510/1‑150 (MQ=38)
                                                                                                                                                  ccGACCAGACAAAATTGTTGACTGATCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTaa     >  2:3440705/1‑150 (MQ=38)
                                                                                                                                                    gATCAGACAACAATGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCAAAGTCCATCCGTGACAGAAAGAAGAGATAAACGTAAGt   >  1:1694232/1‑150 (MQ=255)
                                                                                                                                                     aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt  >  2:1197813/1‑150 (MQ=255)
                                                                                                                                                     aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt  >  2:2139359/1‑150 (MQ=255)
                                                                                                                                                     aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt  >  2:2230677/1‑150 (MQ=255)
                                                                                                                                                     aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt  >  2:2139525/1‑150 (MQ=255)
                                                                                                                                                     aTCAGACAACAATTTTGTCTGCACGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTTTGGCATGGCATGACAGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAgat  >  1:2794525/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CAAGACCGAAGTCGAGCTCGGCGCAAAGGTCGTCAAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTT  >  NC_002947/444949‑445246

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: