Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,647,605 C→T 5.7% Q185* (CAG→TAG) ‡ PP_4107 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,647,6050CT 5.7% 1297.8 / 40.8 401Q185* (CAG→TAG) ‡PP_4107LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (183/195);  new base T (14/9);  total (197/204)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 408 total aligned reads displayed.

ATTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCATCTCC  >  NC_002947/4647457‑4647754
                                                                                                                                                    |                                                                                                                                                     
atacGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGACGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCa                                                                                                                                                      >  2:1758328/4‑150 (MQ=255)
aTTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCa                                                                                                                                                      >  1:6055431/1‑150 (MQ=255)
 ttCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAg                                                                                                                                                     >  1:4128478/1‑150 (MQ=255)
 ttCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAg                                                                                                                                                     >  2:4179867/1‑150 (MQ=255)
 ttCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAg                                                                                                                                                     >  2:3951343/1‑150 (MQ=255)
 ttCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAg                                                                                                                                                     >  2:6406967/1‑150 (MQ=255)
     aCCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGgcgc                                                                                                                                                 >  1:2282953/1‑150 (MQ=255)
        tGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgt                                                                                                                                              <  2:6302990/150‑2 (MQ=255)
        tGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgt                                                                                                                                              <  2:5985343/150‑2 (MQ=255)
        tGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctga                                                                                                                                              >  2:1865223/1‑150 (MQ=255)
         ggTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgtc                                                                                                                                             <  2:665077/150‑3 (MQ=255)
          gTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGTCAGCGCCTGGGCGctgtcg                                                                                                                                            >  1:5371411/1‑147 (MQ=255)
          gTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgact                                                                                                                                            >  1:2486504/1‑150 (MQ=255)
          gTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgact                                                                                                                                            <  2:7221969/150‑1 (MQ=255)
          gTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgact                                                                                                                                            <  2:5891428/150‑1 (MQ=255)
          gTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgact                                                                                                                                            <  1:3987739/150‑1 (MQ=255)
            aCAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgactga                                                                                                                                          <  1:2582951/150‑1 (MQ=255)
            aCAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgactga                                                                                                                                          <  2:7292706/150‑1 (MQ=255)
            aCAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgactga                                                                                                                                          <  1:5739061/150‑1 (MQ=255)
             cAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGTCAGCGCCTGGGCGctgtcgaaa                                                                                                                                         <  2:5371411/150‑7 (MQ=255)
              aGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGc                                                                                                                                        >  1:7266282/1‑150 (MQ=255)
                gCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCAc                                                                                                                                      >  1:7201198/1‑150 (MQ=255)
                gCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCAc                                                                                                                                      >  1:4575784/1‑150 (MQ=255)
                      gCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTg                                                                                                                                <  2:911118/150‑1 (MQ=255)
                        ggACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGtt                                                                                                                              <  1:3020018/150‑1 (MQ=255)
                        ggACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGtt                                                                                                                              <  2:803708/150‑1 (MQ=255)
                         gACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTg                                                                                                                             >  2:52461/1‑150 (MQ=255)
                               gCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCTGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGc                                                                                                                       >  1:1464590/1‑150 (MQ=255)
                                 cTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCgt                                                                                                                     >  1:3395866/1‑150 (MQ=255)
                                  tttttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCgtg                                                                                                                    >  2:1124583/1‑150 (MQ=255)
                                   ttttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCgtgt                                                                                                                   <  1:7128152/150‑1 (MQ=255)
                                      tcttcGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTcgc                                                                                                                <  2:1305056/150‑1 (MQ=255)
                                        ttcGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAg                                                                                                                                         <  2:4101395/123‑1 (MQ=255)
                                        ttcGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCgg                                                                                                              <  2:3130282/150‑1 (MQ=255)
                                         tcGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGc                                                                                                             <  1:7418999/150‑1 (MQ=255)
                                            gCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCgg                                                                                                          >  2:5831092/1‑150 (MQ=255)
                                                ggtgcccTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGc                                                                                                      >  1:6753095/6‑150 (MQ=255)
                                                gggCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGc                                                                                                      <  2:2747969/150‑1 (MQ=255)
                                                 ggCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGcc                                                                                                     >  1:3718317/1‑150 (MQ=255)
                                                   cACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTg                                                                                                   <  2:6506831/150‑1 (MQ=255)
                                                      cTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCa                                                                                                <  2:515625/150‑1 (MQ=255)
                                                      cTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCa                                                                                                >  1:5665096/1‑150 (MQ=255)
                                                        ggCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATg                                                                                              >  1:3465973/1‑150 (MQ=255)
                                                         gCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGa                                                                                             <  2:2993871/150‑1 (MQ=255)
                                                          cccGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGaa                                                                                            <  2:2259245/150‑1 (MQ=255)
                                                           ccGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAAt                                                                                           >  2:5051746/1‑150 (MQ=255)
                                                             ggggCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATgg                                                                                         >  1:4897446/1‑150 (MQ=255)
                                                                     cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTggg                                                                                 >  1:4815236/1‑150 (MQ=255)
                                                                     cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTggg                                                                                 >  2:1443055/1‑150 (MQ=255)
                                                                           cGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTgg                                                                           <  1:5048237/150‑1 (MQ=255)
                                                                               ttgttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTg                                                                       <  1:4698529/150‑1 (MQ=255)
                                                                                 gttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCa                                                                     <  2:4016504/150‑1 (MQ=255)
                                                                                   tCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATa                                                                   <  1:7053798/150‑1 (MQ=255)
                                                                                   tCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCACAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATa                                                                   <  2:471205/150‑1 (MQ=255)
                                                                                     aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                 >  1:6315720/1‑150 (MQ=255)
                                                                                     aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                 >  1:4013512/1‑150 (MQ=255)
                                                                                        gaagaaGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTg                                                                                                   <  1:795944/113‑1 (MQ=255)
                                                                                        gaagaaGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAg                                                              <  1:99825/150‑1 (MQ=255)
                                                                                        gaagaaGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCATTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTgg                                                                           <  2:242122/137‑1 (MQ=255)
                                                                                          agaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGaa                                                            <  1:1809839/150‑1 (MQ=255)
                                                                                          agaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCACCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGaa                                                            <  2:6832090/150‑1 (MQ=255)
                                                                                           gaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGtt                                                                        <  1:6718617/137‑1 (MQ=255)
                                                                                            aagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACa                                                          <  1:1155109/150‑1 (MQ=255)
                                                                                            aagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACa                                                          >  1:2399585/1‑150 (MQ=255)
                                                                                            aagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACa                                                          <  1:1370199/150‑1 (MQ=255)
                                                                                                  tgCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTg                                                                                   >  2:3306430/1‑119 (MQ=255)
                                                                                                            cAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCcgcg                                          >  1:7341832/1‑150 (MQ=255)
                                                                                                            cAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCcgcg                                          <  2:1545719/150‑1 (MQ=255)
                                                                                                            cAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCcgcg                                          >  1:4484322/1‑150 (MQ=255)
                                                                                                             aGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCgcgc                                         >  2:4428633/1‑150 (MQ=255)
                                                                                                                ccGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTc                                      <  1:3536946/150‑1 (MQ=255)
                                                                                                                ccGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTc                                      >  1:2219299/1‑150 (MQ=255)
                                                                                                                ccGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTc                                      <  2:7248829/150‑1 (MQ=255)
                                                                                                                   cAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGac                                   <  1:6221470/150‑1 (MQ=255)
                                                                                                                     gctgctGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGacac                                 <  1:5876381/150‑1 (MQ=255)
                                                                                                                       tgctgGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGt                               <  1:1758328/150‑1 (MQ=255)
                                                                                                                       tgctgGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGt                               <  2:551630/150‑1 (MQ=255)
                                                                                                                        gctgGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGtt                              <  2:2032666/150‑1 (MQ=255)
                                                                                                                         ctgGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTc                             >  1:4979298/1‑150 (MQ=255)
                                                                                                                            gccagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTAc                          <  2:7126625/150‑1 (MQ=255)
                                                                                                                            gccagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTAc                          <  2:5568885/150‑1 (MQ=255)
                                                                                                                            gccagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTAc                          >  1:1022891/1‑150 (MQ=255)
                                                                                                                               agccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtg                       <  2:4575784/150‑1 (MQ=255)
                                                                                                                               agccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtg                       >  2:3518935/1‑150 (MQ=255)
                                                                                                                               agccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtg                       >  1:3211706/1‑150 (MQ=255)
                                                                                                                                gccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGGCGCGCTTCGACACGTTCTACCtgc                      <  1:2895278/150‑2 (MQ=255)
                                                                                                                                gccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                      <  1:4188558/150‑1 (MQ=255)
                                                                                                                                gccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                      <  2:2552023/150‑1 (MQ=255)
                                                                                                                                  caCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTgtgt                    >  2:6012950/1‑150 (MQ=255)
                                                                                                                                      gTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATcc                >  1:6596301/1‑150 (MQ=255)
                                                                                                                                         tCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCcgcg             <  2:4817153/150‑1 (MQ=255)
                                                                                                                                          cGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCcgcgc            >  2:7362561/1‑150 (MQ=255)
                                                                                                                                          agtcAGCGCCTGGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCTTAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCcgcgc            >  2:6955952/4‑150 (MQ=255)
                                                                                                                                           gACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCcgcgcg           >  1:3325259/1‑150 (MQ=255)
                                                                                                                                           gACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCcgcgcg           >  1:1975106/1‑150 (MQ=255)
                                                                                                                                            tcAGCGCCTGGGCGCTGTCTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCgg          >  2:3770159/2‑150 (MQ=255)
                                                                                                                                            aCAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCgg          >  2:1613109/1‑150 (MQ=255)
                                                                                                                                               gcgcCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCa       >  1:5233629/1‑150 (MQ=255)
                                                                                                                                                   ccAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGc                                                                      <  1:358234/83‑1 (MQ=255)
                                                                                                                                                   ccAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGc                                                                      >  2:358234/1‑83 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCATCTCC  >  NC_002947/4647457‑4647754

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: