Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,316,740 A→G 5.3% S263P (TCT→CCT)  nusA ← transcription termination/antitermination L factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,316,7400AG 5.3% 1033.4 / 30.9 301S263P (TCT→CCT) nusAtranscription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base A (145/140);  new base G (8/8);  total (153/148)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 309 total aligned reads displayed.

CCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCA  >  NC_000913/3316592‑3316889
                                                                                                                                                    |                                                                                                                                                     
ccGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGac                                                                                                                                                      <  1:3798468/150‑1 (MQ=255)
ccGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGac                                                                                                                                                      <  2:2390442/150‑1 (MQ=255)
ccGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGac                                                                                                                                                      >  2:2906006/1‑150 (MQ=255)
 cGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGaca                                                                                                                                                     >  2:764828/1‑150 (MQ=255)
 cGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGAGCCGGGTTATCATCCCACAGGACGAAATCGATACGCTCGCCACCCAGTTCAGTAGaca                                                                                                                                                     >  1:772931/1‑150 (MQ=255)
   gCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGGCAcc                                                                                                                                                   >  2:605246/1‑150 (MQ=255)
     ttCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGc                                                                                                                                                 >  1:1936886/1‑150 (MQ=255)
     ttCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGc                                                                                                                                                 <  1:2178577/150‑1 (MQ=255)
       cAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCt                                                                                                                                               >  1:2524292/1‑150 (MQ=255)
           ggCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAac                                                                                                                                           <  1:1826520/150‑1 (MQ=255)
             cGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAacac                                                                                                                                         >  2:1617275/1‑150 (MQ=255)
                 atCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACAcgcgc                                                                                                                                     >  2:2018994/1‑150 (MQ=255)
                  tCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACAcgcgcg                                                                                                                                    >  1:1523210/1‑150 (MQ=255)
                    cagggTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACCGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGcc                                                                                                                                  >  2:2175483/4‑150 (MQ=255)
                        gtgagattATCTTCATCCACCACGATAGAAGCAACCTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGCACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGc                                                                                                                              <  2:1000942/144‑1 (MQ=255)
                            gtTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGGCACCGCCTGAACACGCGCGCCACGCATAc                                                                                                                          >  2:2546764/1‑150 (MQ=255)
                            gtTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGGCACCGCCTGAACACGCGCGCCACGCATAc                                                                                                                          >  2:298556/1‑150 (MQ=255)
                            gtTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATAc                                                                                                                          <  2:3945200/150‑1 (MQ=255)
                             tttATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATAcc                                                                                                                         >  1:2802561/1‑150 (MQ=255)
                                  cTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGc                                                                                                                    <  2:1323977/150‑1 (MQ=255)
                                  cTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGc                                                                                                                    >  2:3401743/1‑150 (MQ=255)
                                    tCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCaa                                                                                                                  <  1:1254820/150‑1 (MQ=255)
                                    tCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCaa                                                                                                                  <  1:661052/150‑1 (MQ=255)
                                     cATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAg                                                                                                                 >  1:2016929/1‑150 (MQ=255)
                                            cacGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTAc                                                                                                          <  2:44456/150‑1 (MQ=255)
                                                 tAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGAt                                                                                                     <  2:2992233/150‑1 (MQ=255)
                                                 tAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGAt                                                                                                     <  2:3439956/150‑1 (MQ=255)
                                                 tAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGAt                                                                                                     <  1:1675282/150‑1 (MQ=255)
                                                      gCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATa                                                                                                >  2:686756/1‑150 (MQ=255)
                                                      gCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATa                                                                                                >  2:164699/1‑150 (MQ=255)
                                                       cAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATAc                                                                                               >  1:213361/1‑150 (MQ=255)
                                                       cAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATAc                                                                                               <  1:626192/150‑1 (MQ=255)
                                                        aaCGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACg                                                                                              <  2:565297/150‑1 (MQ=255)
                                                        aaCGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACg                                                                                              >  2:52665/1‑150 (MQ=255)
                                                         aCGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGt                                                                                             <  1:823853/150‑1 (MQ=255)
                                                           gTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGGCACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGttt                                                                                           >  2:3511445/1‑150 (MQ=255)
                                                           gTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGttt                                                                                           >  2:1888869/1‑150 (MQ=255)
                                                            tCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTa                                                                                          <  1:1770465/150‑1 (MQ=255)
                                                              tGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATc                                                                                        <  1:1070575/150‑1 (MQ=255)
                                                              tGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATc                                                                                        <  1:1636545/150‑1 (MQ=255)
                                                              tGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATc                                                                                        <  1:300157/150‑1 (MQ=255)
                                                              tGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATc                                                                                        <  1:3331698/150‑1 (MQ=255)
                                                              tGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATc                                                                                        <  1:1391873/150‑1 (MQ=255)
                                                              tGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATAGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATc                                                                                        >  2:2068489/1‑150 (MQ=255)
                                                              tGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGAGATCGATACGCTCGCCACCCAGTTCAGTTGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTtctc                                                                                        <  2:1981848/150‑1 (MQ=255)
                                                                ccGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGGCACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGt                                                                                      >  2:3828954/1‑150 (MQ=255)
                                                                ccGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGt                                                                                      >  2:2808342/1‑150 (MQ=255)
                                                                  ggTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTg                                                                                    <  2:3034581/150‑1 (MQ=255)
                                                                     gCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGtt                                                                                 >  1:2890375/1‑150 (MQ=255)
                                                                      ccATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGttt                                                                                >  1:3565261/1‑150 (MQ=255)
                                                                         ttGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCa                                                                             <  1:3153418/150‑1 (MQ=255)
                                                                          tGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCAc                                                                            >  1:3099550/1‑150 (MQ=255)
                                                                             gTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcgc                                                                         >  2:3298289/1‑150 (MQ=255)
                                                                                 aTCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGAAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt                                                                     >  1:3543744/1‑150 (MQ=255)
                                                                                   cACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtttt                                                                   <  1:1015806/150‑1 (MQ=255)
                                                                                   cACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtttt                                                                   >  1:2130486/1‑150 (MQ=255)
                                                                                    aCGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTc                                                                  >  2:3969179/1‑150 (MQ=255)
                                                                                    aCGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTc                                                                  <  2:1415820/150‑1 (MQ=255)
                                                                                    aCGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTc                                                                  >  2:3627583/1‑150 (MQ=255)
                                                                                    aCGAACTGCGCCGGGTTATCATCCCACAGGACGATATAGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTc                                                                  >  2:3641847/1‑150 (MQ=255)
                                                                                      gAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGc                                                                <  1:3826911/150‑1 (MQ=255)
                                                                                          tGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACga                                                            >  1:2898620/1‑150 (MQ=255)
                                                                                              ccGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAc                                                        <  1:1907370/150‑1 (MQ=255)
                                                                                              ccGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAc                                                        <  2:213361/150‑1 (MQ=255)
                                                                                                 ggTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCg                                                     <  2:335014/150‑1 (MQ=255)
                                                                                                  gTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCgg                                                    <  2:915888/150‑1 (MQ=255)
                                                                                                  gTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCgg                                                    <  1:3960971/150‑1 (MQ=255)
                                                                                                  gTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCgg                                                    <  2:477634/150‑1 (MQ=255)
                                                                                                    tATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGAt                                                  <  2:3556745/150‑1 (MQ=255)
                                                                                                       catcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcgc                                               >  1:4016639/1‑150 (MQ=255)
                                                                                                        atcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcgcg                                              <  2:2802561/150‑1 (MQ=255)
                                                                                                         tcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGa                                             <  1:530787/150‑1 (MQ=255)
                                                                                                           ccACAGGACGCTATCGCTCCGCTCGCCACCCCGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCGACGCACGCACCTCCCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGc                                           <  2:2244752/150‑1 (MQ=255)
                                                                                                           ccACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGc                                           <  2:3708999/150‑1 (MQ=255)
                                                                                                           ccACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGc                                           >  1:1136230/1‑150 (MQ=255)
                                                                                                               aGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgct                                       >  2:4071297/1‑150 (MQ=255)
                                                                                                               aGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgct                                       >  1:2453455/1‑150 (MQ=255)
                                                                                                               aGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgct                                       >  1:362472/1‑150 (MQ=255)
                                                                                                                 gACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgctgc                                     <  2:4093045/150‑1 (MQ=255)
                                                                                                                 gACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCAATTTTCGCACGAGAACCCGGATCGCGAGCCgctgc                                     <  2:3463236/150‑1 (MQ=255)
                                                                                                                      tatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTaa                                <  1:2039630/150‑1 (MQ=255)
                                                                                                                      tatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTaa                                >  2:883231/1‑150 (MQ=255)
                                                                                                                            tACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCaa                          >  1:202613/1‑150 (MQ=255)
                                                                                                                            tACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCaa                          >  1:3059826/1‑150 (MQ=255)
                                                                                                                             aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAAt                         <  1:2844526/150‑1 (MQ=255)
                                                                                                                                cTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCAc                      >  1:2642892/1‑150 (MQ=255)
                                                                                                                                cTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCAc                      <  2:2334088/150‑1 (MQ=255)
                                                                                                                                 tCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCAct                     >  1:1433924/1‑150 (MQ=255)
                                                                                                                                 tCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCAct                     >  2:273977/1‑150 (MQ=255)
                                                                                                                                      aCCCAGTTCAGTAGACACCGCCTGAACTCGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACttcttc                <  2:247583/150‑1 (MQ=255)
                                                                                                                                      aCCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACttcttc                <  2:1791570/150‑1 (MQ=255)
                                                                                                                                      aCCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACttcttc                <  1:3524130/150‑1 (MQ=255)
                                                                                                                                        ccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGc              <  1:3758686/150‑1 (MQ=255)
                                                                                                                                        ccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGc              >  2:123191/1‑150 (MQ=255)
                                                                                                                                           gTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGa           >  1:1718158/1‑150 (MQ=255)
                                                                                                                                            ttCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTGTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAt          >  2:1863622/1‑150 (MQ=255)
                                                                                                                                                gTAGGCACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCt      <  1:3198845/150‑1 (MQ=255)
                                                                                                                                                gTAGGCACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCt      <  1:4080158/150‑1 (MQ=255)
                                                                                                                                                gTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCTAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCt      <  1:3401743/150‑1 (MQ=255)
                                                                                                                                                  aGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTgg    >  1:3688405/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGGCTTCAACGGCGATATCCATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCA  >  NC_000913/3316592‑3316889

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: