Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,393 T→A 6.4% intergenic (+13/+129) PP_3812 → / ← PP_5632 lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,338,3930TA 6.4% 223.9 / 3.2 72intergenic (+13/+129)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base A (3/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTC  >  NC_002947/4338251‑4338539
                                                                                                                                                 |                                                                                                                                                   
gACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGAccc                                                                                                                                               >  6:45765/1‑149 (MQ=255)
gACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGAccc                                                                                                                                               >  6:53073/1‑149 (MQ=255)
gACAGTACCTCGGTGCGCAACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGGTGCGCCAAGTGCCGGCAGTAAACGCCGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATGTGCAGTGCGGCGATGAGCAATAGAcccc                                                                                                                                              >  6:122523/1‑148 (MQ=255)
   aGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGcc                                                                                                                                            >  3:202249/1‑149 (MQ=255)
         cTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCg                                                                                                                                      >  2:253670/1‑149 (MQ=255)
            ggTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCggg                                                                                                                                    <  7:94305/148‑1 (MQ=255)
            ggTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCgg                                                                                                                                     <  7:231834/147‑1 (MQ=255)
                 gcTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGTAGATGGGGCAGGCGGGCAACGAAATTAAATTGAAATTGCATTGAGGGGATTACCAAAAAACCCCCCTGGGCgggggt                                                                                                                                 >  1:292001/1‑146 (MQ=255)
                  cTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTg                                                                                                                             <  8:333691/149‑1 (MQ=255)
                  cTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTg                                                                                                                             <  5:53073/149‑1 (MQ=255)
                  cTACAACGGAATGATCCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAAAAGACCCGCCCAGGCGGGTCTATTg                                                                                                                             <  5:277846/149‑1 (MQ=255)
                   tACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGc                                                                                                                            >  2:9519/1‑149 (MQ=255)
                               aTTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCgg                                                                                                                <  5:122523/149‑1 (MQ=255)
                                 tCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGa                                                                                                               >  3:217678/1‑148 (MQ=255)
                                   aTGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGaaa                                                                                                             >  7:197979/1‑148 (MQ=255)
                                   aTGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGaaa                                                                                                             >  6:328181/1‑148 (MQ=255)
                                    tGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACa                                                                                                           <  2:361264/149‑1 (MQ=255)
                                      aCTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGc                                                                                                         <  5:328181/149‑1 (MQ=255)
                                         gcacttggctCAACGGGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCg                                                                                                      >  6:368472/8‑149 (MQ=255)
                                             gctgctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAgg                                                                                                  >  2:339973/1‑149 (MQ=255)
                                              ctgctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGc                                                                                                 <  5:45765/149‑1 (MQ=255)
                                                 ctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCAGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTc                                                                                              <  3:262418/149‑1 (MQ=255)
                                                  tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAAAAGACCCGCCCAGGCGGGTCtatt                                                                                                                              >  5:145795/1‑116 (MQ=39)
                                                  tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAAAAGACCCGCCCAGGCGGGTCtatt                                                                                                                              <  6:145795/116‑1 (MQ=39)
                                                       ccGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGaa                                                                                        <  4:217678/149‑1 (MQ=255)
                                                       ccGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGaa                                                                                        <  1:9519/149‑1 (MQ=255)
                                                             aGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGc                                                                                   <  3:149948/148‑1 (MQ=255)
                                                             aGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGc                                                                                   <  3:144587/148‑1 (MQ=255)
                                                                  aGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCgg                                                                             >  8:105980/1‑149 (MQ=255)
                                                                   gTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCtttt                                                                                                                              >  6:31919/1‑99 (MQ=255)
                                                                   gTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCtttt                                                                                                                              <  5:31919/99‑1 (MQ=255)
                                                                   gTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGt                                                                            >  2:339770/1‑149 (MQ=255)
                                                                      ccGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCg                                                                         >  8:60431/1‑149 (MQ=255)
                                                                        ggCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCgg                                                                        >  4:133898/1‑148 (MQ=255)
                                                                        ggCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGc                                                                       <  7:60431/149‑1 (MQ=255)
                                                                         gCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCgg                                                                        >  2:48587/1‑147 (MQ=255)
                                                                           aGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCTGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAgg                                                                    <  4:202249/149‑1 (MQ=255)
                                                                              aaaGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTcc                                                                 >  7:339875/1‑149 (MQ=255)
                                                                                aGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCgg                                                               >  7:274723/1‑149 (MQ=255)
                                                                                aGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCgg                                                               >  1:191062/1‑149 (MQ=255)
                                                                                aGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCgg                                                               <  1:339973/149‑1 (MQ=255)
                                                                                 ggcggc‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGc                                                              >  7:316235/1‑149 (MQ=255)
                                                                                              cGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                   <  5:148724/149‑1 (MQ=255)
                                                                                              cGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                   <  8:197979/149‑1 (MQ=255)
                                                                                              cGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                   <  6:137336/149‑1 (MQ=255)
                                                                                              cGGCAGGCGGGCAACGAACTTAAAGGGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                   <  3:350922/149‑1 (MQ=255)
                                                                                                gCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGt                                                 <  7:105980/149‑1 (MQ=255)
                                                                                                      gggCAACGAACTTAGAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                   >  7:304908/1‑141 (MQ=255)
                                                                                                      gggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCgg                                                   <  8:304908/141‑1 (MQ=255)
                                                                                                       ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcg                                          <  5:4188/149‑1 (MQ=255)
                                                                                                            cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                     <  2:292001/149‑1 (MQ=255)
                                                                                                            cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                     <  8:274723/149‑1 (MQ=255)
                                                                                                            cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                     <  5:368472/149‑1 (MQ=255)
                                                                                                             gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAg                                    <  3:133898/149‑1 (MQ=255)
                                                                                                               aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                   <  3:219389/148‑1 (MQ=255)
                                                                                                               aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                   >  4:219389/1‑148 (MQ=255)
                                                                                                                 ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGcc                                >  7:67512/1‑149 (MQ=255)
                                                                                                                 ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAGCCCGCGACAggggc                                >  4:92752/1‑147 (MQ=255)
                                                                                                                  tAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTACCCGGGAAAACCGACAAAAGGccc                                >  7:300901/1‑149 (MQ=255)
                                                                                                                   aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCgg                              >  5:287541/1‑149 (MQ=255)
                                                                                                                      gTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCttt                                                         <  7:224595/120‑1 (MQ=255)
                                                                                                                      gTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCttt                                                         >  8:224595/1‑120 (MQ=255)
                                                                                                                      gTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCa                            >  5:291384/1‑148 (MQ=255)
                                                                                                                               cAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTg                     <  8:339875/146‑1 (MQ=255)
                                                                                                                                  tGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAgg               <  8:316235/149‑1 (MQ=255)
                                                                                                                                   gAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTGGCGGGGTATCCCGCGCAAGGGGCGGCAAATCCTGGTCAggg              >  5:208393/1‑149 (MQ=255)
                                                                                                                                    aGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGt                                                                            >  3:93828/1‑87 (MQ=255)
                                                                                                                                    aGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGt                                                                            <  4:93828/87‑1 (MQ=255)
                                                                                                                                    aGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGc             >  7:143205/1‑149 (MQ=255)
                                                                                                                                      gCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGTATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCt            >  2:224101/1‑148 (MQ=255)
                                                                                                                                           aaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAcc      >  3:119470/2‑149 (MQ=255)
                                                                                                                                               aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTAGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTACAACCCtta  >  4:57273/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGC‑‑‑GCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTC  >  NC_002947/4338251‑4338539

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