Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,469,844 G→A 12.2% intergenic (+13/+22) 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,8440GA12.2% 236.5 / 17.5 90intergenic (+13/+22)vacJ/PP_2164lipoprotein VacJ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/41);  new base A (4/7);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CCTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTC  >  NC_002947/2469696‑2469991
                                                                                                                                                      |                                                                                                                                                       
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                                                                                            <  6:319281/147‑1 (MQ=255)
                                                             aTCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgct                                                                                            >  5:319281/1‑147 (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                                                                             >  5:241006/1‑149 (MQ=255)
                                                                        cAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑Tg                                                                             >  3:341531/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                                                 <  8:312011/149‑1 (MQ=255)
                                                                                                    aaGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTg                                                 <  6:241006/149‑1 (MQ=255)
                                                                                                    aaGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTg                                                 <  1:311536/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                                 <  7:223436/144‑1 (MQ=255)
                                                                                                                           gAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCa                                 >  8:223436/1‑144 (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                           >  1:79502/3‑133 (MQ=38)
                                                                                                                                            attgAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                           <  2:79502/131‑1 (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)
                                                                                                                                                     ggggCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCACCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTc  >  4:148741/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
CCTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTC  >  NC_002947/2469696‑2469991

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: