Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,158 C→T 23.0% intergenic (+109/‑204) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1580CT23.0% 269.5 / 82.7 157intergenic (+109/‑204)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (64/55);  new base T (18/16);  total (84/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

Only 100 of 158 total aligned reads displayed.

CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCCGCCA  >  NC_002947/955010‑955307
                                                                                                                                                     |                                                                                                                                                       
ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCAACAAGGTTCCt                                                                                                                                                        >  2:3285569/1‑150 (MQ=255)
 cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑g                                                                                                                                                      <  2:582778/150‑1 (MQ=255)
 cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑g                                                                                                                                                      <  1:3915083/150‑1 (MQ=255)
      cGCTGCGGATGTGCTGATGACCTACGCACGCTGAACCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGa                                                                                                                                                 <  2:6525622/150‑1 (MQ=255)
      cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGa                                                                                                                                                 <  2:6525733/150‑1 (MQ=255)
      cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGa                                                                                                                                                 <  2:3504839/150‑1 (MQ=255)
       gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTT‑CTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTT‑GTAGGAG‑c                                                                                                                                              <  1:4823756/150‑1 (MQ=25)
       gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGG‑TTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCTCCAACAAGGTTCCT‑GTAGGAG‑c                                                                                                                                              >  2:5947353/1‑150 (MQ=255)
        cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTT‑GTAGGAG‑c                                                                                                                                              <  2:2973896/150‑1 (MQ=38)
        cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTT‑GTAGGAG‑c                                                                                                                                              <  1:6189638/150‑1 (MQ=38)
        cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑c                                                                                                                                              <  2:3192945/150‑1 (MQ=255)
         tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTT‑GTAGGAG‑Cg                                                                                                                                             <  2:1036228/150‑1 (MQ=38)
          gCGTATGTGCTTATGACCTACGCACGATGAAACTGCTGCCTGTTTCAAGAAAACCCCGGCCACTGTGCCTGGTTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑Cgt                                                                                                                                            <  1:6536110/150‑2 (MQ=38)
          gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGCACCTT‑GTAGGAG‑Cgg                                                                                                                                            <  1:5598469/150‑1 (MQ=38)
          gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTT‑GTAGGAG‑Cgg                                                                                                                                            <  2:2534014/150‑1 (MQ=38)
          gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTT‑GTAGGAG‑Cgg                                                                                                                                            <  2:1512091/150‑1 (MQ=38)
          gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTT‑GTAGGAG‑Cgg                                                                                                                                            <  2:5087240/150‑1 (MQ=38)
          gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑Cgg                                                                                                                                            <  2:6543070/150‑1 (MQ=255)
          gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑Cgg                                                                                                                                            <  2:1675872/150‑1 (MQ=255)
          gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCGGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTT‑GTAGGAG‑Cgg                                                                                                                                            <  2:6464889/150‑1 (MQ=37)
          gCGGAGGTGCTGATGACCTACGCACGCTGAAACTGCGGCCTGTTTCACGAAACCCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑Cgg                                                                                                                                            <  1:3662330/150‑1 (MQ=255)
           cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑Cggg                                                                                                                                           <  1:3412775/150‑1 (MQ=255)
           cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTT‑CTGTGCTGGCCTCTTCACGGGTAAACCCGCTCCTACAGGAACCTT‑GTAGGAG‑CGGGt                                                                                                                                          <  1:2147007/150‑1 (MQ=25)
              aTGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGttt                                                                                                                                        <  1:2108466/150‑1 (MQ=255)
                gtgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTT‑GTAGGAG‑CGGGtggtc                                                                                                                                      >  1:2194068/1‑146 (MQ=37)
                gtgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTAc                                                                                                                                      <  2:755361/150‑1 (MQ=255)
                     gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGGCTCTTCGCGGGTAAAACCCCACCGACAAGGTTCGG‑GGTGGGG‑CGGGTTTTGCcgcg                                                                                                                                 >  1:446686/1‑150 (MQ=255)
                       tgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTTTCCT‑GTAGGAG‑CGGGTTTACCCGCGaa                                                                                                                               <  2:679372/150‑1 (MQ=255)
                         aCCTCCGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACGGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTGCCTGTGCCGGCCTCTTCGCGGGGAACCCCGCTCCTACCAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAaga                                                                                                                             <  2:5530742/150‑1 (MQ=38)
                            tACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCGACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGc                                                                                                                          >  2:5873832/1‑150 (MQ=255)
                              cgccggcTGAACCTGCGGCCTGTTACACGAACCCCCCGGCCACGGTGCGGGGGTTTTTTTTG‑GGGGGTTTTTTTCGTGGGCGGGCGTCTTCGCGGGTAAACCGGCTCC‑CCCCGGACCCTGGTCGGCG‑CGGGTTTACCCGCGAAGAGGccg                                                                                                                        <  2:6403431/145‑2 (MQ=12)
                              cccccgcTGAACCTGCGGCCTGTGTCACGAAAACCCCGGCCCCTGTGCGGGGGATTTTTTTT‑GGGGGTATTTTACGTGGGCCGGCCTCTTCGCGGGTACCCCCGCTCCGACAAGGGACCG‑GTCGGAG‑CGGGTTTACCCGCGAAGAGGccg                                                                                                                        <  2:4151939/146‑2 (MQ=21)
                                     tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCCAGGTTCCT‑GTAGGAG‑CGGGGTTACCCGCGAAGAGGCCAGCACAGa                                                                                                                 >  2:5551960/1‑150 (MQ=255)
                                     tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAg                                                                                                                                                   <  2:1971994/117‑1 (MQ=255)
                                                    tttCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCCGCACAGAAAAATAAAAATAccg                                                                                                  >  2:5883996/1‑149 (MQ=255)
                                                        aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGAtt                                                                                              >  2:3355951/1‑150 (MQ=255)
                                                         cGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTAGCGGGTTAAACCCCGCCGACCAGGTTGCG‑GGGGGCGCCGGGTTTA‑CCGCCAAGAGGGCCGCCCCGCAAAATAAAAATAACGGAttt                                                                                             >  2:2314391/1‑150 (MQ=255)
                                                           aaaaCCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑Cggg                                                                                                                                           >  1:552738/1‑102 (MQ=255)
                                                           aaaaCCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑Cggg                                                                                                                                           <  2:552738/102‑1 (MQ=255)
                                                            accccccGCCCAC‑GGGCCGGGGTTTTTTTTTGGGGGGTATTTTCCCGGGGCCGCCCCCGCCGGGGGAAACCCCCCCCCAACACGGTCCCT‑GCTGGCG‑CGGGTTCACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCa                                                                                          <  1:1469556/147‑1 (MQ=255)
                                                                   ggCCACTGTGCGGGGTTTTTTTTTT‑GGGGTTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATa                                                                                   <  1:6217403/150‑1 (MQ=255)
                                                                    gCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCGCCAACAAGGTTGCG‑GGAGGAG‑CGGGGTTACCCCCGCCGAGGCCCGCACACAAAAAAAAAAAAAACAGATTTTCAAAAAATAg                                                                                  >  2:2543344/1‑150 (MQ=255)
                                                                      cACGGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt                                                                                <  2:5653458/150‑1 (MQ=255)
                                                                                ggtttttttttt‑GGGGTTATTTTGCCTGGGCCGGCCTCTGCGGGGGAAACCCCGCTCCTACAGGGTTCCT‑GTTGTTG‑GGTTTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGtt                                                                      <  1:4330097/147‑1 (MQ=255)
                                                                                  gtattttttt‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGAGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCt                                                                    <  2:6778838/147‑1 (MQ=255)
                                                                                  gTTTTTTTTT‑GGGGGTATTTTTCCTGGGCCGGCCTCTGCGCGGGTAAACCCGCTCCTCCACGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCt                                                                    <  2:4879277/150‑1 (MQ=255)
                                                                                    tttttttt‑GGGGTTTTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa                                                                  <  1:2226506/150‑1 (MQ=255)
                                                                                           t‑GGGGTTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACAGGGTTCGT‑GTTGGAG‑GGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                           <  1:2600705/150‑1 (MQ=38)
                                                                                           t‑GGGGGTATTTTACCTGTGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACCAGGTTCCT‑GTTGGCG‑CGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                           <  1:3285569/150‑1 (MQ=38)
                                                                                                ggTATTTTTCCTGTGCCGGCCTCTTCGCGGTTAAACCCGATCCTACAAGGTTCCT‑GTAGGTG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaa                                                       <  1:2697009/150‑1 (MQ=255)
                                                                                                             tgCCGGCATCTTGGCGGTAAACCCCGCTCCTCAATGGTTCCT‑GTTGGTG‑CGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa                                          <  1:1496866/150‑1 (MQ=255)
                                                                                                             ggCCGGCCCCTGCGCGGGTAAACCCGCACCAACCAGGTTCCT‑GTTGTAG‑CGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa                                          <  1:442749/149‑1 (MQ=255)
                                                                                                                  gCCTCTTCGCGGGTAAAACCGCTCCTACAAGGTTCCG‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCAGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATACTTTACAGGCAt                                     >  2:1594386/1‑150 (MQ=255)
                                                                                                                       ttCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAAAAAAACAGGAAATTATCCTTTACAGGCAGCCTGt                                >  2:5240903/1‑150 (MQ=255)
                                                                                                                         cgcgGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACCGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCGAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                              >  2:3979732/1‑150 (MQ=255)
                                                                                                                             ggTAAACCCGCTCCTACACGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                          <  1:3920698/150‑1 (MQ=255)
                                                                                                                             ggTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                          >  2:3499869/1‑150 (MQ=255)
                                                                                                                               tAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc                        >  2:2794039/1‑150 (MQ=255)
                                                                                                                                aaccccGCTCCTCCACGTTTCGT‑GTAGGCG‑CGGGTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACt                       <  1:3549158/147‑1 (MQ=39)
                                                                                                                                  aCCCGCTCCTACAGGAACCTT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTGCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt                     >  1:166320/1‑150 (MQ=37)
                                                                                                                                  aCCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt                     >  2:4017283/1‑150 (MQ=255)
                                                                                                                                   cccGCTCCTACCAGGTTCCT‑GTAGGAG‑CGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                    <  2:4249873/150‑1 (MQ=255)
                                                                                                                                   cccGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                    >  1:6736174/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCTACAGGAACCTT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  2:4430716/1‑150 (MQ=38)
                                                                                                                                    ccGCTCCTACAGGAACCTT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  1:4539755/1‑150 (MQ=38)
                                                                                                                                    ccGCTCCTACAGGAACCTT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  1:1700782/1‑150 (MQ=38)
                                                                                                                                    ccGCTCCTACAGGAACCTT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  2:2678149/1‑150 (MQ=38)
                                                                                                                                    ccGCTCCTACAGGAACCTT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  1:4388752/1‑150 (MQ=38)
                                                                                                                                    ccGCTCCTACAGGAACCTT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  1:3018178/1‑150 (MQ=38)
                                                                                                                                    ccGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  2:684376/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTGTACAGGCATCCTGTTTCGGCACTAtct                   >  2:3055127/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATAATCCTTTACAGGCATCCTGTTTCGGCACTAtct                   >  2:1280899/1‑150 (MQ=255)
                                                                                                                                     cGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt                  >  2:3465196/1‑150 (MQ=255)
                                                                                                                                      gCTCCTCCTAGGTCCCT‑GTTGGAG‑GGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc                 <  1:1820906/150‑1 (MQ=255)
                                                                                                                                      gCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                             <  2:5584659/138‑1 (MQ=255)
                                                                                                                                      gCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                             >  1:5584659/1‑138 (MQ=255)
                                                                                                                                       cTCCTACACGGTTCCT‑GTCGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttct                <  1:3303943/150‑1 (MQ=255)
                                                                                                                                        tCCTACCAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTg               >  2:2463691/1‑150 (MQ=255)
                                                                                                                                        tCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTg               >  1:3987468/1‑150 (MQ=255)
                                                                                                                                         ccTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa              <  1:2794039/150‑1 (MQ=255)
                                                                                                                                          cTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc                 >  2:2161469/1‑146 (MQ=255)
                                                                                                                                           tACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt            >  2:6537817/1‑150 (MQ=255)
                                                                                                                                           tACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt            >  2:6164745/1‑150 (MQ=255)
                                                                                                                                           tACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt            >  2:785678/1‑150 (MQ=255)
                                                                                                                                           tACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt            >  2:901948/1‑150 (MQ=255)
                                                                                                                                           tACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt            >  2:5320652/1‑150 (MQ=255)
                                                                                                                                           tACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt            >  2:3075817/1‑150 (MQ=255)
                                                                                                                                           tACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt            >  1:5126826/1‑150 (MQ=255)
                                                                                                                                           tACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGG‑CAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc           >  1:5674881/1‑150 (MQ=255)
                                                                                                                                             cAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa          >  2:47857/1‑150 (MQ=255)
                                                                                                                                             cAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa          >  1:104470/1‑150 (MQ=255)
                                                                                                                                              aaGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa          <  1:5789330/149‑1 (MQ=255)
                                                                                                                                              aaGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg         <  2:1828402/150‑1 (MQ=255)
                                                                                                                                              aaGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg         >  2:1783190/1‑150 (MQ=255)
                                                                                                                                              aaGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg         >  1:5739071/1‑150 (MQ=255)
                                                                                                                                              aaGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg         >  1:3118461/1‑150 (MQ=255)
                                                                                                                                                ggTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgc       <  1:5067785/150‑1 (MQ=255)
                                                                                                                                                 gTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgcc      >  2:4055794/1‑150 (MQ=255)
                                                                                                                                                  ttCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa                                                      >  1:1407308/1‑101 (MQ=255)
                                                                                                                                                     cT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCCGCca  >  2:1475377/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTT‑GGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCT‑GTAGGAG‑CGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCCGCCA  >  NC_002947/955010‑955307

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: ‑