Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 445,123 G→A 17.5% intergenic (+93/‑195) 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,1230GA17.5% 219.1 / 34.8 95intergenic (+93/‑195)PP_0367/PP_0368hypothetical protein/acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (50/27);  new base A (10/6);  total (60/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

AAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGC  >  NC_002947/444983‑445267
                                                                                                                                             |                                                                                                                                                
aaGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGGGTGCCCGCGCATGCGGCCCAGCCAACAAAAAAATTTTATGGTGGGGccacgtagc                                                                                                                                         >  1:3372144/1‑143 (MQ=255)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                       <  2:2574791/150‑1 (MQ=255)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                       <  1:172720/150‑1 (MQ=255)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                       <  1:666944/150‑1 (MQ=255)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                       <  1:4537837/150‑1 (MQ=255)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                       <  1:4256776/150‑1 (MQ=255)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCgg                                                                                                                                       <  2:4600174/150‑1 (MQ=255)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                       <  1:1207970/150‑1 (MQ=18)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                       <  2:1803560/150‑1 (MQ=18)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                       <  2:1150259/150‑1 (MQ=18)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                       <  1:930759/150‑1 (MQ=18)
  ggCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTGGCGTCATCCCCTATTTCCGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCgg                                                                                                                                       <  2:1150278/150‑1 (MQ=17)
     aTCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCa                                                                                                                                    <  2:3042089/150‑1 (MQ=255)
      tCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCaa                                                                                                                                   <  1:3413524/150‑1 (MQ=255)
      tCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCaa                                                                                                                                   <  2:3601502/150‑1 (MQ=255)
      tCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCaa                                                                                                                                   <  1:2521992/150‑1 (MQ=17)
          acgaAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAAAAAATTTGTCTGGTCGGGCCGCTTAGCGGGCAAGccc                                                                                                                               >  2:1679111/1‑150 (MQ=255)
              aaTGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGGATGCGTCCGATCAGAAAAACAATTTTTCTGGTCGGGCCTCGTAGCGGGCAAGCCCGCTc                                                                                                                           >  2:3392637/1‑150 (MQ=255)
                      ggcGGATCTTCGCCACCTCCGCCTGAATAGTTGCGTCATCCCCTATTCCAGCGCCTGCCTTGTAGGAGGGGGCTTGCCCGCGAATGCGTCCGATCAGACAACATTTTTGTCTGGTGGGGCATCTTCGGGGGCACGCCCGCTCCTACAAgg                                                                                                                   <  2:4661650/148‑1 (MQ=255)
                                      cccccccTTAATCCGTGCGTCAGCCCCTCTTTCCGGGCCGGCCTTGTGGGCGGGGGCTGGCCCGCGTCTGGGACCTAACACAAAACTTTTTTGTGTGGGGGGGCCCCGGCGGGGGAACCCCCGCTCCTACAAGGGTTTGTggcatggcat                                                                                                   <  2:689101/145‑1 (MQ=255)
                                                      gCGTCATCCCCTATTTCCGCGCCTGCGTTGTAGGCGCGGGCTTGCCCGCGAATGCGTCCGACCAGACAACAATTTTGACGGGTGGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCAttt                                                                                   <  2:4211511/150‑1 (MQ=38)
                                                         tCATCCCCTTTTACCGCGCCGGCGTTGTTGGCGGGGGCTTCCCGGCTATTGGGTCCGATCAGACAACAATTTGTTGGGGGGGGGCCTCTTCGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTc                                                                                <  2:2693252/150‑1 (MQ=255)
                                                                cTATTTCAGCGCCTGCCTTGTAGGAGCGGGCGTGCCCCCCAAAGCGTCCGCTCACACAACAAATTTTTCTGGTCGGGCCGCGTCGCGGGCAAGCCCCCGCCTACAACGGTTTTTTGCCTGGCATGACTGATTTTGCAATTTTAGTGGCgg                                                                         >  2:2030608/1‑149 (MQ=255)
                                                                             tgcctTGTCGGCGGGGGCTTGCCCGCTACTGCGTCCGACCAGACAAAAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATa                                                            <  2:2041390/150‑1 (MQ=38)
                                                                               cgttgttgGCGGGGGCTGGCCGGCGATTGCGCCCGATCAGCAAACTTTTTTGTGTGTTGGGGCCTCTCCGGGGGAAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGt                                                          <  2:4401228/148‑1 (MQ=255)
                                                                                    tcggtgggggCTGGCCCGCTACGGCGGCCTATCACCCAACTTTTTTGTCGGGTGGGGCCTCGGCGCGGGCCCGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATc                                                     <  1:3389136/143‑1 (MQ=255)
                                                                                       gAGGGGGCTTGCCGGCTACTGCTTCCTATATGACAACTTTTTTGTGTGGTGGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATACtt                                                  <  2:2402119/150‑1 (MQ=255)
                                                                                         gCGGGCTTGCCCGCGACTGCGTCCGATCAGACAACATTTTTGTCTGGGCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                                <  2:4588703/150‑1 (MQ=255)
                                                                                         gCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTGTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTTTGGCATGGCATGTCTGATTTGGCATTTTTGGTGGCGTAATTGGCACAATAGTCAATCCTTgt                                                >  2:3037697/1‑149 (MQ=38)
                                                                                         gCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGa                                                >  2:1103995/1‑150 (MQ=255)
                                                                                               ttGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                          >  2:4768524/1‑150 (MQ=255)
                                                                                               ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                          >  2:1644434/1‑150 (MQ=17)
                                                                                               ttGCCCGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGaaa                                          >  1:4357375/1‑150 (MQ=17)
                                                                                                  cccGCGAAGAGGCCCGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCt                                       >  2:4402487/1‑150 (MQ=18)
                                                                                                   ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  2:3221015/1‑150 (MQ=255)
                                                                                                   ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  2:1306432/1‑150 (MQ=255)
                                                                                                   ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTg                                      >  2:1550485/1‑150 (MQ=255)
                                                                                                   ccGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCAAGGCATGACTGATATGCAATTTTTCGTGGCGTAAATGGAACCATAGTCCATCCTTGACAGAAAGCTg                                      >  1:4309837/1‑150 (MQ=255)
                                                                                                    cgctAATGCGTCCTATCAGAAAACAATTTTGTCTGGTCGGCCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                     <  2:4826484/150‑1 (MQ=255)
                                                                                                    cgcgAATGCGTCCGATCAGATAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                     >  1:568558/1‑150 (MQ=255)
                                                                                                    cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCTTAGTCCATCCGTGACAGAACGCTGa                                     >  2:4927133/1‑150 (MQ=255)
                                                                                                    cgcgAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGa                                     >  2:4675954/1‑150 (MQ=255)
                                                                                                          gaggcccGACCAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTaaa                               >  1:156420/6‑150 (MQ=21)
                                                                                                            cGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAAct                             >  2:98202/1‑149 (MQ=255)
                                                                                                             gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTg                                                                             <  2:3750749/101‑1 (MQ=255)
                                                                                                             gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTg                                                                             >  1:3750749/1‑101 (MQ=255)
                                                                                                             gTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGt                            >  1:1694559/1‑150 (MQ=255)
                                                                                                              tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                           >  2:4286242/1‑150 (MQ=255)
                                                                                                              tCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTa                           >  2:4286272/1‑150 (MQ=255)
                                                                                                                 gaccAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt                        >  2:1346875/4‑150 (MQ=25)
                                                                                                                 gaccAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt                        >  1:2708039/4‑150 (MQ=25)
                                                                                                                 gATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt                        >  2:1614955/1‑150 (MQ=255)
                                                                                                                 gATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt                        >  1:2030638/1‑150 (MQ=255)
                                                                                                                 gATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt                        >  2:1015941/1‑150 (MQ=255)
                                                                                                                 gATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGt                        >  2:44338/1‑150 (MQ=255)
                                                                                                                  accAGACAAAATTGTTGTCTGATCGGACCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                       >  2:2911995/3‑150 (MQ=38)
                                                                                                                  accAGACAAAATTGTTGTCTGATCGGACCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                       >  2:354020/3‑150 (MQ=38)
                                                                                                                  aTCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGtt                       >  2:4117861/1‑150 (MQ=255)
                                                                                                                   tCAGACAAACATTTTGTGTGGTCGGGCCTCTTCGGGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTACACGTAAGttt                      <  2:4054290/150‑1 (MQ=255)
                                                                                                                    cAGACAAAATTGTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc                     >  1:4249955/1‑150 (MQ=255)
                                                                                                                    cAGACAAAATTGTTGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc                     >  1:515933/1‑150 (MQ=37)
                                                                                                                     aGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCa                    >  2:3628506/1‑150 (MQ=255)
                                                                                                                      gACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaa                   >  2:4258751/1‑150 (MQ=255)
                                                                                                                        acaaaaTTTTTGTGTGGTCGGGCCTCTGCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCaaa                  <  2:559440/150‑1 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACa                                              >  1:4121264/1‑120 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACa                                              <  2:4121264/120‑1 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc                     <  2:1415356/145‑1 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTc                     >  1:1415356/1‑145 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  1:2418884/1‑150 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  1:4861589/1‑150 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  1:3225240/1‑150 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  2:556086/1‑150 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  1:4535454/1‑150 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  1:4649429/1‑150 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  2:138848/1‑150 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  2:3049029/1‑150 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaaca                >  2:962035/1‑150 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCCCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGc                                                                  >  1:2409011/1‑100 (MQ=255)
                                                                                                                          aacaaTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCCCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGc                                                                  <  2:2409011/100‑1 (MQ=255)
                                                                                                                            catttttGTCTGTTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac              <  2:3141835/147‑1 (MQ=255)
                                                                                                                            caaTTTTGTCTGGTGGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac              <  2:3372144/150‑1 (MQ=255)
                                                                                                                            aattgttGTCTGATCGGACGCATTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAaacaac              >  1:3194112/6‑150 (MQ=38)
                                                                                                                             aaTTTTGTCTGCTCGGGCCTCTTCGCGGGCAACCCCGCTCCTACAAGGGTTTGTGGCATGGCAAGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCAACCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACt             >  2:633162/1‑150 (MQ=255)
                                                                                                                              aTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTg            >  2:4280394/1‑150 (MQ=255)
                                                                                                                              aTTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTg            >  2:3526007/1‑150 (MQ=255)
                                                                                                                               ttttGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGt           >  2:3164485/1‑150 (MQ=255)
                                                                                                                                  tGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTg        >  2:3695840/1‑150 (MQ=255)
                                                                                                                                    tcgggTCGGGCCTCGTGGCGGGCACGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGAt      <  2:2621713/147‑1 (MQ=255)
                                                                                                                                      tGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg    >  2:2515113/1‑150 (MQ=255)
                                                                                                                                      tGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg    >  2:2274166/1‑150 (MQ=255)
                                                                                                                                      tGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg    >  1:771524/1‑150 (MQ=255)
                                                                                                                                      tGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg    >  1:605728/1‑150 (MQ=255)
                                                                                                                                      tGGTCGGCCCTCTACGCGGGCAACCCCGCTCCTACAAGGGATTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCg    <  2:1206050/150‑1 (MQ=255)
                                                                                                                                       ggTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCgg   >  1:90353/1‑150 (MQ=255)
                                                                                                                                        gTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGc  >  2:1998470/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                
AAGGCATCCGACGAAATGCTTGGCCGGATCATCGACACCTACGCCTGAATACTTGCGTCATCCCCTATTTCAGCGCCTGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCGTCCGATCAGACAACAATTTTGTCTGGTCGGGCCTCTTCGCGGGCAAGCCCGCTCCTACAAGGGTTTGTGGCATGGCATGACTGATTTGCCATTTTTCGTGGCGTAAATGGCACCATAGTCCATCCTTGACAGAAAGCTGAGCTAAACGTAAGTTTCAAACAACTGTTTGATCGGC  >  NC_002947/444983‑445267

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: