Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,469,842 T→C 10.9% intergenic (+11/+24) 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,8420TC10.9% 242.0 / 16.3 93intergenic (+11/+24)vacJ/PP_2164lipoprotein VacJ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/42);  new base C (3/7);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACC  >  NC_002947/2469694‑2469989
                                                                                                                                                      |                                                                                                                                                       
gACCTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCt                                                                                                                                                         >  5:370444/1‑149 (MQ=255)
  ccTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTgg                                                                                                                                                       <  6:370444/149‑1 (MQ=255)
    tGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTgggg                                                                                                                                                     >  7:138070/1‑149 (MQ=255)
     gATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATT‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGc                                                                                                                                                    >  2:294198/1‑149 (MQ=255)
     gATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTggg                                                                                                                                                      <  5:59697/147‑1 (MQ=255)
     gATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTggg                                                                                                                                                      >  6:59697/1‑147 (MQ=255)
     gATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGc                                                                                                                                                    >  5:35070/1‑149 (MQ=255)
      aTCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGcc                                                                                                                                                   >  1:221649/1‑149 (MQ=255)
            acCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAg                                                                                                                                             >  4:130927/1‑149 (MQ=255)
                     aGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGGTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGc‑‑cccag                                                                                                                                  <  5:301748/149‑5 (MQ=255)
                     aGCCTGCTCTCGGCAGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGc‑‑cccag                                                                                                                                  <  7:351981/149‑5 (MQ=255)
                      gCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCggg                                                                                                                                 <  2:221649/149‑1 (MQ=255)
                       ccTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGt                                                                                                                                <  6:35070/149‑1 (MQ=255)
                        cTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGtt                                                                                                                               >  6:187031/1‑149 (MQ=255)
                        cTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGtt                                                                                                                               <  3:112919/149‑1 (MQ=255)
                           ctctCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCac                                                                                                                            <  5:175709/149‑2 (MQ=255)
                              tcGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAAt                                                                                                                           <  5:264365/147‑1 (MQ=255)
                              tcGCCCGAAAAGCTGATTCAGGGTGACAAGTACATCCTCATTCGCAATGCCTACCTGCAAAACCTCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCTAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAAt                                                                                                                           <  3:272812/147‑1 (MQ=255)
                                ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGCCACGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CGCGGGTTCaa                                                                                                                            <  1:115962/146‑1 (MQ=255)
                                ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCacagat                                                                                                                        <  4:60701/148‑6 (MQ=38)
                                ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCacagat                                                                                                                        <  7:92505/148‑6 (MQ=38)
                                    gAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGc                                                                                                                    >  1:227113/1‑148 (MQ=255)
                                            gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCAt                                                                                                                       <  1:50645/137‑1 (MQ=255)
                                            gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCAt                                                                                                                       >  2:50645/1‑137 (MQ=255)
                                            gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCg                                                                                                           >  4:82351/1‑149 (MQ=255)
                                             aTTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATcc                                                                                                            <  2:279730/147‑1 (MQ=255)
                                             aTTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATcc                                                                                                            <  2:285323/147‑1 (MQ=255)
                                                 aGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGcc                                                                                                      >  7:63784/1‑149 (MQ=255)
                                                   ggTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCg                                                                                                    <  5:296505/149‑1 (MQ=255)
                                                      gACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCaaa                                                                                                                                              >  6:395587/1‑106 (MQ=255)
                                                      gACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCaaa                                                                                                                                              <  5:395587/106‑1 (MQ=255)
                                                         aaGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTc                                                                                              >  2:325529/1‑149 (MQ=255)
                                                             aCAACTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAACGTGAACGATG‑‑GCGACGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgct                                                                                          >  3:326756/1‑149 (MQ=255)
                                                              cATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgct                                                                                          <  6:143287/148‑1 (MQ=255)
                                                               atcatcCTTCGCAATGCCTCCCTGCAAACCCGCGAGTTCAAGGGGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCtgctgc                                                                                        <  4:212939/149‑1 (MQ=255)
                                                               aTCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgct                                                                                          >  5:319281/1‑147 (MQ=255)
                                                               aTCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgct                                                                                          <  6:319281/147‑1 (MQ=255)
                                                                tCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcc                                                                                       >  4:260736/1‑149 (MQ=255)
                                                                     aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCtgctgc                                                                                        >  2:275743/1‑143 (MQ=255)
                                                                     aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCtgctgc                                                                                        <  1:275743/143‑1 (MQ=255)
                                                                     aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgccc                                                                                  <  3:82351/149‑1 (MQ=255)
                                                                     aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgccc                                                                                  >  4:276937/1‑149 (MQ=255)
                                                                      ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGcccgcccg                                                                                 >  5:130776/1‑149 (MQ=255)
                                                                        cGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCgg                                                                                <  2:365550/148‑1 (MQ=255)
                                                                         gCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATT‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCATTGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑t                                                                            <  1:294198/149‑1 (MQ=255)
                                                                         gCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑t                                                                            >  1:393510/1‑149 (MQ=255)
                                                                          cAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑Tg                                                                           >  3:341531/1‑149 (MQ=255)
                                                                          cAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑Tg                                                                           >  5:241006/1‑149 (MQ=255)
                                                                            aTGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGt                                                                         <  4:136010/149‑1 (MQ=255)
                                                                                cTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCgg                                                                      >  8:43160/1‑148 (MQ=255)
                                                                                   ccTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCa                                                                    <  3:276937/147‑1 (MQ=255)
                                                                                    cTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATc                                                                 <  4:331411/149‑1 (MQ=255)
                                                                                      gCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCa                                                               >  2:236003/1‑149 (MQ=255)
                                                                                       cAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCAt                                                              <  2:393510/149‑1 (MQ=255)
                                                                                         aaaCCGCGAGTTCTTGGTGTAGGAGG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACGTGGGGGGGCATTGCGGC‑‑CGCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGGCCGCCCGGC‑‑TGGTCGGCAATCCATAc                                                            <  4:326756/149‑1 (MQ=255)
                                                                                         aaaCCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATAc                                                            <  8:63784/149‑1 (MQ=255)
                                                                                         aaaCCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGACAGAATCCGCAGCCCGAACTTCTGCTGCCCGCACGGC‑‑TGGTCGGCAATCCATAc                                                            >  7:221239/1‑149 (MQ=255)
                                                                                                 aGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCgg                                                    <  5:187031/149‑1 (MQ=255)
                                                                                                   ttCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTg                                                  >  4:312133/1‑149 (MQ=255)
                                                                                                      aaGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCGGCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGc                                                 >  1:120726/1‑149 (MQ=255)
                                                                                                      aaGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTg                                               <  1:311536/149‑1 (MQ=255)
                                                                                                      aaGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTg                                               <  8:312011/149‑1 (MQ=255)
                                                                                                      aaGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTg                                               <  6:241006/149‑1 (MQ=255)
                                                                                                          tGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCAt                                           >  4:223068/1‑149 (MQ=255)
                                                                                                          tGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCAt                                           <  3:260736/149‑1 (MQ=255)
                                                                                                            aaGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATg                                          <  1:394289/148‑1 (MQ=255)
                                                                                                              ggATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAg                                       >  5:281427/1‑149 (MQ=255)
                                                                                                                 tG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGccc                                    <  1:236003/149‑1 (MQ=255)
                                                                                                                          gTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAgc                             >  6:87414/1‑149 (MQ=255)
                                                                                                                           tCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAgcg                            >  3:13500/1‑149 (MQ=255)
                                                                                                                           tCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAgcg                            >  8:51534/1‑149 (MQ=255)
                                                                                                                             gAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCa                               >  8:223436/1‑144 (MQ=255)
                                                                                                                             gAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCa                               <  7:223436/144‑1 (MQ=255)
                                                                                                                                     cTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGAc                  >  6:277423/1‑149 (MQ=255)
                                                                                                                                     cTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGAc                  >  2:297177/1‑149 (MQ=255)
                                                                                                                                      ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATg                                                                                                                      >  8:10650/1‑50 (MQ=255)
                                                                                                                                      ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATg                                                                                                                      <  7:10650/50‑1 (MQ=255)
                                                                                                                                      ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                 >  8:54089/1‑149 (MQ=255)
                                                                                                                                      ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                 >  7:992/1‑149 (MQ=255)
                                                                                                                                       tCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGcccc                                   <  2:306178/130‑1 (MQ=255)
                                                                                                                                       tCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGcccc                                   >  1:306178/1‑130 (MQ=255)
                                                                                                                                        cTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                 >  4:80620/1‑147 (MQ=255)
                                                                                                                                         catgcattgAACCCGAGGCCGCTTTGCGGA‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAg              >  7:371498/8‑149 (MQ=38)
                                                                                                                                           aTCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGa               <  6:130776/146‑1 (MQ=255)
                                                                                                                                            tctggtcaCCTGGGGGCCGAAACCGGG‑‑CGCGGGGTCAATGGATTGCCATTAACCGCATCCCGAACTTCTGCCGCCCGCCCGGCCTTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAg              <  1:325529/141‑1 (MQ=255)
                                                                                                                                              attgAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                         <  2:79502/131‑1 (MQ=38)
                                                                                                                                              attgAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                         >  1:79502/3‑133 (MQ=38)
                                                                                                                                              attgAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt         >  3:24633/3‑149 (MQ=38)
                                                                                                                                              attgAACCAGAGGCCGCTTTGCGGA‑‑CCCACGTTCAATCCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCCGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTCGCATGCAGCCCCTTCAGCGCAGGATTGACGGAGTCGAt         >  8:221269/3‑149 (MQ=255)
                                                                                                                                                tgAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTc            <  4:13500/144‑1 (MQ=255)
                                                                                                                                                tgAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATcc       <  7:43160/149‑1 (MQ=255)
                                                                                                                                                 gAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATcc       <  4:24633/148‑1 (MQ=255)
                                                                                                                                                   acatGGGGCCGGAATGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCTCAGCGCGAACTTCTGCTGCACGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGt     <  8:221239/145‑1 (MQ=255)
                                                                                                                                                    ccTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCGATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAc   <  8:138070/149‑1 (MQ=255)
                                                                                                                                                     cTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAcc  <  5:277423/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
GACCTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACC  >  NC_002947/2469694‑2469989

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