Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,515 G→T 5.5% R134S (CGT→AGT) ‡ amaB ← L‑piperidine‑6‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,5150GT 5.5% 790.9 / 24.3 236R134S (CGT→AGT) ‡amaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (106/117);  new base T (4/9);  total (110/126)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 237 total aligned reads displayed.

AGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGC  >  NC_002947/6005367‑6005664
                                                                                                                                                    |                                                                                                                                                     
aGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGc                                                                                                                                                      <  2:4894276/150‑1 (MQ=255)
     ggTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCAtgtg                                                                                                                                                 <  2:2611284/150‑1 (MQ=255)
       tCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTAGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCAggtggt                                                                                                                                               >  1:3881009/1‑150 (MQ=255)
        cAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATgtggtg                                                                                                                                              <  2:3852019/150‑1 (MQ=255)
         aGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGtggtgg                                                                                                                                             <  1:3993049/150‑1 (MQ=255)
           cGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGcc                                                                                                                                           >  1:480031/1‑150 (MQ=255)
             gAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCg                                                                                                                                         >  2:1860572/1‑150 (MQ=255)
              aGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCCTGTGGTGGCCCgg                                                                                                                                        <  2:4056032/150‑1 (MQ=255)
              aGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGAGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCgg                                                                                                                                        >  2:2636964/1‑150 (MQ=255)
                tcttctCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGgc                                                                                                                                      >  2:524727/1‑150 (MQ=255)
                tcttctCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGgc                                                                                                                                      >  2:3023304/1‑150 (MQ=255)
                  ttctCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGgcgc                                                                                                                                    >  2:503025/1‑150 (MQ=255)
                        gACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGaa                                                                                                                              <  1:1574389/150‑1 (MQ=255)
                         aCGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAg                                                                                                                             <  2:38007/150‑1 (MQ=255)
                           ggTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCg                                                                                                                           >  1:1762400/1‑150 (MQ=255)
                            gTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGa                                                                                                                          <  1:4411181/150‑1 (MQ=255)
                             tttCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGAt                                                                                                                         >  1:3373986/1‑150 (MQ=255)
                              ttCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATg                                                                                                                        <  1:3210675/150‑1 (MQ=255)
                               tCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATgg                                                                                                                       <  2:2580759/150‑1 (MQ=255)
                                 cATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTc                                                                                                                     >  2:4495770/1‑150 (MQ=255)
                                  aTACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCa                                                                                                                    >  2:3433719/1‑150 (MQ=255)
                                   tACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAg                                                                                                                   <  2:4300585/150‑1 (MQ=255)
                                      caccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGcc                                                                                                                >  2:3354034/1‑150 (MQ=255)
                                      caccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGcc                                                                                                                >  1:3952146/1‑150 (MQ=255)
                                       accGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCg                                                                                                               >  1:2936382/1‑150 (MQ=255)
                                        ccGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGt                                                                                                              <  1:52005/150‑1 (MQ=255)
                                         cGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCa                                                                                                                    <  2:3447489/143‑1 (MQ=255)
                                            gTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCATCCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTActg                                                                                                          <  2:2927481/150‑3 (MQ=255)
                                              tGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCt                                                                                                        <  1:3438912/150‑1 (MQ=255)
                                              tGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCt                                                                                                        <  1:1289805/150‑1 (MQ=255)
                                               gCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTg                                                                                                       <  2:3971036/150‑1 (MQ=255)
                                                ccGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTgg                                                                                                      >  1:2724596/1‑150 (MQ=255)
                                                 cGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGgc                                                                                                     >  2:198348/1‑150 (MQ=255)
                                                         agggccagggcGGTGTTCCAGTCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAgg                                                                                             <  2:4309737/150‑1 (MQ=255)
                                                         agggccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAgg                                                                                             <  1:2891912/150‑1 (MQ=255)
                                                             ccagggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCaa                                                                                         <  2:1742531/150‑1 (MQ=255)
                                                                 ggcGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGc                                                                                     <  1:1119247/150‑1 (MQ=255)
                                                                    ggTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGaa                                                                                  <  1:4879315/150‑1 (MQ=255)
                                                                    ggTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGaa                                                                                  <  2:1496724/150‑1 (MQ=255)
                                                                    ggTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGaa                                                                                  <  2:2758596/150‑1 (MQ=255)
                                                                              gCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATg                                                                        >  1:1066263/1‑150 (MQ=255)
                                                                               cccACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGt                                                                       <  1:3433712/150‑1 (MQ=255)
                                                                                 cacaCGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCg                                                                     >  2:4174233/1‑150 (MQ=255)
                                                                                   cacGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGAt                                                                   <  1:2537286/150‑1 (MQ=255)
                                                                                         aaCCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTc                                                             <  2:2689243/150‑1 (MQ=255)
                                                                                         aaCAGGGAAGTTGAAGGCCGAGATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTc                                                             <  2:4472492/150‑1 (MQ=255)
                                                                                            cGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTg                                                          >  1:605423/1‑150 (MQ=255)
                                                                                            cGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTg                                                          <  2:4805304/150‑1 (MQ=255)
                                                                                            cGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTg                                                          <  2:3739986/150‑1 (MQ=255)
                                                                                            cGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTg                                                          >  1:3162165/1‑150 (MQ=255)
                                                                                            agggAAGTTGAAGGCCGAGATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTg                                                          >  1:557770/2‑150 (MQ=255)
                                                                                                   ttGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCa                                                   <  1:98885/150‑1 (MQ=255)
                                                                                                     gAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCAcc                                                 >  1:1866340/1‑150 (MQ=255)
                                                                                                     gAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCAcc                                                 >  2:3122243/1‑150 (MQ=255)
                                                                                                       aGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCa                                               <  2:2544156/150‑1 (MQ=255)
                                                                                                        ggCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAg                                              >  1:2449985/1‑150 (MQ=255)
                                                                                                         gcgcTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAgg                                             >  1:3363956/1‑150 (MQ=255)
                                                                                                           gcTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGcc                                           >  2:1063328/1‑150 (MQ=255)
                                                                                                            cTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCt                                          >  2:1911033/1‑150 (MQ=255)
                                                                                                             tGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCtt                                         >  1:2611043/1‑150 (MQ=255)
                                                                                                             tGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCtt                                         <  1:4054158/150‑1 (MQ=255)
                                                                                                             agATAACGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCtt                                         <  1:3798194/149‑1 (MQ=255)
                                                                                                                 cACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGTGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTg                                     <  1:4590594/150‑1 (MQ=255)
                                                                                                                  aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                    <  2:3307230/150‑1 (MQ=255)
                                                                                                                  aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                    <  1:766306/150‑1 (MQ=255)
                                                                                                                  aCGCCGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTccgaa                                    <  2:4778589/150‑4 (MQ=255)
                                                                                                                   cGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTggg                                   <  2:4716350/150‑1 (MQ=255)
                                                                                                                    gCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGt                                  >  2:4242978/1‑150 (MQ=255)
                                                                                                                      cGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGa                                <  2:1395940/150‑1 (MQ=255)
                                                                                                                      cGACAACGCCAAGCGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTccgaaggaa                                <  2:1132350/150‑8 (MQ=255)
                                                                                                                        aCCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATc                              <  2:4271826/150‑1 (MQ=255)
                                                                                                                        aCCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATc                              <  2:1117762/150‑1 (MQ=255)
                                                                                                                        aCCACACCCACCGGGTGCCAGGTTTCACGCAAGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTAAAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGAGCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATc                              >  2:2859687/1‑150 (MQ=255)
                                                                                                                          cacaCCCAGCGGGTGCCTGGTTTGACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGACGTACAGCTGGCGCTACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCTGGCCTTCCTGGGTGATCtt                            <  2:1315717/150‑1 (MQ=255)
                                                                                                                          cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                            <  1:1654478/150‑1 (MQ=255)
                                                                                                                          cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                            >  1:3146112/1‑150 (MQ=255)
                                                                                                                          cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                            >  2:4182793/1‑150 (MQ=255)
                                                                                                                          cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                            <  2:2036565/150‑1 (MQ=255)
                                                                                                                             acCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGcc                         >  1:3652105/1‑150 (MQ=255)
                                                                                                                             acCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGcc                         >  1:4199431/1‑150 (MQ=255)
                                                                                                                             acCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGcc                         <  1:2209798/150‑1 (MQ=255)
                                                                                                                               ccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCgg                       >  1:2257532/1‑150 (MQ=255)
                                                                                                                                cAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGc                      <  1:1455182/150‑1 (MQ=255)
                                                                                                                                 aGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCt                     >  2:1394771/1‑150 (MQ=255)
                                                                                                                                 aGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCt                     <  1:4142522/150‑1 (MQ=255)
                                                                                                                                  gCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCtt                    <  1:4495770/150‑1 (MQ=255)
                                                                                                                                     ggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcga                 <  1:2028698/150‑1 (MQ=255)
                                                                                                                                     ggTGCCAGCTTTCACGCATGTGGTGGCCCGCGCGCTCGCAAGCGATGGTCAGGCCGTACAGCTGGCGCGACACGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCACGCCTTCCTGGGTGAACTTGCCGGCTtcga                 <  2:4225076/150‑1 (MQ=255)
                                                                                                                                     gatGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTtcga                 >  1:705349/3‑150 (MQ=255)
                                                                                                                                     gatGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTtcga                 <  2:360364/148‑1 (MQ=255)
                                                                                                                                      gTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgat                <  2:2170059/150‑1 (MQ=255)
                                                                                                                                       tGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatc               <  2:619250/150‑1 (MQ=255)
                                                                                                                                         ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga             <  1:806372/150‑1 (MQ=255)
                                                                                                                                         ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga             >  2:1611133/1‑150 (MQ=255)
                                                                                                                                          cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag            <  2:4371448/150‑1 (MQ=255)
                                                                                                                                            ggTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGAc          <  2:3514562/150‑1 (MQ=255)
                                                                                                                                              tttCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCa        >  2:1729238/1‑150 (MQ=255)
                                                                                                                                                 cACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCt     <  2:5070899/150‑1 (MQ=255)
                                                                                                                                                  aCGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTc    >  1:4093920/1‑150 (MQ=255)
                                                                                                                                                    gCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGc  >  2:2414348/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGGCGGTCAGCGGAGTCTTCTCGGACGGTTTCCATACCACCGAGTTGCCGGCCACCAGGGCCAGGGCGGTGTTCCAGGCCCACACGGCAACCGGGAAGTTGAAGGCGCTGATCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGC  >  NC_002947/6005367‑6005664

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: