Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,164,568 T→A 6.1% intergenic (+121/+80) 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,5680TA 6.1% 263.9 / 3.9 86intergenic (+121/+80)PP_4544/PP_4545hypothetical protein/beta‑ketodecanoyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base T (28/52);  new base A (4/2);  total (32/54)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGTAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAG  >  NC_002947/5164420‑5164717
                                                                                                                                                    |                                                                                                                                                     
gCAGGTAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATc                                                                                                                                                      <  1:2027121/150‑1 (MQ=255)
 cAGGTAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATca                                                                                                                                                     <  2:3935287/150‑1 (MQ=255)
 cAGGTAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATca                                                                                                                                                     <  2:1333627/150‑1 (MQ=255)
    gTAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGGCCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCGTGGCCGCTCCCACAGTATcaccc                                                                                                                                                  >  2:3773589/1‑149 (MQ=255)
    gTAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCTTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATcacca                                                                                                                                                  <  1:3006071/150‑1 (MQ=255)
     tAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATcaccac                                                                                                                                                 <  2:3363238/150‑1 (MQ=255)
      acaCCTGCAACACAATTTATAAGCAACACCCTGCACAGCTACACCTCAAGCGGCAAACGCCCACCCCAACCGGCGTTTCCCGCTTCTTATTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACcaca                                                                                                                                                <  1:1088640/148‑1 (MQ=255)
          atgacctggccaaaggctcgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGAcc                                                                                                                                            <  2:1267085‑M1/40‑1 (MQ=255)
            ggAAGCCAATTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCtt                                                                                                                                          <  2:1024479/150‑1 (MQ=255)
            ggAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCtt                                                                                                                                          <  2:1024197/150‑1 (MQ=255)
            ggAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCtt                                                                                                                                          <  2:1429377/150‑1 (MQ=255)
            ggAAGCCAAGTTATAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCtt                                                                                                                                          <  1:4436110/150‑1 (MQ=255)
              aaGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTga                                                                                                                                        <  2:2041190/150‑1 (MQ=255)
                 gggcacaggctcgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGaa                                                                                                                                     >  2:4124849‑M1/104‑150 (MQ=255)
                 ggccaaaggctcgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGaa                                                                                                                                     >  1:1080689‑M1/104‑150 (MQ=255)
                 ggccaaaggctcgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGaa                                                                                                                                     >  1:109‑M1/104‑150 (MQ=255)
                 ggccaaaggctcgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGaa                                                                                                                                     >  2:5039256‑M1/104‑150 (MQ=255)
                 ccAAGTTATAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGaa                                                                                                                                     <  1:3497972/150‑1 (MQ=255)
                   aaGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGTAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACt                                                                                                                                   <  2:1856908/150‑1 (MQ=255)
                   aaGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACt                                                                                                                                   <  2:2022918/150‑1 (MQ=255)
                    caaaggctcgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcccCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACtg                                                                                                                                  <  1:4580641‑M1/50‑1 (MQ=255)
                    aGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACtg                                                                                                                                  <  1:5078272/150‑1 (MQ=255)
                          ctcgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGGTACAGCACAGACCTTGAGAACTGTGCTGt                                                                                                                            >  1:3692048‑M1/95‑150 (MQ=255)
                           aGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTa                                                                                                                           <  2:3765499/150‑1 (MQ=255)
                            gCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTAc                                                                                                                          <  2:4520246/150‑1 (MQ=255)
                            cgttgggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTAc                                                                                                                          >  1:3333070‑M1/93‑150 (MQ=255)
                              aaCACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCt                                                                                                                        <  2:4780449/150‑1 (MQ=255)
                              aaCACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCATAGCGGGCTTGCCCCCTACCACAGTATCACCACACACCTTGAGAACTGTGCTGTaacg                                                                                                                        >  2:5095709/1‑147 (MQ=255)
                                cacCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCtgt                                                                                                                      <  1:4786226/150‑1 (MQ=255)
                                 ggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCtgtg                                                                                                                     >  1:3997997‑M1/88‑150 (MQ=255)
                                 ggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCtgtg                                                                                                                     >  1:1267085‑M1/88‑150 (MQ=255)
                                 ggcgtgttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctggtacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCtgtg                                                                                                                     >  1:4748115‑M1/88‑150 (MQ=255)
                                    ctgctctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGa                                                                                                                  <  2:246348/150‑1 (MQ=255)
                                      gttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGc                                                                                                                >  1:4425879‑M1/83‑150 (MQ=255)
                                      gttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGc                                                                                                                >  1:1847947‑M1/83‑150 (MQ=255)
                                       ttgagctcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcccCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCg                                                                                                               >  2:2986929‑M1/82‑150 (MQ=255)
                                       ctctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAGCg                                                                                                               <  2:2477093/150‑1 (MQ=18)
                                        tctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCgg                                                                                                              <  2:2362557/150‑1 (MQ=255)
                                        tctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCgg                                                                                                              <  1:349911/150‑1 (MQ=255)
                                        cctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCGGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCgg                                                                                                              <  2:3489667/149‑1 (MQ=255)
                                         ctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCggg                                                                                                             <  1:4350969/150‑1 (MQ=255)
                                             tcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                    <  2:2652381‑M1/64‑1 (MQ=255)
                                             tcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcctCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTgg                                                                                                                    >  1:2652381‑M1/76‑139 (MQ=255)
                                             tcgtttggtgcggggtattcgattggtagtgtgattctgcgcaagcgctgattcatcaggcgcctgttacggcccCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAg                                                                                                         >  2:4580641‑M1/76‑150 (MQ=255)
                                                         gggAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCGTGCCCCCTACCACAGTATCACCACCGACCGTTATGACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAgg                                                                                             >  2:3247512/1‑150 (MQ=38)
                                                         gggAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCGTGCCCCCGCCCCCAGGATCCCCCCAGCCCATGAGGACTGGGCGGTAACGGGGGGGGCGGGCCACCCCGCCAAGAgg                                                                                             >  2:3247610/1‑150 (MQ=255)
                                                              aCGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCACAGACCATGAGAACTGTGCGGTACCGGGGGGGGCGGGCAAGCCCCCGAAGAGGGCgga                                                                                        >  2:1642691/1‑148 (MQ=38)
                                                              aCGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCGTGCCCGCGCCCCCAGGATAACCCCAGACCGTGAGAACGGTGCTGTAACGGTGGGCGCGGGCCAGCCCGCGAAGAGGCCgga                                                                                        >  2:1642763/1‑148 (MQ=255)
                                                                      cccATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCACCCACAGTAACACCCCAGACCGTGAGAAATGTTCTGTAACTGTGGGAGCGGGCCAGCCCCCGAACAGGCCGGAACAGgcgc                                                                                >  2:1164186/1‑150 (MQ=255)
                                                                      cccATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCATGCCCGCGCCCACCGGATCACCACACACCATGGGAACTGTGCGGTACCGGTGGGGGCGGGCCAGCCCCCGCAGAGGCCGGGACAGgcgc                                                                                >  2:1164200/1‑150 (MQ=255)
                                                                       ccATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTACGCGGGCTTGCCCGCTCCCACAGTATCACCCCAGACCTTGAGACCTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGcc                                                                               <  2:4219899/150‑1 (MQ=255)
                                                                           cgggCTTTTCCCGCTTCTTTTTGTCTGTGTTAATTGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACTCACATTGATAACTGTGCTGTACCTGTGGGAGCGGGTAAGCCCGCGAAGAGGCCGTAACAGGCGCCtgat                                                                           <  2:1729498/148‑1 (MQ=37)
                                                                           cgggCGTTTCCCGCTTATTTTTGTCTGTGTGAATTGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGTAAGCCCGCGAAGAGGCCGTAACAGGCGCCtgat                                                                           <  2:1729540/148‑1 (MQ=255)
                                                                           ccGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCGGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCtgat                                                                           <  2:1232452/150‑1 (MQ=255)
                                                                              gggTTTCCCGCTTATTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGaa                                                                        <  2:3557822/148‑1 (MQ=255)
                                                                                   tCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCCCAGTATCACCACAGACCTTGAGAACGGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAgc                                                                   <  2:4693575/150‑1 (MQ=255)
                                                                                   ccccGCTTCTTTTTGTCGGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCGCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAgc                                                                   <  2:2942721/149‑1 (MQ=255)
                                                                                     ccGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCCCGCCCACACAAAAAACACACAACATGAGAACACTGCTGTAACTCGTGGAGCCGGCAACCCCCCGCAAACGCCCAAACACGCGCCAGATTAATCaaccc                                                                 >  2:4834924/1‑146 (MQ=255)
                                                                                               ttGTCTGTGTGAAAGGCGCCGGCCCCTGCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGCGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGaa                                                       <  2:3654818/150‑1 (MQ=255)
                                                                                                    tgtgtgAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCt                                                                                                                        <  2:5060934/80‑1 (MQ=255)
                                                                                                    tgtgtgAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCt                                                                                                                        >  1:5060934/1‑80 (MQ=255)
                                                                                                       gtgAAATGCGCCGGCCCCTTCGCGGGCTTTCCCGCTCCCACACTATAAACACACACCTTGAGAAATGTGCTGTAACTGTGGGTGCGGGCAAGCCCCCGAAGAGGCCCGAACAGGCGCCTGATGAATAAGCGCGTGCGCAGTATCaaacaa                                               >  2:1178275/1‑148 (MQ=255)
                                                                                                        tgAAAGGCGCCGCCCCCTGCGCGGGCTTGCCCGCTCCCACCTTCTCACCCCAGACCTTGAGAACGGTGCTGTACCTGTGGTAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTAc                                              <  2:2242079/150‑1 (MQ=38)
                                                                                                            aTGCGCCGGCCCCTTCGCGGGCGTGCCCGCGCCCCCAGTATCACCACAGACCGAGAGGACAGTGCTGTACCCGTGGGAGCCGGCCAGCCCGCGAAGAGGCCGGAACACGCCCCTGCAGAAAAAGCGCATGCGCAGAATAACacaaacaaa                                          >  2:4665965/1‑149 (MQ=255)
                                                                                                                         ttgggGGGTTTGCCGGCACCCACAGTAACACCAATCACTTTGATAAGTGTGCTGTACCTGTGGTAGCGGGCACGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCa                             >  1:4665965‑M1/6‑150 (MQ=255)
                                                                                                                               GGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCt                                                                                                                              <  2:2838015‑M1/47‑1 (MQ=255)
                                                                                                                               GGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCt                                                                                                                              >  1:2838015‑M1/1‑47 (MQ=255)
                                                                                                                                  TTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAg                    <  1:4811148‑M1/150‑1 (MQ=255)
                                                                                                                                  TTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAg                    <  2:4388916‑M1/150‑1 (MQ=255)
                                                                                                                                  TTGCCCGCTCCCACAGGTACAGCACAGTTCTCAAGGTCTGTGGTGATACTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAg                    <  1:1960917‑M1/150‑1 (MQ=255)
                                                                                                                                     CCCGCTCCCACTGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTc                 >  2:3843710‑M1/1‑150 (MQ=255)
                                                                                                                                      CCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa                <  2:502524‑M1/150‑1 (MQ=255)
                                                                                                                                      CCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa                <  2:668329‑M1/150‑1 (MQ=255)
                                                                                                                                      CCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa                <  2:2560001‑M1/150‑1 (MQ=255)
                                                                                                                                      CCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa                <  2:3958374‑M1/150‑1 (MQ=255)
                                                                                                                                        GCTCCAACTTTATCACCACAGACTTTGATAAGTTTGCTTTACGTTTGGGAGGGGCAAACCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAac              >  1:4722621‑M1/1‑150 (MQ=255)
                                                                                                                                          TCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAacac            <  1:213609‑M1/150‑1 (MQ=255)
                                                                                                                                            CCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc          <  1:2516996‑M1/150‑1 (MQ=255)
                                                                                                                                            CCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc          <  1:2231698‑M1/150‑1 (MQ=255)
                                                                                                                                            CCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc          <  2:2153758‑M1/150‑1 (MQ=255)
                                                                                                                                             CACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGcc         <  2:4424607‑M1/150‑1 (MQ=255)
                                                                                                                                             CACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGcc         <  2:3766308‑M1/150‑1 (MQ=255)
                                                                                                                                             CACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGcc         <  1:1645015‑M1/150‑1 (MQ=255)
                                                                                                                                              ACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGccc        <  2:2959677‑M1/150‑1 (MQ=255)
                                                                                                                                                 GTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAAc     >  1:3766308‑M1/1‑150 (MQ=255)
                                                                                                                                                  tATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACg    >  2:686080/1‑150 (MQ=255)
                                                                                                                                                    tCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAg  >  2:1202105/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCAGGTAGACCTGGAAGCCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGTATCACCACAGACCTTGAGAACTGTGCTGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAG  >  NC_002947/5164420‑5164717

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: