Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,469,860 T→C 5.5% intergenic (+29/+6) vacJ → / ← PP_2164 lipoprotein VacJ/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,469,8600TC 5.5% 371.6 / 8.0 127intergenic (+29/+6)vacJ/PP_2164lipoprotein VacJ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (64/56);  new base C (3/4);  total (67/60)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-02

Only 100 of 128 total aligned reads displayed.

GCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGAT  >  NC_002947/2469712‑2469999
                                                                                                                                                     |                                                                                                                                           
gcCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCt                                                                                                                                             <  1:155876/149‑1 (MQ=255)
   aGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCGTACGTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCgg                                                                                                                                          <  6:195664/149‑1 (MQ=255)
   aGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCgg                                                                                                                                          <  1:59994/149‑1 (MQ=255)
    gCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCggg                                                                                                                                         >  4:263973/1‑149 (MQ=255)
     ccTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAAGCTGGGGCCGCAAAGCGGCCTCGGGt                                                                                                                                        <  7:98117/149‑1 (MQ=255)
     ccTGCTCTAGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGt                                                                                                                                        <  6:250138/149‑1 (MQ=255)
        gCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGttt                                                                                                                                      >  7:50109/1‑147 (MQ=255)
        gCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGt                                                                                                                                        <  8:50109/146‑1 (MQ=255)
        gCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTc                                                                                                                                      >  7:162543/1‑148 (MQ=255)
         ctctCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGATTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGTTCa                                                                                                                                     <  1:133551/148‑1 (MQ=255)
         ctctCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGTTCa                                                                                                                                     <  7:147437/148‑1 (MQ=255)
          tctcGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAAt                                                                                                                                   >  7:130418/1‑149 (MQ=255)
              ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGc                                                                                                                                 <  3:27212/147‑1 (MQ=255)
                 cGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGc                                                                                                                            >  7:271776/1‑149 (MQ=255)
                   aaaaGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCa                                                                                                                          <  8:271776/149‑1 (MQ=255)
                     catcTGATTCAGGGTGACAAGTACATCTTCATTCGCTATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGa                                                                                                                        <  4:39706/146‑1 (MQ=255)
                      aGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGaa                                                                                                                       >  3:16029/1‑149 (MQ=255)
                       gCTGATTCAGGGTGACAAGTACATCTGTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGaa                                                                                                                       >  1:171249/1‑149 (MQ=255)
                       gCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAAt                                                                                                                      >  5:321438/1‑149 (MQ=255)
                          gATTCAGGGTGACAAGTACATCTAATTAGCAAAGCCTACCTGAAAAACCACAAATTAAAGGTAAAGGATGGCAAGGTAGAAAACGAATTATGAAATGTGAACATGGGGCCGCAAAGCGGCCTCGGGTTAAATGCATGGCCAGAATCCg                                                                                                                   >  7:204434/1‑148 (MQ=255)
                           aTTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGCCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATcc                                                                                                                    <  3:200445/147‑1 (MQ=255)
                                 ggTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGccc                                                                                                             >  8:221272/1‑148 (MQ=255)
                                    gACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAAc                                                                                                         >  3:227877/1‑149 (MQ=255)
                                       aaGTACATCTTCATTCGCAATGCCTACCTTCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTc                                                                                                      <  5:292918/149‑1 (MQ=255)
                                          tACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgc                                                                                                   <  2:171249/149‑1 (MQ=255)
                                            cATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgctg                                                                                                 >  3:183516/1‑149 (MQ=255)
                                               cTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCAt                                                                                                                               <  6:42774/116‑1 (MQ=255)
                                               cTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCAt                                                                                                                               >  5:42774/1‑116 (MQ=255)
                                               cTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgccc                                                                                              >  5:123724/1‑149 (MQ=255)
                                               cTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgccc                                                                                              >  6:243127/1‑149 (MQ=255)
                                                    aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgccc                                                                                          <  3:82884/149‑1 (MQ=255)
                                                     ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCCAGGTCGAAGACGCCTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCCGCCCGACCTTCTGCTGcccgcccg                                                                                         >  2:18526/1‑149 (MQ=255)
                                                       cGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCg                                                                                                                   <  3:173610/121‑1 (MQ=255)
                                                            tGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTc                                                                                  >  6:103254/1‑149 (MQ=255)
                                                             gCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCg                                                                                 >  8:169112/1‑149 (MQ=255)
                                                              ccTACCTGCAACACCGCGAGTTCAAGGTGAAGGATGGCGACGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCgg                                                                                <  4:73582/149‑1 (MQ=255)
                                                              ccTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCgg                                                                                <  4:110818/149‑1 (MQ=255)
                                                              ccTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCgg                                                                                <  7:221272/149‑1 (MQ=255)
                                                                tACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATATGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGc                                                                               >  7:55089/1‑148 (MQ=255)
                                                                 aCCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgc                                                                                                   <  7:37186/127‑1 (MQ=255)
                                                                 aCCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgc                                                                                                   >  8:37186/1‑127 (MQ=255)
                                                                 aCCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCaa                                                                             >  6:232517/1‑149 (MQ=255)
                                                                    tGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATc                                                                           >  4:31781/1‑148 (MQ=255)
                                                                       aaaaCCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATa                                                                       <  6:227175/149‑1 (MQ=255)
                                                                       aaaaCCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATa                                                                       >  8:155392/1‑149 (MQ=255)
                                                                        aaaCCGCGAGTTCAAGGTGAAGGATGGTGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATAc                                                                      <  6:1571/149‑1 (MQ=255)
                                                                           ccGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGCCGGCAATCCATACCAc                                                                   >  4:67277/1‑149 (MQ=255)
                                                                            cgcgAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCgg                                                                                <  4:174684/135‑1 (MQ=255)
                                                                            cgcgAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCgg                                                                                >  3:174684/1‑135 (MQ=255)
                                                                            cgcgAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACt                                                                  >  4:14340/1‑149 (MQ=255)
                                                                             gcgAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGCTCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCAGGTCGGCACTCCATACCACtt                                                                 >  2:87869/1‑149 (MQ=255)
                                                                                 gTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGt                                                             >  7:273612/1‑149 (MQ=255)
                                                                                  ttCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGt                                                             <  4:183516/148‑1 (MQ=255)
                                                                                  ttCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTg                                                            <  8:204434/149‑1 (MQ=255)
                                                                                       ggTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTgg                                                        >  2:92899/1‑148 (MQ=255)
                                                                                           aaGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCCCTTCGGTGCTGGCATg                                                    >  5:177408/1‑148 (MQ=255)
                                                                                             ggATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAg                                                 >  5:231307/1‑149 (MQ=255)
                                                                                             ggATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAg                                                 >  8:152180/1‑149 (MQ=255)
                                                                                             ggATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCGGGTCGGCAATCCATCCCACTTCGGTGCTGGCATGCAg                                                 >  4:46473/1‑149 (MQ=255)
                                                                                             ggATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGATTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATc                                                                           <  5:207787/123‑1 (MQ=255)
                                                                                             ggATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGATTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATc                                                                           >  6:207787/1‑123 (MQ=255)
                                                                                                tGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGccc                                              >  4:122936/1‑149 (MQ=255)
                                                                                                    gAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGAATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCtt                                           >  3:266617/1‑148 (MQ=255)
                                                                                                          gAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGcc                                    >  5:150893/1‑149 (MQ=255)
                                                                                                          gAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTCCTGGCATGCAGCCCCTTCAGCGcc                                    >  6:52591/1‑149 (MQ=255)
                                                                                                           aaGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                                   <  7:169112/149‑1 (MQ=255)
                                                                                                            aGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCgg                                  <  3:234286/149‑1 (MQ=255)
                                                                                                              acgacTTCTGATCTGTGAACCTGGGGGCGCAAAGCGTCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAAGTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTTCTGGCATGCAGCCCCTTCAGCGCCTGAt                                <  8:55089/149‑1 (MQ=255)
                                                                                                              acgacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                                   <  5:112125/146‑1 (MQ=255)
                                                                                                              acgacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                                   >  6:112125/1‑146 (MQ=255)
                                                                                                              acgacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGAt                                <  5:243127/149‑1 (MQ=255)
                                                                                                               cgacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATg                               <  2:164596/149‑1 (MQ=255)
                                                                                                                gacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATg                               >  2:14368/1‑148 (MQ=255)
                                                                                                                  cTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGAc                            >  1:136440/1‑149 (MQ=255)
                                                                                                                  cTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGAc                            <  1:14368/149‑1 (MQ=255)
                                                                                                                  cTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGAc                            <  6:11889/149‑1 (MQ=255)
                                                                                                                      tGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCTGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAg                        >  4:128790/1‑149 (MQ=255)
                                                                                                                       gATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGa                             >  7:74411/1‑143 (MQ=255)
                                                                                                                       atgcattgAACCCGAGGCCGCTTTGCGGCCCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGt                       >  8:152408/7‑149 (MQ=38)
                                                                                                                        aTCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGt                       <  4:227877/148‑1 (MQ=255)
                                                                                                                           attgAACCCGAGGCCGCTTTGCGGCCCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt                   >  1:34626/3‑149 (MQ=38)
                                                                                                                             tgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                           <  7:152180/139‑1 (MQ=255)
                                                                                                                             tgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                           <  4:16029/139‑1 (MQ=255)
                                                                                                                             tgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATcc                 <  6:127076/149‑1 (MQ=255)
                                                                                                                             tgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATcc                 >  5:127076/1‑149 (MQ=255)
                                                                                                                             tgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATc                  <  7:155392/148‑1 (MQ=255)
                                                                                                                             tgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATc                  <  1:92899/148‑1 (MQ=255)
                                                                                                                               aaCCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATCGACTCCGTCCATCCgg               <  2:9980/149‑2 (MQ=255)
                                                                                                                                  cTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAcc            >  1:189987/1‑149 (MQ=255)
                                                                                                                                  cTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAcc            >  6:67485/1‑149 (MQ=255)
                                                                                                                                   tGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGCTGGACg                           <  6:164667/133‑1 (MQ=255)
                                                                                                                                   tGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                           >  5:164667/1‑133 (MQ=255)
                                                                                                                                    ggggCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTc          >  5:131309/1‑149 (MQ=255)
                                                                                                                                     gggCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTc          <  5:52591/148‑1 (MQ=255)
                                                                                                                                      ggCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGa        <  5:67485/149‑1 (MQ=255)
                                                                                                                                       gCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGAt       <  6:155143/149‑1 (MQ=255)
                                                                                                                                         cGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCt     >  6:305030/1‑149 (MQ=255)
                                                                                                                                          gCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCt     >  1:269900/1‑148 (MQ=255)
                                                                                                                                           cAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGa   >  5:54531/1‑149 (MQ=255)
                                                                                                                                             aaGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGAt  >  1:182941/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGAT  >  NC_002947/2469712‑2469999

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: