Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,487,011 A→G 6.5% K21E (AAG→GAG)  PP_2181 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,487,0110AG 6.5% 174.6 / 2.6 61K21E (AAG→GAG) PP_2181Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base G (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGCGAAGAGCTCGGTGTATCCACAGCCTACTTCTACAA  >  NC_002947/2486868‑2487159
                                                                                                                                               |                                                                                                                                                    
cTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCaa                                                                                                                                                     >  8:92684/1‑145 (MQ=255)
cTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCaa                                                                                                                                                     <  7:92684/145‑1 (MQ=255)
     gCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGa                                                                                                                                             >  1:139121/1‑148 (MQ=255)
      cAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGaca                                                                                                                                           >  5:301149/1‑149 (MQ=255)
              tttGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGc                                                                                                                                     >  5:157551/1‑147 (MQ=255)
               ttGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGa                                                                                                                                   <  2:43241/148‑1 (MQ=255)
                 gCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCGAGGGCCTGACACTGGGCGAGCt                                                                                                                                >  2:295391/1‑149 (MQ=255)
                 gCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCGAGGGCCTGACACTGGGCGAGCt                                                                                                                                >  2:187351/1‑149 (MQ=255)
                     tgtgAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGca                                                                                                                            >  1:70878/1‑149 (MQ=255)
                       tgAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGca                                                                                                                            >  3:101813/1‑147 (MQ=255)
                        gAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAg                                                                                                                         >  4:294978/1‑149 (MQ=255)
                                  cTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGcc                                                                                                               >  2:64871/1‑149 (MQ=255)
                                   tGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCg                                                                                                              >  3:61735/1‑149 (MQ=255)
                                   tGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCg                                                                                                              >  1:302737/1‑149 (MQ=255)
                                     tACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGtt                                                                                                            <  3:230329/149‑1 (MQ=255)
                                       cTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCg                                                                                                          <  2:302737/149‑1 (MQ=255)
                                        ttttCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCgg                                                                                                         >  6:261135/1‑149 (MQ=255)
                                        ttttCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTAAAGGGCCTGACACTGGGCGATCTGGCACAGCGCATCTGCCGTTCgg                                                                                                         >  7:179177/1‑149 (MQ=255)
                                          ttCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGt                                                                                                        <  1:64871/148‑1 (MQ=255)
                                           tCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTg                                                                                                       <  2:228977/148‑1 (MQ=255)
                                                gCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCttt                                                                                                 <  6:301149/149‑1 (MQ=255)
                                                gCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCttt                                                                                                 <  4:101813/149‑1 (MQ=255)
                                                  aTCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGATGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCt                                                                                               >  6:8167/1‑149 (MQ=255)
                                                     cATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGt                                                                                             <  7:47868/148‑1 (MQ=255)
                                                     cATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTc                                                                                            >  1:103382/1‑149 (MQ=255)
                                                          cccGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAg                                                                                        <  6:157551/148‑1 (MQ=255)
                                                                  aaaGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAAc                                                                                 <  2:103382/147‑1 (MQ=255)
                                                                  aaaGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGACCCA‑GTCGAAc                                                                                 <  2:70878/146‑1 (MQ=255)
                                                                   aaGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCGAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAAcgcg                                                                              <  1:187351/149‑1 (MQ=255)
                                                                   aaGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCGAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAAcgcg                                                                              <  1:295391/149‑1 (MQ=255)
                                                                    aGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACgcgg                                                                             <  3:294978/149‑1 (MQ=255)
                                                                     gTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACgcggc                                                                            >  6:99785/1‑149 (MQ=255)
                                                                        gcgATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCgtg                                                                         >  2:246130/1‑149 (MQ=255)
                                                                             tCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCCCg                                                                    >  4:255404/1‑149 (MQ=255)
                                                                               aaCCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGTCGTGTCACGcc                                                                  >  8:292676/1‑149 (MQ=255)
                                                                               aaCCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGcc                                                                  >  1:246516/1‑149 (MQ=255)
                                                                                    tGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGAcc                                                             >  4:121356/1‑149 (MQ=255)
                                                                                     gAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTgg                                                                                                      <  4:54624/107‑1 (MQ=255)
                                                                                         gTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTTGAACGCGGCGTGTCACGCCCGACCGTGGc                                                        <  2:139121/149‑1 (MQ=255)
                                                                                             aCGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGAt                                                    >  4:109685/1‑149 (MQ=255)
                                                                                             aCGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGAt                                                    <  6:39999/149‑1 (MQ=255)
                                                                                               gaagaaATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGAtct                                                  >  5:100868/1‑149 (MQ=255)
                                                                                                aagaaATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGAtct                                                  >  1:31693/1‑148 (MQ=255)
                                                                                                 agaaATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGAt                                                    <  1:145043/145‑1 (MQ=255)
                                                                                                 agaaATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGAt                                                    >  2:145043/1‑145 (MQ=255)
                                                                                                   aaaTCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTcaca                                              <  5:99785/149‑1 (MQ=255)
                                                                                                        gAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCa                                          <  5:8167/148‑1 (MQ=255)
                                                                                                        gAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCAt                                         <  2:200042/149‑1 (MQ=255)
                                                                                                         aaGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATc                                        <  5:261135/149‑1 (MQ=255)
                                                                                                          aGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCa                                       <  2:31693/149‑1 (MQ=255)
                                                                                                            gCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGc                                     <  3:121356/149‑1 (MQ=255)
                                                                                                                 gCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGAtc                                                   >  8:256755/1‑130 (MQ=255)
                                                                                                                 gCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGAtc                                                   <  7:256755/130‑1 (MQ=255)
                                                                                                                 gCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGCGaa                                  <  2:246516/147‑1 (MQ=255)
                                                                                                                 gCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGCGAaga                                <  6:127763/149‑1 (MQ=255)
                                                                                                                        tGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTca                                                >  3:240224/1‑126 (MQ=255)
                                                                                                                        tGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTca                                                <  4:240224/126‑1 (MQ=255)
                                                                                                                        tGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGCGAAGAGCTCGgt                         >  1:143/1‑149 (MQ=255)
                                                                                                                          aTTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCCCAGCCATCAGCGaa                                  <  7:262144/138‑1 (MQ=255)
                                                                                                                          aTTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGCGaa                                  >  8:262144/1‑138 (MQ=255)
                                                                                                                           ttCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGCGAAGAGCTCGGTGTa                      >  8:182927/1‑149 (MQ=255)
                                                                                                                                             tCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGCGAAGAGCTCGGTGTATCCACAGCCTACTTCTa     >  8:25275/1‑148 (MQ=255)
                                                                                                                                              cAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGCGAAGAGCTCTGTGTATCCACAGCCTACTTCTAc    <  2:143/148‑1 (MQ=255)
                                                                                                                                               aaGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGCGAAGAGCTCGGTGTATCCACAGCCTACTTCTAc    >  3:299705/1‑147 (MQ=255)
                                                                                                                                               aaGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCAAAGCCATCAGCGAAGAGCTCGGTGTATCCACAGCCTACTTCTACaa  <  7:292676/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAAGTGCACGAAGAAATCGAAGGCCTGGCAGTACTGATTCGCGACCTGCGCAAGTTCAAGGGCCTGACACTGGGCGAGCTGGCACAGCGCATCGGCCGTTCGGTGGGCTTTCTGTCCCAGGTCGAACGCGGCGTGTCACGCCCGACCGTGGCCGATCTCACAGCCATCAGCGAAGAGCTCGGTGTATCCACAGCCTACTTCTACAA  >  NC_002947/2486868‑2487159

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