Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,839,742 T→A 7.1% intergenic (+63/+40) yheT → / ← sseA putative hydrolase/3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,839,7420TA 7.1% 405.8 / 16.4 128intergenic (+63/+40)yheT/sseAputative hydrolase/3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (59/59);  new base A (4/5);  total (63/65)
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.65e-01

Only 100 of 130 total aligned reads displayed.

TTTGTCGATGGTAGCCTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAG  >  NC_002947/5839602‑5839886
                                                                                                                                            |                                                                                                                                                
tttGTCGATGGTAGCCTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCAATGCCGGTgt                                                                                                                                         <  1:362310/150‑2 (MQ=255)
      gATGGTAGCCTGCGAAATCCGGGGTACGGCTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTGGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCGCCCTGTTGGAGCGGGCGTGCCGGCGAACACCGGCATTGCGCGTGccaccca                                                                                                                                   >  1:1345690/1‑150 (MQ=255)
        tGGTAGCCTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCcacccacc                                                                                                                                 >  1:1683750/1‑150 (MQ=255)
        tGGTAGCCTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCcacccacc                                                                                                                                 >  1:1683829/1‑150 (MQ=255)
             gCCTGCGAAACTCGGGGTAATACTTGGAACGACGGATTCCGCAGTGGCTGGTAGACGGTAGCTAGCCCGCGCATTGGTAAACGCTGCGCGCCCCAAGTGGGAGCGGGCATGCCCGCGAACACCGGCATTCCCGGTGCCACCCACCGCggg                                                                                                                            <  2:1345690/150‑1 (MQ=38)
              ccTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCgggg                                                                                                                           <  2:1683829/150‑1 (MQ=255)
               cTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCggggg                                                                                                                          >  2:2359440/1‑150 (MQ=255)
                tGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGt                                                                                                                         <  2:1778986/150‑1 (MQ=255)
                          ggggTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGt                                                                                                               >  2:1635027/1‑150 (MQ=255)
                          ggggTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGt                                                                                                               <  1:581321/150‑1 (MQ=255)
                           gggTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTc                                                                                                              <  2:2370047/150‑1 (MQ=255)
                               actactTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCt                                                                                                          <  1:2270377/150‑1 (MQ=255)
                                 tacaTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTa                                                                                                        <  1:1612041/150‑1 (MQ=255)
                                   ctTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACt                                                                                                      >  2:1586460/1‑150 (MQ=255)
                                    ttGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCTCGGGGGTTCAACCCGGTCTGCTTACTc                                                                                                     >  1:386270/1‑150 (MQ=255)
                                    ttGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTc                                                                                                     >  2:330558/1‑150 (MQ=255)
                                          cgacgGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGc                                                                                                                                               <  1:1926424/102‑1 (MQ=255)
                                          cgacgGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGc                                                                                                                                               >  2:1926424/1‑102 (MQ=255)
                                           gacgGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGc                                                                                                                                               >  2:1914963/1‑101 (MQ=255)
                                           gacgGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGc                                                                                                           >  1:735730/1‑137 (MQ=255)
                                            acgGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTgg                                                                                             >  1:1952189/1‑150 (MQ=255)
                                             cgGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGc                                                                                            >  2:478565/1‑150 (MQ=255)
                                              ggATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGcc                                                                                           >  1:1926577/1‑150 (MQ=255)
                                                 ttCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCAcc                                                                                        >  1:1477381/1‑150 (MQ=255)
                                                    cGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTg                                                                                              <  1:486193/141‑1 (MQ=255)
                                                    cGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCa                                                                                     >  2:612692/1‑150 (MQ=255)
                                                     gCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCAt                                                                                    <  2:1926577/150‑1 (MQ=255)
                                                          ggCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCAtgctgc                                                                               <  2:978984/150‑1 (MQ=255)
                                                           gCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGc                                                                                                           <  2:1875173/121‑1 (MQ=255)
                                                           gCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGc                                                                                                           >  1:1875173/1‑121 (MQ=255)
                                                           gCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCt                                                                                                          <  1:2322640/122‑1 (MQ=255)
                                                             tggtggACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCggg                                                                            >  2:1366899/1‑150 (MQ=255)
                                                               gtggACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGgtt                                                                          <  1:1785694/150‑1 (MQ=255)
                                                               gtggACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGCGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGgtt                                                                          <  1:2329647/150‑1 (MQ=255)
                                                                tggACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGgttg                                                                         <  2:559812/150‑1 (MQ=255)
                                                                    cGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGa                                                                     <  1:2359440/150‑1 (MQ=255)
                                                                     ggTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGAt                                                                    >  1:256899/1‑150 (MQ=255)
                                                                       tagctagcCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCg                                                                              >  1:2147430/1‑138 (MQ=255)
                                                                       tagctagcCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCg                                                                              <  2:2147430/138‑1 (MQ=255)
                                                                        agctagcCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCca                                                                 <  2:7335/150‑1 (MQ=255)
                                                                          ctagcCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCcact                                                               >  2:931509/1‑150 (MQ=255)
                                                                              cccGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGACCACCGGCATTGCCGGTGCCACCCACCGCGGGGGtt                                                                                                                        <  1:2371844/89‑1 (MQ=255)
                                                                              cccGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGtt                                                                                                                        >  2:2371844/1‑89 (MQ=255)
                                                                               ccGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGc                                                                                                           >  2:496543/1‑101 (MQ=255)
                                                                               ccGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGc                                                                                                           <  1:496543/101‑1 (MQ=255)
                                                                               ccGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCggg                                                                            <  2:419660/132‑1 (MQ=255)
                                                                               ccGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCactcactc                                                          >  1:682088/1‑150 (MQ=255)
                                                                               ccGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCactcactc                                                          <  1:847159/150‑1 (MQ=255)
                                                                               ccGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACAGCATGCTGCggg                                                                            >  1:419660/1‑132 (MQ=255)
                                                                                cgcgAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTcc                                                         <  1:478565/150‑1 (MQ=255)
                                                                                cgcgAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTcc                                                         >  2:814925/1‑150 (MQ=255)
                                                                                         gTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCg                                                <  1:1428781/150‑1 (MQ=255)
                                                                                         gTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCg                                                >  2:1331705/1‑150 (MQ=255)
                                                                                          tAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCgg                                                                                                                             >  1:1459020/1‑72 (MQ=255)
                                                                                          tAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCgg                                                                                                                             <  2:1459020/72‑1 (MQ=255)
                                                                                          tAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCAc                                                       >  2:1783566/1‑142 (MQ=255)
                                                                                          tAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCAc                                                       <  1:1783566/142‑1 (MQ=255)
                                                                                               gCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTa                                                                                                        <  1:1838235/88‑1 (MQ=255)
                                                                                               gCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACa                                          >  2:2095767/1‑150 (MQ=255)
                                                                                                    gcgcCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTg                                     >  1:1282254/1‑150 (MQ=255)
                                                                                                     cgcCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTTCACGACCCCGCATACAGCTTGc                                    <  2:2021890/150‑1 (MQ=255)
                                                                                                     cgcCCCCTGGGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGc                                    <  2:566293/150‑1 (MQ=255)
                                                                                                      gcCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACACCTTGcc                                   >  1:2214241/1‑150 (MQ=255)
                                                                                                       cccccTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGc                                                                                                   <  1:652582/85‑1 (MQ=255)
                                                                                                       cccccTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGc                                                                                                   >  2:652582/1‑85 (MQ=255)
                                                                                                         cccTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCgg                                                                                                                             <  1:20442/57‑1 (MQ=255)
                                                                                                         cccTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCgg                                                                                                                             >  2:20442/1‑57 (MQ=255)
                                                                                                          ccTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGc                               <  1:1635027/150‑1 (MQ=255)
                                                                                                           cTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCcactcac                                                            <  1:2409291/120‑1 (MQ=255)
                                                                                                           cTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCcactcac                                                            >  2:2409291/1‑120 (MQ=255)
                                                                                                            tgtgGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACa                                          <  2:1899920/137‑1 (MQ=255)
                                                                                                            tgtgGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACa                                          >  1:1899920/1‑137 (MQ=255)
                                                                                                                ggAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTg                                     <  1:2269880/138‑1 (MQ=255)
                                                                                                                ggAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTg                                     >  2:2269880/1‑138 (MQ=255)
                                                                                                                ggAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTAc                         <  2:1703436/150‑1 (MQ=255)
                                                                                                                ggAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCAGTTCAGGCCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGGCCAGCGGGTAc                         >  1:1354103/1‑150 (MQ=38)
                                                                                                                          tGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCcgccagcg               >  2:1318085/1‑150 (MQ=255)
                                                                                                                            gacacgAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCcgccagcgcc             >  2:33169/5‑150 (MQ=255)
                                                                                                                            cccGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTg                                                                                              <  1:150972/69‑1 (MQ=255)
                                                                                                                             ccGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGAcccc                                                 <  2:726085/113‑1 (MQ=255)
                                                                                                                             ccGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGAcccc                                                 >  1:726085/1‑113 (MQ=255)
                                                                                                                             ccGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCcgccagcgcca            >  2:841386/1‑150 (MQ=255)
                                                                                                                              cgcgAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGAccc                                                  >  1:2132390/1‑111 (MQ=255)
                                                                                                                               gcgAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTc                                                                                                              <  2:2037247/50‑1 (MQ=255)
                                                                                                                                   gcacCGGCAATGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCg                                                                                                 <  2:96756/58‑1 (MQ=255)
                                                                                                                                   gcacCGGCAATGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCactcactc                                                          >  2:1515598/2‑98 (MQ=255)
                                                                                                                                   gcacCGGCAATGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCactcactc                                                          <  1:1515598/97‑1 (MQ=255)
                                                                                                                                   acacCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCaa      >  1:1809122/1‑150 (MQ=255)
                                                                                                                                      ccGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACa   <  1:1318085/150‑1 (MQ=255)
                                                                                                                                      ccGGCAATGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCggg                                                                            >  2:1485016/1‑77 (MQ=255)
                                                                                                                                      ccGGCAATGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCggg                                                                            >  2:1485498/1‑77 (MQ=255)
                                                                                                                                      ccGGCAATGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCggg                                                                            <  1:1485498/77‑1 (MQ=255)
                                                                                                                                       cGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGgttg                                                                         >  1:258878/1‑79 (MQ=255)
                                                                                                                                       cGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGgttg                                                                         <  2:258878/79‑1 (MQ=255)
                                                                                                                                       cGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCcgccagc                <  1:783664/136‑1 (MQ=255)
                                                                                                                                       cGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCcgccagc                >  2:783664/1‑136 (MQ=255)
                                                                                                                                       cGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAg  <  2:434027/150‑1 (MQ=255)
                                                                                                                                           aTTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTAcccc                      <  2:231972/126‑1 (MQ=255)
                                                                                                                                            ttGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCAtgctgc                                                                               >  1:1931194/1‑68 (MQ=255)
                                                                                                                                            ttGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCAtgctgc                                                                               <  2:1931194/68‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
TTTGTCGATGGTAGCCTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAG  >  NC_002947/5839602‑5839886

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