Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,469,862 G→A 10.9% intergenic (+31/+4) 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,8620GA10.9% 241.3 / 14.9 91intergenic (+31/+4)vacJ/PP_2164lipoprotein VacJ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/43);  new base A (4/6);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

GCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATCACCCAC  >  NC_002947/2469716‑2470006
                                                                                                                                                      |                                                                                                                                                  
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                                                                                                                                         <  7:92505/148‑6 (MQ=38)
          ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGC‑‑CCCAGGTTCacagat                                                                                                                                         <  4:60701/148‑6 (MQ=38)
              gAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGc                                                                                                                                     >  1:227113/1‑148 (MQ=255)
                      gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCAt                                                                                                                                        >  2:50645/1‑137 (MQ=255)
                      gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCAt                                                                                                                                        <  1:50645/137‑1 (MQ=255)
                      gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCg                                                                                                                            >  4:82351/1‑149 (MQ=255)
                       aTTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATcc                                                                                                                             <  2:285323/147‑1 (MQ=255)
                       aTTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATcc                                                                                                                             <  2:279730/147‑1 (MQ=255)
                           aGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGcc                                                                                                                       >  7:63784/1‑149 (MQ=255)
                             ggTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCg                                                                                                                     <  5:296505/149‑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                                                                                                         <  1:275743/143‑1 (MQ=255)
                                               aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCtgctgc                                                                                                         >  2:275743/1‑143 (MQ=255)
                                               aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgccc                                                                                                   >  4:276937/1‑149 (MQ=255)
                                               aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgccc                                                                                                   <  3:82351/149‑1 (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                                                                <  1:311536/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                                                                                                                                       <  7:10650/50‑1 (MQ=255)
                                                                                                                ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATg                                                                                                                                       >  8:10650/1‑50 (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)
                                                                                                                                    gCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGAt              <  3:130927/149‑1 (MQ=255)
                                                                                                                                      cGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCTGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCAGGATGGACTGAGTCGATCCGGAACTCGAtat            >  2:291673/1‑147 (MQ=255)
                                                                                                                                      cGCAAAGCGGC‑‑CTCGGGTTCAATGCAAGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCt            <  2:346946/149‑1 (MQ=255)
                                                                                                                                        cAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGa          >  7:151078/1‑149 (MQ=255)
                                                                                                                                          aaGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATc        <  3:80620/149‑1 (MQ=255)
                                                                                                                                              ggC‑‑CCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGTATGGACGTAGTCGATCCGTACCTCGATCTGATcaccc    <  7:221269/149‑1 (MQ=255)
                                                                                                                                                    tCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                                          >  1:390682/1‑107 (MQ=255)
                                                                                                                                                    tCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                                          <  2:390682/107‑1 (MQ=255)
                                                                                                                                                     cGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATcacccac  >  2:7755/1‑146 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATG‑‑GCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGC‑‑CTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGC‑‑TGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATCACCCAC  >  NC_002947/2469716‑2470006

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: ‑