Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,265,951 C→G 5.6% D170H (GAC→CAC)  PP_3738 ← transcriptional regulator VanR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,265,9510CG 5.6% 226.3 / 3.8 71D170H (GAC→CAC) PP_3738transcriptional regulator VanR
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base G (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

TCGGCATAGCGCAAGGTGGCATTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAGGTCGTGATAGACCTGCAAGTCGTCTTCGCTC  >  NC_002947/4265803‑4266096
                                                                                                                                                    |                                                                                                                                                   
tCGGCATAGCGCAAGGTGGCATTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGACAGATCATGGCGGTc                                                                                                                                                     >  8:93208/1‑149 (MQ=255)
   gCATAGCGCAAGGTGGCATTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCAc                                                                                                                                                  >  3:295170/1‑149 (MQ=255)
     aTAGCGCAAGGTGGCATTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCg                                                                                                                                                >  7:183690/1‑149 (MQ=255)
         cgcAAGGTGGCATTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAg                                                                                                                                            >  6:316172/1‑149 (MQ=255)
          gcAAGGTGGCATTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCAGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGt                                                                                                                                           >  5:329648/1‑149 (MQ=255)
          gcAAGGTGGCATTGGCGTGCTCGCGCATGATCGCTGCGGCGCGGGCACCCTGGCGGTTGACCCGGGCGTCGAACACGGCGGGGGGCTGCATATGGGCGCAGTTGACGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCCCCGCCAGt                                                                                                                                           >  5:329653/1‑149 (MQ=255)
           cAAGGTGGCATTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCTGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTg                                                                                                                                          <  4:156869/149‑1 (MQ=255)
            aaGGTGGCATTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTGGACCGCCatgat                                                                                                                                         <  8:174152/149‑5 (MQ=255)
               gTGGCATTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGt                                                                                                                                      >  4:348919/1‑149 (MQ=255)
                    aTTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGAcc                                                                                                                                  >  6:150350/1‑148 (MQ=255)
                       ggCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAg                                                                                                                               <  5:274545/148‑1 (MQ=255)
                        gCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAgg                                                                                                                              <  6:216885/148‑1 (MQ=255)
                        gCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGc                                                                                                                             <  5:131160/149‑1 (MQ=255)
                        gCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGc                                                                                                                             <  5:316172/149‑1 (MQ=255)
                         cGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGc                                                                                                                             >  3:368028/1‑148 (MQ=255)
                          gTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑Ga                                                                                                                          <  4:295170/149‑1 (MQ=255)
                            gCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAAt                                                                                                                        <  4:25482/149‑1 (MQ=255)
                               cgcgCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGGGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGc                                                                                                                     <  5:374562/149‑1 (MQ=255)
                                gcgcATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATgg                                                                                                                      >  2:221801/1‑147 (MQ=255)
                                 cgcATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAg                                                                                                                   <  2:295902/149‑1 (MQ=255)
                                     tGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTGGACCGCCa                                                                                                                                             >  8:131214/1‑120 (MQ=255)
                                     tGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTGGACCGCCa                                                                                                                                             <  7:131214/120‑1 (MQ=255)
                                      gATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGAt                                                                                                              >  2:172589/1‑149 (MQ=255)
                                          gCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGtt                                                                                                          <  7:93208/149‑1 (MQ=255)
                                          gCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGt                                                                                                           <  3:185850/148‑1 (MQ=255)
                                           cTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTg                                                                                                         <  1:172589/149‑1 (MQ=255)
                                            ttCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTgc                                                                                                        >  4:357730/1‑149 (MQ=255)
                                            ttCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTgc                                                                                                        >  6:2140/1‑149 (MQ=255)
                                               ggCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTgcgc                                                                                                      <  6:299826/148‑1 (MQ=255)
                                                       cACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCATGTGATCGTTGCGCG‑CCAGGGc                                                                                             >  7:55902/1‑148 (MQ=255)
                                                         cccTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAAc                                                                                          >  3:41837/1‑149 (MQ=255)
                                                          ccTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCaatg                                                                                         >  7:267092/1‑147 (MQ=255)
                                                          ccTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACg                                                                                         >  2:83468/1‑149 (MQ=255)
                                                           cTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACgg                                                                                        <  1:221801/149‑1 (MQ=255)
                                                            tGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACgg                                                                                        >  5:230381/1‑148 (MQ=255)
                                                             ggCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGc                                                                                       >  3:176873/1‑148 (MQ=255)
                                                               cGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAAt                                                                                    <  4:176873/149‑1 (MQ=255)
                                                                  ttGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGc                                                                                 <  5:150350/149‑1 (MQ=255)
                                                                     aCCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATCCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCgg                                                                              <  5:272508/149‑1 (MQ=255)
                                                                      ccAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCggg                                                                             <  3:357730/149‑1 (MQ=255)
                                                                           gCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTtgtg                                                                        >  3:300119/1‑149 (MQ=255)
                                                                                      ggCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCga                                                             <  4:368028/149‑1 (MQ=255)
                                                                                            gtgCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCa                                                       >  2:311181/1‑149 (MQ=255)
                                                                                               ctgcATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAAt                                                                                    <  2:366589/117‑1 (MQ=255)
                                                                                                 gcATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCGCCCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCtg                                                   <  1:83468/149‑1 (MQ=255)
                                                                                                 gcATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGCTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCAGCtga                                                  >  5:201818/1‑149 (MQ=255)
                                                                                                  cATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGcc                                                                                <  6:281370/118‑1 (MQ=255)
                                                                                                  cATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGcc                                                                                >  5:281370/1‑118 (MQ=255)
                                                                                                  cATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCtgat                                                 >  4:103717/1‑149 (MQ=255)
                                                                                                  cATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCtga                                                  <  8:183690/148‑1 (MQ=255)
                                                                                                  cATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCAAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCtgat                                                 >  5:345125/1‑149 (MQ=255)
                                                                                                     atGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCgg                                                                                                                                       <  4:306750/62‑1 (MQ=255)
                                                                                                     atGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCgg                                                                                                                                       >  3:306750/1‑62 (MQ=255)
                                                                                                         gCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAgcg                                          <  3:103717/149‑1 (MQ=255)
                                                                                                          cGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAgcgc                                         <  3:348919/149‑1 (MQ=255)
                                                                                                            aaGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCAt                                       >  8:151653/1‑149 (MQ=255)
                                                                                                              gTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCg                                                                 >  5:125682/1‑121 (MQ=255)
                                                                                                                tGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCATTGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGtt                                    >  3:67624/1‑148 (MQ=255)
                                                                                                                tGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCATTGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGtt                                    >  3:72598/1‑148 (MQ=255)
                                                                                                                 gAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTc                                   <  7:151653/148‑1 (MQ=255)
                                                                                                                  aaGCGCCGGTTTTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGGGATCGTGGCGCG‑CCAGGGCAACGGAAATGGCCGGGTTGGGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAg                                 <  6:329653/149‑1 (MQ=255)
                                                                                                                  aaGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCGAGTGCGGTGACCGAGGC‑GAATGGCATGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAg                                 <  6:329648/149‑1 (MQ=255)
                                                                                                                    gcgcCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAGGt                               >  6:188036/1‑149 (MQ=255)
                                                                                                                       ccGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAGGTCGt                            >  8:329546/1‑149 (MQ=255)
                                                                                                                       ccGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAGGTCGt                            <  6:230381/149‑1 (MQ=255)
                                                                                                                                gcgAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAATGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAGGTCGTGATAGACCt                   <  8:267092/149‑1 (MQ=255)
                                                                                                                                gcgAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAGGTCGTGATAGACCt                   >  1:333540/1‑149 (MQ=255)
                                                                                                                                gcgAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGAGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCACGTCGAGATAGACCt                   >  4:178396/1‑149 (MQ=255)
                                                                                                                                        gATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGCAGAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAGGTCGTGATAGACCTGCAAgtc            >  1:67555/1‑149 (MQ=255)
                                                                                                                                           actGGCGGTGGACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGAGGATCACCTGATGAAAGCGCATGTTCAGGTCGTGATAGACCTGCAAGTCGTCt        >  3:19530/3‑149 (MQ=255)
                                                                                                                                             tGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAGGTCGTGATAGACCTGCAAGTCGTCTTc      <  2:333540/149‑1 (MQ=255)
                                                                                                                                                 ggTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAGGTCGTGATAGACCTGCAAGTCGTCTTCGCTc  <  5:2140/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGGCATAGCGCAAGGTGGCATTGGCGTGCTCGCGCATGATCGCTTCGGCGCGGGCACCCTGGCGGTTGACCAGGGCGTCGAACACGGCGTGGTGCTGCATATGGGCGAAGTTGAAGCGCCGGTATTCGCGAGCCAGATCATGGCGGTCCACCGCCAGTGCGGTGACCGAGGC‑GAATGGCAGGTGATCGTTGCGCG‑CCAGGGCAACGGCAATGGCCGGGTTGTGGCTGGCGTCGATGATCACCTGATGAAAGCGCATGTTCAGGTCGTGATAGACCTGCAAGTCGTCTTCGCTC  >  NC_002947/4265803‑4266096

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