Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,164,567 A→T 5.9% intergenic (+120/+81) PP_4544 → / ← PP_4545 hypothetical protein/beta‑ketodecanoyl‑ACP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,164,5670AT 5.9% 490.8 / 14.3 152intergenic (+120/+81)PP_4544/PP_4545hypothetical protein/beta‑ketodecanoyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base A (49/94);  new base T (2/7);  total (51/101)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 153 total aligned reads displayed.

TGCAGGTAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGA  >  NC_002947/5164419‑5164716
                                                                                                                                                    |                                                                                                                                                     
tGCAGGTAGCCCTGGAAGCCAAGTTCTACGCAACACCCTGCTCTGCTTCACCTGAGGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTAt                                                                                                                                                      <  1:998590/150‑1 (MQ=255)
     gTAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATcacca                                                                                                                                                 <  2:4960519/150‑1 (MQ=255)
      tAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATcaccac                                                                                                                                                <  2:2146253/150‑1 (MQ=255)
         aCCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATACGGCGTTTGCCGCTTATTTTTTTGTGTGTGTAATGGGCCGGCCCCGTAGCGGGGTTGGCGGCTACCACAGGATGAACACAGa                                                                                                                                             >  2:6554965/1‑150 (MQ=255)
             ggAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCtt                                                                                                                                         <  1:663332/150‑1 (MQ=255)
             ggAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCGGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCtt                                                                                                                                         <  1:7331893/150‑1 (MQ=255)
               aaGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGACGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTga                                                                                                                                       <  1:1223283/150‑1 (MQ=255)
               aaGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTga                                                                                                                                       <  2:796542/150‑1 (MQ=255)
                ctggccaaaggctcgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTgag                                                                                                                                      <  2:971623‑M1/45‑1 (MQ=255)
                aGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTgag                                                                                                                                      <  1:2047498/150‑1 (MQ=255)
                  ccAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGaa                                                                                                                                    <  2:1750420/150‑1 (MQ=255)
                      gTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACtgt                                                                                                                                <  2:1536636/150‑1 (MQ=255)
                      gTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACtgt                                                                                                                                <  2:1536490/150‑1 (MQ=255)
                       ttCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAACCGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACtgtg                                                                                                                               <  1:1170064/150‑1 (MQ=255)
                       ttCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGACCtgtg                                                                                                                               <  1:2153690/150‑1 (MQ=255)
                       ttCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCGTGCCCCCGCCCACACTATAAACACAGACCAGGAGAAcaggg                                                                                                                               >  2:8316697/1‑146 (MQ=255)
                         cTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCt                                                                                                                             <  2:2429174/150‑1 (MQ=255)
                            tcgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGGTACACCACAGACCTTGAGAACTGTGCTGTa                                                                                                                          >  2:8640213‑M1/94‑150 (MQ=255)
                            tcgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGGTACACCACAGACCTTGAGAACTGTGCTGTa                                                                                                                          >  1:7693344‑M1/94‑150 (MQ=255)
                             gCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTAc                                                                                                                         <  2:4305174/150‑1 (MQ=255)
                               ttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCt                                                                                                                       >  2:363400‑M1/91‑150 (MQ=255)
                               ttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCt                                                                                                                       >  1:6332843‑M1/91‑150 (MQ=255)
                                acacCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCtg                                                                                                                      <  1:8567338/150‑1 (MQ=255)
                                 cacCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCtgt                                                                                                                     <  2:3718960/150‑1 (MQ=255)
                                  ggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggctcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTTCTGTACCtgtg                                                                                                                    >  2:3619984‑M1/88‑150 (MQ=255)
                                  ggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGTGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCtgtg                                                                                                                    >  1:6531976‑M1/88‑150 (MQ=255)
                                  acCCTGCTCTGCTTCACCTGAAGGGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCtgtg                                                                                                                    <  1:5771044/150‑1 (MQ=255)
                                   gcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                   >  1:7563391‑M1/87‑150 (MQ=255)
                                   gcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                   >  2:7565528‑M1/87‑150 (MQ=255)
                                   cccTGCTCTGCTTCACCTGACGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                   <  2:224662/150‑1 (MQ=255)
                                   cccTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                   <  2:6135175/150‑1 (MQ=255)
                                   cccTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                   <  2:8904570/150‑1 (MQ=255)
                                   cccTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                   <  2:2195328/150‑1 (MQ=255)
                                   cccTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                   <  2:6818839/150‑1 (MQ=255)
                                    ccTGCTCTGCTTCACCTGAAGCGGGAAACGCCGCCCCCATACGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCCCAGTATCACCACACACCATGGGTACAGTGCTGTACCTGTggg                                                                                                                  >  2:9089672/1‑150 (MQ=255)
                                     gtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGa                                                                                                                 >  1:971633‑M1/85‑150 (MQ=255)
                                     gtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcccCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGa                                                                                                                 >  1:4137472‑M1/85‑150 (MQ=255)
                                     ctgctctgctTCACCTGAAGCGGGAAACGCCGACCCCACCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGa                                                                                                                 <  2:5749705/150‑1 (MQ=255)
                                      tgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcccCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAg                                                                                                                >  2:7305179‑M1/84‑150 (MQ=255)
                                      tgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcccCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAg                                                                                                                >  2:7305266‑M1/84‑150 (MQ=255)
                                       gttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGc                                                                                                               >  1:8454710‑M1/83‑150 (MQ=255)
                                       gttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcccCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGc                                                                                                               >  2:6724398‑M1/83‑150 (MQ=255)
                                       gttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcccCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGc                                                                                                               >  2:4812092‑M1/83‑150 (MQ=255)
                                         tgagctcgtttggtgcgaggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCgg                                                                                                             >  1:7646219‑M1/81‑150 (MQ=255)
                                         tctgctTCACCTGACGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAGCgg                                                                                                             <  1:5614594/150‑1 (MQ=17)
                                         tctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCgg                                                                                                             <  2:7775226/150‑1 (MQ=255)
                                         tctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAGCgg                                                                                                             <  1:6648352/150‑1 (MQ=18)
                                             ctTCACCTGCAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCATTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCaa                                                                                                         <  2:6413428/150‑1 (MQ=255)
                                             ctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCaa                                                                                                         <  2:1068790/150‑1 (MQ=255)
                                                 tttggtggggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACATTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGccc                                                                                                     >  1:3316379‑M1/73‑150 (MQ=255)
                                                 gtttgggggggggttttcgtttgtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGccc                                                                                                     >  1:6048289‑M1/73‑150 (MQ=255)
                                                   tggtggggggtattcgattggtagtttgattctgcgcatgcgctgattcatctggcgccttttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACGTGTGGTAGCGGGCAAGCCcgc                                                                                                   >  2:837104‑M1/71‑150 (MQ=255)
                                                    tGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCCCAGTATCACCACAGACCTGGGGGACTGTGCTGTACCTGGGGGAGCGGGCCAGCCcgcg                                                                                                  >  2:8266028/1‑150 (MQ=39)
                                                       cgCGGTAAACGCCGACCCCAGCGGGGGTTCCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACATACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAag                                                                                               <  1:3721036/149‑1 (MQ=37)
                                                                    gACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCCCTCCCACAGTATAACCACAGACCATTATAACATTGCTGTACCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAgg                                                                                  >  2:8303046/1‑150 (MQ=38)
                                                                        ccATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCCCTCCCACAGTATCACCACAGACCTGGAGAACTGTGCTGTACCGGGGGGAGCGGGCAAGCCCGCGAAGAGGCCGGAACAGGCGcc                                                                              >  2:8150283/1‑150 (MQ=38)
                                                                           gcggggttttCCCGCTTATTTTTTTCTGTGTTACCTGCGCCGGCCCCTGCGGGGGCTGGCCCGCACCCACATTAACACCACCTACATTGATAAGTGTGCTGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCtga                                                                           <  1:8187983/143‑1 (MQ=19)
                                                                                 gTTTCCCGCTTCTTTTTGTCTGTGTGAATTGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACTGACCTTTATACCTGTGGTGTACCTGTGGGCGGGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATc                                                                     <  2:3183717/150‑1 (MQ=38)
                                                                                  gtccccGCTTCTTTTGGTCGGTGTGAATTGCGCCGGCCCCTTCGCGGGCTGGCCCGCTCCCACATTAACCCCACTGACCTTGAGACCTGGGCGGTACGTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCa                                                                    <  1:8316697/147‑1 (MQ=25)
                                                                                       cGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACACTAACACCCCAGCACATGAGAACTGTGCTGTAACAGTGGGTGCGGGCAAGCCCCCGAAGAGGCCGTAAACCGCGCCAGATGAAACAGCGCt                                                               >  2:8791873/1‑150 (MQ=25)
                                                                                        gCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                   <  1:3743980/97‑1 (MQ=255)
                                                                                        gCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                   >  2:3743980/1‑97 (MQ=255)
                                                                                             tttttGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCCCCCCCAATATAACCACAGACCTTGAGAACAGTGCAGTACCAGTGGGAGCGGGCAAGCCCCCGAAGAGGCCGGAACAGGCGCCAGATGAAAAAGCCCGTGCGCa                                                         >  2:7474496/1‑150 (MQ=255)
                                                                                                  gTCTGTGTGAAATGCGCCGGCCCCTGCGCGGGCTGGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATc                                                    <  2:8940912/150‑1 (MQ=255)
                                                                                                    cTGTGTGAAAGGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATcac                                                  <  2:8336240/150‑1 (MQ=255)
                                                                                                        gtgAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCAACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTa                                              <  2:4377146/150‑1 (MQ=255)
                                                                                                         tgaacggcgcCGGCCCCTTCGCGGGCTGGCCCGCCCCCCCCGTATCACCCCAGACCTTGAGAACTGTGCTGTACCTGGGGGAGCGGGCACGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTAc                                             <  2:6424531/144‑1 (MQ=255)
                                                                                                         tgAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCACCCCCAGTACCACCCCAGACCTTGAGAACTGTGGTGTACCTGTGGGAGGGGGAACGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTAc                                             <  2:6424530/150‑1 (MQ=255)
                                                                                                            aaTGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCCCAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCaa                                          <  2:7265473/150‑1 (MQ=255)
                                                                                                                 gcCGGCCCCTTCGGGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGaa                                     <  2:4308858/150‑1 (MQ=255)
                                                                                                                   cGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATa                                   <  2:5237298/150‑1 (MQ=255)
                                                                                                                      ccccTTCGCGGGCTTGCCCGCACCCCCAGTAGAACCACAGACCTTGAGACCGGTGCTGTACCTGGGGGCGCGGGCACGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATAccc                                <  1:9026507/150‑1 (MQ=255)
                                                                                                                               gggttggcccgcccCCACATTATCACCACACACTTTGATATGGGGGGTGTACCGGGGGGTGGGGGCACGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCaaa                       <  1:7754675/144‑1 (MQ=255)
                                                                                                                               GGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATAACACTACCAATCGAATACCCCGCACCaaa                       <  1:7281251‑M1/150‑1 (MQ=255)
                                                                                                                                GGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAAc                      >  2:9141591‑M1/1‑150 (MQ=255)
                                                                                                                                 GCTTGCCCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACg                     <  2:3712276‑M1/150‑1 (MQ=255)
                                                                                                                                     GCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCt                 >  2:5508882‑M1/1‑150 (MQ=255)
                                                                                                                                     GCCCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCt                 <  2:3880071‑M1/150‑1 (MQ=255)
                                                                                                                                      CCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTc                >  2:8491715‑M1/1‑150 (MQ=255)
                                                                                                                                      CCCGCTCCAACAGTATCACCACAGACCTTGTGAACTGTGCTGTACCTGTGGGAGCGTGCATGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTc                >  2:6432295‑M1/1‑150 (MQ=255)
                                                                                                                                       CCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa               >  2:9051252‑M1/1‑150 (MQ=255)
                                                                                                                                       CCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa               <  1:3271522‑M1/150‑1 (MQ=255)
                                                                                                                                       CCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa               <  1:7801879‑M1/150‑1 (MQ=255)
                                                                                                                                       CCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa               <  1:7139239‑M1/150‑1 (MQ=255)
                                                                                                                                        CGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCaa              <  2:1189434‑M1/150‑1 (MQ=255)
                                                                                                                                        CGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACACCGCACCAAACGAGCTCaa              <  2:7639153‑M1/150‑1 (MQ=255)
                                                                                                                                         GCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCaaa                       >  1:1958459‑M1/1‑140 (MQ=255)
                                                                                                                                         GCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCaaa                       <  2:1958459‑M1/140‑1 (MQ=255)
                                                                                                                                           TCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAacac           <  2:1669586‑M1/150‑1 (MQ=255)
                                                                                                                                           TCCCACAATATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAacac           <  1:6546764‑M1/150‑1 (MQ=255)
                                                                                                                                             CCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc         <  1:6347028‑M1/150‑1 (MQ=255)
                                                                                                                                             CCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc         <  1:5355818‑M1/150‑1 (MQ=255)
                                                                                                                                              CACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGcc                                                                                          >  1:964332‑M1/1‑68 (MQ=255)
                                                                                                                                              CACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGcc        <  1:2962561‑M1/150‑1 (MQ=255)
                                                                                                                                              CACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGcc        <  2:262630‑M1/150‑1 (MQ=255)
                                                                                                                                              CACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGcc        <  2:4512013‑M1/150‑1 (MQ=255)
                                                                                                                                               ACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGccc       <  1:8242958‑M1/150‑1 (MQ=255)
                                                                                                                                               ACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGccc       <  2:7054313‑M1/150‑1 (MQ=255)
                                                                                                                                                 AGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGAGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCaa     <  1:5110598‑M1/150‑1 (MQ=255)
                                                                                                                                                    aTCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGa  <  1:2614403/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCAGGTAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGA  >  NC_002947/5164419‑5164716

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: