Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,469,862 G→A 14.0% 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,8620GA14.0% 133.8 / 12.4 57intergenic (+31/+4)vacJ/PP_2164lipoprotein VacJ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/22);  new base A (3/5);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.96e-01

CAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATCACCCACCAG  >  NC_002947/2469714‑2470009
                                                                                                                                                    |                                                                                                                                                   
cAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCg                                                                                                                                                     >  2:208707/1‑149 (MQ=255)
cAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACAGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGCCTTCTGCTCTGTGAACCTGGGGCCGCAAAGCGGCCTCg                                                                                                                                                     >  3:94057/1‑149 (MQ=255)
  gCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCggg                                                                                                                                                   <  5:277731/149‑1 (MQ=255)
        tctcGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAAt                                                                                                                                             >  7:51575/1‑149 (MQ=255)
        tctcGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTATGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAAt                                                                                                                                             >  7:56188/1‑149 (MQ=255)
          tcGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGCAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGTTCaca                                                                                                                                             <  7:169672/147‑3 (MQ=38)
          tcGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATg                                                                                                                                            >  5:214906/1‑148 (MQ=255)
          tcGGCCGAAAAGCTGATTCAGGGGGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGTTCaca                                                                                                                                             <  5:292742/147‑3 (MQ=38)
            ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCAt                                                                                                                                         >  4:803/1‑149 (MQ=255)
            ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGTTCacaga                                                                                                                                           <  5:330181/147‑5 (MQ=38)
            ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGTTCac                                                                                                                                              >  3:141296/1‑143 (MQ=38)
            ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGCCCCAGGTTCac                                                                                                                                              <  4:141296/144‑2 (MQ=38)
             gCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCTTTGCGGCCCCAGGTTCacagatca                                                                                                                                        <  4:209175/149‑8 (MQ=255)
              ccGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATgg                                                                                                                                       <  8:82205/149‑1 (MQ=255)
              ccGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATgg                                                                                                                                       <  6:239243/149‑1 (MQ=255)
                 aaaaGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGc                                                                                                                                      >  1:219137/1‑147 (MQ=255)
                    aGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGaa                                                                                                                                 <  2:181716/149‑1 (MQ=255)
                     gCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGAAGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAAt                                                                                                                                >  3:56322/1‑149 (MQ=255)
                         aTTCAGGGTTACAATTACATCTTCATTCTCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCTGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCg                                                                                                                             <  8:56188/148‑1 (MQ=255)
                         aTTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCg                                                                                                                             <  8:51575/148‑1 (MQ=255)
                         aTTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGc                                                                                                                            <  1:314543/149‑1 (MQ=255)
                                   aCAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAAc                                                                                                                   <  2:219137/148‑1 (MQ=255)
                                          cATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGAGCCGCAAACCGGCCTCGGGTTc                                                                                                                                                >  8:351177/1‑112 (MQ=255)
                                              ttCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccg                                                                                                       >  1:199708/1‑149 (MQ=255)
                                                  ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCgg                                                                                                                                                    <  5:235677/100‑1 (MQ=255)
                                                  ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCgg                                                                                                                                                    >  6:235677/1‑100 (MQ=255)
                                                  ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGcccgcccg                                                                                                   >  7:117628/1‑149 (MQ=255)
                                                           ccTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCg                                                                                           <  8:292883/148‑1 (MQ=255)
                                                                  gCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGAGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATcc                                                                                    >  8:116830/1‑148 (MQ=255)
                                                                   cAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCAt                                                                                  <  2:199708/149‑1 (MQ=255)
                                                                   cAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGCGCCGCAAAGCGGCCTCGGGTCCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCa                                                                                   <  7:242107/148‑1 (MQ=255)
                                                                          gcgAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACtt                                                                           >  8:98547/1‑149 (MQ=255)
                                                                          gcgAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACtt                                                                           >  5:125572/1‑149 (MQ=255)
                                                                          gcgAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACtt                                                                           >  4:179580/1‑149 (MQ=255)
                                                                                  aaGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCt                                                                    <  1:208707/148‑1 (MQ=255)
                                                                                      tGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCAt                                                               <  6:125572/149‑1 (MQ=255)
                                                                                            aTGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGcc                                                         >  5:383698/1‑149 (MQ=255)
                                                                                               gCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCGCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGcccc                                                       <  4:88415/148‑1 (MQ=255)
                                                                                                     tCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAgcg                                                >  3:2780/1‑149 (MQ=255)
                                                                                                       gAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGcc                                              <  3:179580/149‑1 (MQ=255)
                                                                                                             gacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATg                                         >  7:271956/1‑148 (MQ=255)
                                                                                                               cTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGa                                       <  8:117628/148‑1 (MQ=255)
                                                                                                                ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTg                                                                   >  4:204596/1‑119 (MQ=255)
                                                                                                                ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTg                                                                   <  3:204596/119‑1 (MQ=255)
                                                                                                                ttCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCAGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg                                     >  5:359644/1‑149 (MQ=255)
                                                                                                                   tGATCTGTGAACCTGGGGCCGCGAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAg                                  >  8:18232/1‑149 (MQ=255)
                                                                                                                        tgtgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt                             >  1:150380/1‑149 (MQ=255)
                                                                                                                        attgAACCCGAGGCCGCTTTGCGGCCCCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt                             >  7:190095/3‑149 (MQ=255)
                                                                                                                        attgAACCCGAGGCCGCTTTGCGGCCCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt                             >  8:219617/3‑149 (MQ=38)
                                                                                                                        attgAACCCGAGGCCGCTTTGCGGCCCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt                             >  3:313339/3‑149 (MQ=38)
                                                                                                                         gtgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATc                            >  1:116064/1‑149 (MQ=255)
                                                                                                                          cgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATcc                           <  7:116830/148‑1 (MQ=255)
                                                                                                                            aaCCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCg                          >  5:267133/1‑148 (MQ=255)
                                                                                                                              ccTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAc                       >  1:63658/1‑149 (MQ=255)
                                                                                                                               cTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTAcc                      >  1:224405/1‑149 (MQ=255)
                                                                                                                                  gggCCGCAACCCGGCCGCGGGCTCAATGCATGGCCAGAATCCGCAGCCCGCACCTCTGCTGCCCGCCCGGCGGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCt                     <  6:214906/147‑1 (MQ=255)
                                                                                                                                     ccGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATc                <  2:63658/149‑1 (MQ=255)
                                                                                                                                       gCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCt               >  6:55520/1‑148 (MQ=255)
                                                                                                                                        cAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGa             <  4:2780/149‑1 (MQ=255)
                                                                                                                                                ccTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATcacccac     <  4:94057/149‑1 (MQ=255)
                                                                                                                                                ccTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCGGATcacccac     >  7:149869/1‑149 (MQ=255)
                                                                                                                                                 cTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATcacccacc    <  3:803/149‑1 (MQ=255)
                                                                                                                                                    gggTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATCACCCACCAg  >  1:46308/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCGTACCTCGATCTGATCACCCACCAG  >  NC_002947/2469714‑2470009

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