Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,637,439 G→A 6.4% intergenic (+4/‑28) PP_4948 → / → PP_4949 acyl‑CoA dehydrogenase family protein/TldD/PmbA family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,637,4390GA 6.4% 185.3 / 2.1 62intergenic (+4/‑28)PP_4948/PP_4949acyl‑CoA dehydrogenase family protein/TldD/PmbA family peptidase
Reads supporting (aligned to +/- strand):  ref base G (32/26);  new base A (2/2);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

TGCTGATCGCGACTTCTACCTCGGCAAGCTGCAAGCTGCGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCGGTTCGCCGCAACGTC  >  NC_002947/5637294‑5637586
                                                                                                                                                  |                                                                                                                                                   
tgctgATCGCGACTTCTACCTCGGCAAGCTGCAAGCTGCGCGTTATTTCCTGACATGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGTTGGTTCTAGTGGGAGc                                                                                                                                                  >  7:48136/1‑149 (MQ=255)
  ctgATCGCGACTTCTACCTCGGCAAGCTGCAAGCTGCGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGcc                                                                                                                                                 <  4:129773/148‑1 (MQ=255)
     aTCGCGACTTCTACCGCGCCACGCTGCAAGCTGCGGGTTATTTCCTTCCCTTGGAAGTACCGGCCTGCCAAATTCAGCTGCCATTGCGCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCAcc                                                                                                                                              <  8:269504/148‑1 (MQ=255)
      tCGCGACTTCTACCTCGGCAAGCTGCAAGCTGCGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCgt                                                                                                                                            <  5:49580/149‑1 (MQ=255)
      tCGCGACTTCTACCTCGGCAAGCTGCAAGCTGCGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCgt                                                                                                                                            >  2:36292/1‑149 (MQ=255)
                tACCTCGGCAAGCTGCAAGCTGCGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAg                                                                                                                                  >  4:319762/1‑149 (MQ=255)
                  ccTCGGCATGGTGCAAGCTGCGCGTTATTTCCTGACCTGGGAATTACCGGTCTTCCATATTTAGCTGGCATTGCTCGAGTCACACGACGATACCTTCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCTCCGTGTACACGGTGGc                                                                                                                                <  6:237806/149‑1 (MQ=255)
                  ccTCGGCAAGCTGCAAGCTGCGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGa                                                                                                                                   <  6:101221/146‑1 (MQ=255)
                  ccTCGGCAAGCTGCAAGCTGCGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAg                                                                                                                                  <  3:319762/147‑1 (MQ=255)
                   cTCGGCAAGCTGCAAGCTGCGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCa                                                                                                                               >  3:207223/1‑149 (MQ=255)
                               caagctgcGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTc                                                                                                                    >  6:107560/1‑148 (MQ=255)
                                aagctgcGCGTTATTTCCTTACCTGCCACTTCCCGGGCTGCAAAAATGAGCTGTCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGa                                                                                                                  <  8:48136/149‑1 (MQ=255)
                                  gctgcGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCTCCGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGAt                                                                                                                 <  1:263427/149‑1 (MQ=255)
                                  gctgcGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGAtt                                                                                                                >  5:138220/1‑149 (MQ=255)
                                      cgcgTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAg                                                                                                            >  4:298570/1‑149 (MQ=255)
                                       gcgTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAgc                                                                                                           >  6:311478/1‑149 (MQ=255)
                                                ccTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCctgctg                                                                                                  >  6:52934/1‑149 (MQ=255)
                                                      cTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAg                                                                                            >  4:9786/1‑149 (MQ=255)
                                                           aaGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCtt                                                                                       <  5:107560/149‑1 (MQ=255)
                                                              tACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTtgc                                                                                    <  5:311478/149‑1 (MQ=255)
                                                               aCCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACATGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACAAGGAGGCATcc                                                                                                                            >  1:323810/1‑108 (MQ=255)
                                                                    gCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTgctgcgct                                                                              >  8:31134/1‑149 (MQ=255)
                                                                        ccATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAgc                                                                                                           <  7:33431/116‑1 (MQ=255)
                                                                        ccATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAgc                                                                                                           >  8:33431/1‑116 (MQ=255)
                                                                        ccATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCtgcgctgcgc                                                                          <  1:124452/149‑1 (MQ=255)
                                                                        ccATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGAAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGctgcgctgcg                                                                           >  7:285924/1‑148 (MQ=255)
                                                                                  cTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGcc                                                                 >  4:310881/1‑148 (MQ=255)
                                                                                       aTTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTgctgcgctg                                                                             >  8:26758/1‑131 (MQ=255)
                                                                                       aTTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTgctgcgctg                                                                             <  7:26758/131‑1 (MQ=255)
                                                                                       aTTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAgttg                                                           >  6:119336/1‑149 (MQ=255)
                                                                                       aTTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAgttg                                                           >  3:141093/1‑149 (MQ=255)
                                                                                       aTTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAggtg                                                           >  8:203336/1‑149 (MQ=255)
                                                                                          gCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAgttgtt                                                         <  7:31134/148‑1 (MQ=255)
                                                                                           cTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGAATCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGtttt                                                       >  6:160493/1‑149 (MQ=255)
                                                                                             cGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGAAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTc                                                      <  8:285924/148‑1 (MQ=255)
                                                                                              gAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTccc                                                    >  6:283680/1‑149 (MQ=255)
                                                                                              gAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTccc                                                    >  6:117073/1‑149 (MQ=255)
                                                                                                   acacGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGAAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTgcgc                                               >  2:181240/1‑149 (MQ=255)
                                                                                                       gacgaTACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCttt                                                                                      >  7:128627/1‑106 (MQ=255)
                                                                                                       gacgaTACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCttt                                                                                      <  8:128627/106‑1 (MQ=255)
                                                                                                       gacgaTACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCAtg                                           >  5:177172/1‑149 (MQ=255)
                                                                                                        acgaTACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCAtgt                                          <  5:52934/149‑1 (MQ=255)
                                                                                                            tACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGaaa                                      >  1:320297/1‑149 (MQ=255)
                                                                                                            tACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGaaa                                      <  4:141093/149‑1 (MQ=255)
                                                                                                             aCCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCt                                                   >  3:145775/1‑135 (MQ=255)
                                                                                                             aCCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCt                                                   <  4:145775/135‑1 (MQ=255)
                                                                                                             aCCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGaaa                                      >  5:201661/1‑148 (MQ=255)
                                                                                                                 gCCT‑CGGCATGCAAGACGAGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTc                                 >  2:216496/1‑149 (MQ=255)
                                                                                                                  ccT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTc                                 <  3:310881/148‑1 (MQ=255)
                                                                                                                    t‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCa                              >  2:79440/1‑149 (MQ=255)
                                                                                                                         cATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAg                          <  4:207223/149‑1 (MQ=255)
                                                                                                                             cAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGgcgc                      >  6:273658/1‑149 (MQ=255)
                                                                                                                              aaGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCt                     >  8:175552/1‑149 (MQ=255)
                                                                                                                                gACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCttt                   <  6:177172/149‑1 (MQ=255)
                                                                                                                                 aCGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTc                  <  5:273658/149‑1 (MQ=255)
                                                                                                                                 aCGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTc                  <  6:138220/149‑1 (MQ=255)
                                                                                                                                  cGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCg                 <  3:298570/149‑1 (MQ=255)
                                                                                                                                  cGGGTGGTTCTAGGGGAAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCg                 <  1:181240/149‑1 (MQ=255)
                                                                                                                                   gggTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCgg                >  8:122510/1‑149 (MQ=255)
                                                                                                                                        gtTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCGGTTcgc           <  3:9786/149‑1 (MQ=255)
                                                                                                                                            tAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCttt                   >  4:59867/1‑137 (MQ=255)
                                                                                                                                            tAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCttt                   <  3:59867/137‑1 (MQ=255)
                                                                                                                                                 ggAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCGGTTCGCCGCAACGTc  >  4:121761/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGCTGATCGCGACTTCTACCTCGGCAAGCTGCAAGCTGCGCGTTATTTCCTGACCTGGGAAGTACCGGGCTGCCATAATGAGCTGGCATTGCTCGAGGCACACGACGATACCTGCCT‑CGGCATGCAAGACGGGTGGTTCTAGGGGGAGCCACCGTGTACACGGAGGCATCCGCATGTTCGATTCCAGCGCCCTGCTGCGCCAGCGCTTTGCTGCGCTGCGCAGTACGGCCGAGTTGTTTTCCCTGCGCCATGTGAAACAGTCGCACCAGGCGCTTTCGGTTCGCCGCAACGTC  >  NC_002947/5637294‑5637586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑