Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,469,844 G→A 7.1% 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,8440GA 7.1% 243.8 / 6.5 84intergenic (+13/+22)vacJ/PP_2164lipoprotein VacJ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/37);  new base A (4/2);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAC  >  NC_002947/2469697‑2469988
                                                                                                                                                    |                                                                                                                                                 
cTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTCACAAGTACATCTTCATTCGCAATGCCTAACTGCAAAACAGCGAGTTCAAGGTGAAGGATGTCGAGGTCGAAGACGACATCTGATCTGTGAACCTggg                                                                                                                                                  >  6:268589/1‑149 (MQ=255)
       acacCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCa                                                                                                                                            >  7:231579/1‑148 (MQ=255)
        cacCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCaaa                                                                                                                                          <  3:50055/149‑1 (MQ=255)
         acCCGCGCCAGCCTGCACTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCAAGTTCAAGGTGAAGGATGGCGAAATCGAAGACGACTTCTGATCTCTGAACCAGGCGCCGCAAAg                                                                                                                                         >  1:145911/1‑149 (MQ=255)
          cccGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGc                                                                                                                                        <  5:171247/149‑1 (MQ=255)
           ccGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGctttgcg                                                                                                                                       >  2:213901/1‑143 (MQ=255)
                    ccTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGt                                                                                                                             <  3:171982/149‑1 (MQ=255)
                     cTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGtt                                                                                                                            <  8:231579/149‑1 (MQ=255)
                     cTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGtt                                                                                                                            <  1:213901/149‑1 (MQ=255)
                       gCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTc                                                                                                                           >  8:247615/1‑148 (MQ=255)
                           tcGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAAt                                                                                                                        <  1:438159/147‑1 (MQ=255)
                             ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGc                                                                                                                      <  1:45078/147‑1 (MQ=255)
                             ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGc                                                                                                                      <  7:56012/147‑1 (MQ=255)
                             ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCAt                                                                                                                    <  3:388456/149‑1 (MQ=255)
                             ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCAt                                                                                                                    >  3:335892/1‑149 (MQ=255)
                                   aaaTCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCa                                                                                                               <  8:152896/148‑1 (MQ=255)
                                         gATTCAGGGTGACAAGTACGTCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGAATGGCCAGAATCCg                                                                                                        >  1:307990/1‑149 (MQ=255)
                                         gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGtt                                                                                                                            >  6:168383/1‑129 (MQ=38)
                                         gATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGtt                                                                                                                            <  5:168383/129‑1 (MQ=38)
                                           ttCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCg                                                                                                        >  6:273558/1‑147 (MQ=255)
                                              aGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGcc                                                                                                   >  8:324763/1‑149 (MQ=255)
                                               gggTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGccc                                                                                                  >  3:7259/1‑149 (MQ=255)
                                               gggTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGcca                                                                                                  <  2:102205/149‑2 (MQ=255)
                                                ggTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGccc                                                                                                  >  1:205110/1‑148 (MQ=255)
                                                      aaGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTc                                                                                           <  1:399997/149‑1 (MQ=255)
                                                        gTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGc                                                                                                                                     >  5:275327/1‑106 (MQ=255)
                                                        gTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGc                                                                                                                                     <  6:275327/106‑1 (MQ=255)
                                                         tACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgc                                                                                        >  8:398959/1‑149 (MQ=255)
                                                           cATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctgct                                                                                       <  2:240962/148‑1 (MQ=255)
                                                               ttCATTAGCAATGCCGACCAGCAAAACCCCGCGGTGAAGGTGAAGGATGGCGAGGTCGAAGACGACGGGTGATCTGTGAACCAGGGGGCGCAAAGCGGCCGCCGGGTGGAAGGCTGGCCAGAAAACGCCCCCCGACCTTCAGGTgcccg                                                                                  >  8:356840/1‑149 (MQ=255)
                                                                  aTTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgccc                                                                               >  2:46408/1‑149 (MQ=255)
                                                                   ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCg                                                                                                 >  7:304930/1‑130 (MQ=255)
                                                                   ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCg                                                                                                 <  8:304930/130‑1 (MQ=255)
                                                                   ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGcccgcccg                                                                              >  8:202502/1‑149 (MQ=255)
                                                                   ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGcccgcccg                                                                              >  6:144693/1‑149 (MQ=255)
                                                                     cGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCgg                                                                             <  4:335892/148‑1 (MQ=255)
                                                                       cAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGAACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCt                                                                           >  1:330947/1‑149 (MQ=255)
                                                                            ccTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCgg                                                                     <  4:7259/149‑1 (MQ=255)
                                                                             cTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTc                                                                                                                                 >  2:235644/1‑88 (MQ=255)
                                                                             cTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTc                                                                                                                                 <  1:235644/88‑1 (MQ=255)
                                                                             cTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGGCgg                                                                     >  6:232650/1‑148 (MQ=255)
                                                                              tACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGc                                                                    >  4:317798/1‑148 (MQ=255)
                                                                                 cTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATc                                                                <  7:247615/149‑1 (MQ=255)
                                                                                   gCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCa                                                              >  5:19975/1‑149 (MQ=255)
                                                                                    cAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCa                                                              >  3:311108/1‑148 (MQ=255)
                                                                                         ccGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCAc                                                        >  1:162062/1‑149 (MQ=255)
                                                                                         ccGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCAc                                                        >  1:271248/1‑149 (MQ=255)
                                                                                           gcgAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACtt                                                      >  3:343289/1‑149 (MQ=255)
                                                                                             gAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCg                                                    <  7:324763/149‑1 (MQ=255)
                                                                                             gAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCg                                                    >  4:96587/1‑149 (MQ=255)
                                                                                                   aaGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTg                                              <  2:330947/149‑1 (MQ=255)
                                                                                                   aaGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTg                                              <  3:317798/149‑1 (MQ=255)
                                                                                                     ggTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGcc                                                                                                                <  2:252827/81‑1 (MQ=255)
                                                                                                     ggTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGcc                                                                                                                >  1:252827/1‑81 (MQ=255)
                                                                                                     ggTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGc                                            <  5:416332/149‑1 (MQ=255)
                                                                                                     ggTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGc                                            <  1:223414/149‑1 (MQ=255)
                                                                                                      gTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCa                                           >  6:336948/1‑149 (MQ=255)
                                                                                                       tGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCAt                                          <  5:273558/149‑1 (MQ=255)
                                                                                                        gAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATg                                         >  2:70482/1‑149 (MQ=255)
                                                                                                         aaGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGTGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGGCGGCAATCGATACCACTTCGGTGCTGGCATGc                                        <  2:307990/149‑1 (MQ=255)
                                                                                                         aaGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATg                                         <  4:343289/148‑1 (MQ=255)
                                                                                                           ggATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAg                                      >  8:299607/1‑149 (MQ=255)
                                                                                                                   aGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCAGAACTTCTGCTGCCCGCCCGGCTGGTCGGCATTCGATACCTCTTGTGTGCTGGCATGCAGCCCCTTCa                              <  2:145911/149‑1 (MQ=255)
                                                                                                                      tCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAgcg                           >  6:18848/1‑149 (MQ=255)
                                                                                                                       cGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGccc                                   <  8:320455/140‑1 (MQ=255)
                                                                                                                       cGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCAAGCAGccc                                   >  7:320455/1‑140 (MQ=255)
                                                                                                                        gAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGcc                         >  5:7157/1‑149 (MQ=255)
                                                                                                                             acgacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGAt                     >  8:139460/1‑149 (MQ=255)
                                                                                                                              cgacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATg                    <  7:139460/149‑1 (MQ=255)
                                                                                                                               gacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATg                    >  8:161176/1‑148 (MQ=255)
                                                                                                                               gacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATg                    >  4:165614/1‑148 (MQ=255)
                                                                                                                                 cTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGAc                 <  5:58620/149‑1 (MQ=255)
                                                                                                                                  ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                <  5:144693/149‑1 (MQ=255)
                                                                                                                                  ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                <  5:18848/149‑1 (MQ=255)
                                                                                                                                    cTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                >  1:225895/1‑147 (MQ=255)
                                                                                                                                       aTCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGt            <  7:299607/148‑1 (MQ=255)
                                                                                                                                          tgtgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt        >  5:439290/1‑149 (MQ=255)
                                                                                                                                          attgAACCCGAGGCCGCTTTGCGGCCCCAGGTTCAATGCGTGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt        >  4:37510/3‑149 (MQ=38)
                                                                                                                                            tgAACCTGGGGCCTCTAAGCGGCCTCGGGTTGAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt        <  5:268589/147‑1 (MQ=255)
                                                                                                                                             gAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCg     >  1:98113/1‑149 (MQ=255)
                                                                                                                                              aaCCTGGGGCCGCAACGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCCGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGt    <  6:253191/149‑1 (MQ=255)
                                                                                                                                              aaCCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATCGACTCCGTCCATCCg     <  7:202502/148‑1 (MQ=255)
                                                                                                                                              aaCCCGAGGCCGCTTTGCGGCCCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCg     >  4:300464/1‑148 (MQ=255)
                                                                                                                                                ccTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAc  >  2:180534/1‑149 (MQ=255)
                                                                                                                                                ccTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAc  >  7:202078/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAC  >  NC_002947/2469697‑2469988

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: