Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,131,154 T→A 7.0% intergenic (+22/+90) PP_3634 → / ← PP_3635 GNAT family acetyltransferase/sulfonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,131,1540TA 7.0% 286.0 / 9.2 100intergenic (+22/+90)PP_3634/PP_3635GNAT family acetyltransferase/sulfonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (49/41);  new base A (3/4);  total (53/47)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 107 total aligned reads displayed.

CACCTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAG  >  NC_002947/4131006‑4131303
                                                                                                                                                    |                                                                                                                                                     
cacCTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCtt                                                                                                                                                      <  2:37270/150‑1 (MQ=255)
cacCTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCtt                                                                                                                                                      >  2:33059/1‑150 (MQ=255)
 acCTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATAGAAGAGGAGAAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTg                                                                                                                                                     >  1:37270/1‑150 (MQ=255)
   cTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTgag                                                                                                                                                   >  2:432067/1‑148 (MQ=255)
   cTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCGGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTGGCGGGCTTGcc                                                                                                                                                   >  1:857030/1‑150 (MQ=255)
      tCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCcgc                                                                                                                                                >  2:1900310/1‑150 (MQ=255)
          cTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGccc                                                                                                                                            >  2:850693/1‑150 (MQ=255)
           ttAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCca                                                                                                                                           >  1:667038/1‑150 (MQ=255)
           ttAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgcgaaca                                                                                                                                           <  1:1366060/150‑5 (MQ=255)
             aaCCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCcaca                                                                                                                                         >  2:1122470/1‑150 (MQ=255)
              aCCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAg                                                                                                                                        <  2:667038/150‑1 (MQ=255)
              aCCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAg                                                                                                                                        <  1:492535/150‑1 (MQ=255)
                 agcagGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGaa                                                                                                                                     >  2:57443/1‑150 (MQ=255)
                 agcagGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCGCAAGCCCAGGAATGTGTGCGCCAGCGCAAGACCGTGCATTGATCCCATGGCCTGTTCGCGGGCGTGCCCGCGCCCACAGGaa                                                                                                                                     >  1:695270/1‑150 (MQ=255)
                    agGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATAc                                                                                                                                  >  2:1877256/1‑150 (MQ=255)
                       gTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGt                                                                                                                               <  1:850693/150‑1 (MQ=255)
                             aTGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGcc                                                                                                                         <  2:327695/150‑1 (MQ=255)
                                 aaTGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAg                                                                                                                     >  2:1238988/1‑150 (MQ=255)
                                     aTCGAAGCGGAGCAGATGCGCGATGGGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCg                                                                                                                 >  1:1484209/1‑150 (MQ=255)
                                         aaGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCGCGCCCATGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGa                                                                                                             <  2:1484209/150‑1 (MQ=255)
                                         aaGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGa                                                                                                             <  2:1796822/150‑1 (MQ=255)
                                          aGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGAt                                                                                                            <  1:1900310/150‑1 (MQ=255)
                                             ggAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTa                                                                                                         <  2:1115167/150‑1 (MQ=255)
                                              gAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACACGAATACCGTACACCTCAAGGCCGGTGATGTAc                                                                                                        >  1:97555/1‑150 (MQ=255)
                                                  aGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCtgt                                                                                                    <  1:57443/150‑1 (MQ=255)
                                                    aTGCGCGATTTTTGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTgg                                                                                                  <  1:1515637/150‑1 (MQ=255)
                                                        gcgATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGc                                                                                              <  1:1687735/150‑1 (MQ=255)
                                                        gcgATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgc                                                                                                                                                >  2:1920102/1‑100 (MQ=255)
                                                        gcgATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCAAGCCcgc                                                                                                                                                <  1:1920102/100‑1 (MQ=255)
                                                          gATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCgg                                                                                            <  2:27116/150‑1 (MQ=255)
                                                              ttcGGCTGTATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGa                                                                                        <  2:205580/148‑1 (MQ=255)
                                                               gctgCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAg                                                                                       >  2:1222323/4‑150 (MQ=255)
                                                                 ggCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGcc                                                                                     <  1:1777621/150‑1 (MQ=255)
                                                                      gATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCga                                                                                >  1:141339/1‑150 (MQ=255)
                                                                       aTGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCgaa                                                                               >  1:325865/1‑150 (MQ=255)
                                                                       aTGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCgaa                                                                               <  1:1877256/150‑1 (MQ=255)
                                                                        tGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCgaag                                                                              >  2:1405907/1‑150 (MQ=255)
                                                                          aGCTGACCGCTCAAGCCCAGGAATGTGTGCGTTAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCgaagaa                                                                            <  1:1151193/150‑1 (MQ=255)
                                                                           gCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCgg                                                                                                                <  2:701875/113‑1 (MQ=255)
                                                                           gCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGc                                                                                      <  2:739274/139‑1 (MQ=255)
                                                                           gCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGc                                                                                      >  1:739274/1‑139 (MQ=255)
                                                                            cTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGc                                                                          >  2:29320/1‑150 (MQ=255)
                                                                             ttaCCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGCCCTGCCACCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcc                                                                         <  2:366042/148‑1 (MQ=255)
                                                                             tGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcc                                                                         <  1:806232/150‑1 (MQ=255)
                                                                              gACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcc                                                                         >  1:570009/1‑149 (MQ=255)
                                                                              gACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcc                                                                         <  2:570009/149‑1 (MQ=255)
                                                                               aCCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCgg                                                                                                                >  1:1256955/1‑109 (MQ=255)
                                                                               aCCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGcccc                                                                       >  2:349796/1‑150 (MQ=255)
                                                                                 cGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATa                                                                                                                                   >  2:1389263/1‑88 (MQ=255)
                                                                                 cGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATa                                                                                                                                   <  1:1389263/88‑1 (MQ=255)
                                                                                   cTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTg                                                                   <  1:1405907/150‑1 (MQ=255)
                                                                                    tCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGCGAGCGGGCGAGCCCGCGAAGAAGCCCCATTgg                                                                  >  2:1576194/1‑150 (MQ=255)
                                                                                          ccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCtcgctcgct                                                                                                                                               <  2:495761/67‑2 (MQ=255)
                                                                                          ccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGTAGCGGGCGAGCCCGCGCAGAAGCCCCATTGGCCTTTc                                                            >  2:83031/1‑150 (MQ=255)
                                                                                          ccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTc                                                            <  2:1473565/150‑1 (MQ=255)
                                                                                          ccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTc                                                            >  2:1315147/1‑150 (MQ=255)
                                                                                          ccAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTCGCCcgct                                                                                                                                               >  1:495761/1‑66 (MQ=255)
                                                                                           cAGGAATGTGTGCGCCAGCGCAAGACCCACCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                           >  1:1790360/1‑150 (MQ=255)
                                                                                           cAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                           <  1:1238988/150‑1 (MQ=255)
                                                                                           cAGGAATGGGTGCGCCAGCGCAAGAACGTGCATTGATGCGATGGCCTGTTCGCACGCAGGGTCGCGCCCACAGGAATAGCGTAAGGCCCAAGGCCGGTGATGTACTTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                           <  2:258499/150‑1 (MQ=255)
                                                                                             ggAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGc                                                         <  1:1759756/150‑1 (MQ=255)
                                                                                             ggAATGTGTGCGCCACCGCAAGAACGCGCATTGATCCGCTGCTCTATTAGCGATTGATGATGGTTTCGCGGGGAGACTGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGc                                                         <  1:432067/150‑1 (MQ=255)
                                                                                              gAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGcc                                                        >  2:743688/1‑150 (MQ=255)
                                                                                               aaTGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCa                                                       <  1:743688/150‑1 (MQ=255)
                                                                                               aaTGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCa                                                       <  1:349796/150‑1 (MQ=255)
                                                                                                aTGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCtgt                                                                                                    >  1:866754/1‑104 (MQ=255)
                                                                                                aTGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCtgt                                                                                                    >  1:866753/1‑104 (MQ=255)
                                                                                                aTGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCtgt                                                                                                    <  2:866753/104‑1 (MQ=255)
                                                                                                aTGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCtgt                                                                                                    <  2:866754/104‑1 (MQ=255)
                                                                                                    gtgCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                           >  1:1102535/1‑141 (MQ=255)
                                                                                                    gtgCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCa                                                           <  2:1102535/141‑1 (MQ=255)
                                                                                                           aGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                               >  1:1684595/1‑82 (MQ=255)
                                                                                                           aGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGt                                                                                                               <  2:1684595/82‑1 (MQ=255)
                                                                                                              gcAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGGGCCCGCGAAGAAGCCCCATTGGCCTTTCAACCATGCACCCAGCGCCGa                                        >  1:298834/1‑150 (MQ=255)
                                                                                                              gcAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGa                                        <  1:1222323/150‑1 (MQ=255)
                                                                                                               cAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGAt                                       <  1:1234194/150‑1 (MQ=255)
                                                                                                                      gTGCATTGAGCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGccc                                >  2:1713731/1‑150 (MQ=255)
                                                                                                                        gCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAgcg                                            >  2:1409677/1‑136 (MQ=255)
                                                                                                                        gCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAgcg                                            <  1:1409677/136‑1 (MQ=255)
                                                                                                                         cATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCAcg                             >  1:1730170/1‑150 (MQ=255)
                                                                                                                          aTTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCACGCACCCAGCGCCGATGCAGCCCAcgc                            >  2:1411927/1‑150 (MQ=255)
                                                                                                                              aTCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGCCCGGTGATGTACCTGTGGGAACGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCa                        >  2:566863/1‑150 (MQ=255)
                                                                                                                                ccGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCaga                      <  1:443615/150‑1 (MQ=255)
                                                                                                                                 cGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCt                                                                                                      >  2:1079914/1‑69 (MQ=255)
                                                                                                                                 cGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCt                                                                                                      <  1:1079914/69‑1 (MQ=255)
                                                                                                                                  gATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGc                    >  2:1639954/1‑150 (MQ=255)
                                                                                                                                      gCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATcc                >  2:131196/1‑150 (MQ=255)
                                                                                                                                       cgtgggcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGc                                  <  2:1558007/125‑1 (MQ=39)
                                                                                                                                       cgtgggcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGc                                  >  1:1558007/7‑131 (MQ=39)
                                                                                                                                       ccTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCCTGCACCCAGCGCCGATGCAGCCCGCGCGCCAGAGCATCCa               >  2:971192/1‑150 (MQ=255)
                                                                                                                                       ccTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCa               >  1:847167/1‑150 (MQ=255)
                                                                                                                                        cTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGAAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAg              >  1:1892766/1‑150 (MQ=255)
                                                                                                                                          gTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGca            <  1:1122470/150‑1 (MQ=255)
                                                                                                                                            gcgcgGGCAAGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGcaca          >  1:1692438/2‑150 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGc                                                                          >  1:1460563/1‑84 (MQ=255)
                                                                                                                                              gcgGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGc                                                                          <  2:1460563/84‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAAccc     <  2:606047/147‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAAccc     >  1:606047/1‑147 (MQ=255)
                                                                                                                                                    ttGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAg  >  1:145294/1‑150 (MQ=255)
                                                                                                                                                    ttGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAg  >  1:685184/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACCTATCTGCTTAACCAGCAGGGTGAGGATGCAATGATCGAAGCGGAGCAGATGCGCGATGAGCGGCTGGATGAGCTGACCGCTCAAGCCCAGGAATGTGTGCGCCAGCGCAAGAACGTGCATTGATCCGATGGCCTGTTCGCGGGCTTGCCCGCGCCCACAGGAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAG  >  NC_002947/4131006‑4131303

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