Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,228 A→T 5.6% D52E (GAT→GAA soxS ← superoxide response regulon transcriptional activator; autoregulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,277,2280AT 5.6% 409.6 / 7.9 134D52E (GAT→GAAsoxSsuperoxide response regulon transcriptional activator; autoregulator
Reads supporting (aligned to +/- strand):  ref base A (64/59);  new base T (5/3);  total (72/62)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

TCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTGATG  >  NC_000913/4277079‑4277377
                                                                                                                                                     |                                                                                                                                                     
tCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCACTAACAGGCGGCGTTGGCGAATGTaa                                                                                                                                                       >  1:915235/1‑150 (MQ=255)
   cTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCg                                                                                                                                                    <  2:1235666/150‑1 (MQ=255)
     gggAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGcc                                                                                                                                                  <  2:31703/150‑1 (MQ=255)
        cgTGCGATCAAACTGCCGACGGAAAACGCGTGAGAAGGTCTGCTGCGATACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCACCGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAg                                                                                                                                               <  1:322740/149‑1 (MQ=255)
        aGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAg                                                                                                                                               <  2:33559/150‑1 (MQ=255)
         gTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGATTACATTCGCCaacg                                                                                                                                              <  1:1690005/150‑3 (MQ=255)
           gCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCACCTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGctt                                                                                                                                            >  1:758563/1‑148 (MQ=255)
           gCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCAATGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGt                                                                                                                                            >  1:1302220/1‑150 (MQ=255)
             gATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCt                                                                                                                                          >  2:973789/1‑150 (MQ=255)
                 aaaCTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATg                                                                                                                                      <  1:367143/150‑1 (MQ=255)
                      gCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGATCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCa                                                                                                                                 <  2:138474/150‑1 (MQ=255)
                       ccGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCAc                                                                                                                                >  1:124616/1‑150 (MQ=255)
                            ggAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTTGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGc                                                                                                                           <  1:1249736/150‑1 (MQ=255)
                            ggAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGc                                                                                                                           >  1:1744141/1‑150 (MQ=255)
                            ggAAAACGCGGGAGAAGGGCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGc                                                                                                                           <  2:1851032/150‑1 (MQ=255)
                             gAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGAGATATCAAAAATCGTACGCTCGGTGGTGCGCAACTCAACGCCGGCCCGTAACAGGCGGCGTTCGCAAATGTAAACGCCAAGCGACTGATGCGTCCACGTGCg                                                                                                                          >  2:966858/1‑150 (MQ=255)
                              aaaaCGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCgg                                                                                                                         <  2:58764/150‑1 (MQ=255)
                                    cgGGAGAAGGTCTGCTGCGATACATAACCAAGGTCCATTGCGATATCAAAAATCGCACGATCGGTGGTGCCCAACTCAACGGCCAACAGTAACAGTCGGCCTTGACGAATGTAATCCCCAAGCGTCTGAGTCGTCACCGTGCCGAACAtt                                                                                                                   >  2:681172/1‑150 (MQ=255)
                                       gagaAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGt                                                                                                                <  1:505706/150‑1 (MQ=255)
                                         gaAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTg                                                                                                              <  1:1378392/150‑1 (MQ=255)
                                            ggTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGGGGACCATTCGTTGCaa                                                                                                           >  1:1163845/1‑150 (MQ=255)
                                               ctgctgCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTa                                                                                                        >  1:799244/1‑150 (MQ=255)
                                               ctgctgCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGCCGGCCAGTTACAGGCGGCGTTGGCGAATGTAATAGCCAAGCGTCTGATGCGTCAAAGTGCGGAACATTAGTTGCAAGTa                                                                                                        >  2:169354/1‑150 (MQ=255)
                                                tgctgcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTAc                                                                                                       >  2:42416/1‑150 (MQ=255)
                                                  ctgcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCa                                                                                                     <  2:268255/150‑1 (MQ=255)
                                                   tgcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCAc                                                                                                    >  1:1145437/1‑150 (MQ=255)
                                                    gcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCACCGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCAGATGCGTCACCGTACGGAACATTCGTAGCAAGTACCACt                                                                                                   >  1:293580/1‑150 (MQ=255)
                                                    gcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACt                                                                                                   >  1:1031890/1‑150 (MQ=255)
                                                     cGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACtt                                                                                                  >  2:1381941/1‑150 (MQ=255)
                                                        gaCATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTAcc                                                                                                      >  2:356844/1‑143 (MQ=255)
                                                        gaCATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTAcc                                                                                                      <  1:356844/143‑1 (MQ=255)
                                                          cATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAAt                                                                                             >  2:1898951/1‑150 (MQ=255)
                                                            tAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCACCGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAg                                                                                           >  1:1494830/1‑150 (MQ=255)
                                                            aaaCCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAg                                                                                           <  1:825056/149‑1 (MQ=255)
                                                             aaCCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGc                                                                                          <  1:130416/150‑1 (MQ=255)
                                                               cccAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCt                                                                                        <  2:432192/150‑1 (MQ=255)
                                                               cccAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCt                                                                                        <  1:1764765/150‑1 (MQ=255)
                                                                ccAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTg                                                                                       <  2:1340415/150‑1 (MQ=255)
                                                                   ggTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAtt                                                                                    >  2:1732777/1‑150 (MQ=255)
                                                                    gTCCATTGCGATATCAAAACTCGGCCGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAttt                                                                                   <  2:1553405/150‑1 (MQ=255)
                                                                    gTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGCTGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAttt                                                                                   >  1:306000/1‑150 (MQ=255)
                                                                    gTCCAATGCGATCTCAAACATCGGACGCGCGGGGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGGCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAttt                                                                                   <  2:909697/150‑1 (MQ=255)
                                                                         ttGCGATATCAAAATTCGGACTCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACATGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTAGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTg                                                                              <  1:1407920/150‑1 (MQ=255)
                                                                         ttGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTg                                                                              <  1:418175/150‑1 (MQ=255)
                                                                         ttGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCGGATGCGTCACCGTGCGGACCATTCGTTGCACGTACCACTTTGAATAGCCTGATTTTTTTg                                                                              >  1:163900/1‑150 (MQ=255)
                                                                          tGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGc                                                                             <  1:794931/150‑1 (MQ=255)
                                                                           gCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCg                                                                            <  2:589455/150‑1 (MQ=255)
                                                                             gATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGAc                                                                          <  1:1696282/150‑1 (MQ=255)
                                                                              atatCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACACGCGGCGTTGGCGAATGTAAACGCCAAGCGTCTGATGCGTCACCGTGCGGCACATTCCTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACt                                                                         >  1:206946/1‑150 (MQ=255)
                                                                               tatCAAAAATCGGACGCTCGGGGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTa                                                                        <  1:404743/150‑1 (MQ=255)
                                                                                 acAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTACTCGCCACGCGTCTGATGCGTCAACGTGCGGAACATTCGTTGCGAGTACCACTTTGAATAGCCTGAAGTTTTTGCGCCTAcc                                                                      >  1:1147498/2‑149 (MQ=255)
                                                                                  cAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTACTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGCCTACAt                                                                     >  2:743732/1‑150 (MQ=255)
                                                                                  cAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACaa                                                                     >  2:1505322/1‑149 (MQ=255)
                                                                                  cAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACAt                                                                     >  1:578162/1‑150 (MQ=255)
                                                                                  cAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACAt                                                                     >  1:1773884/1‑150 (MQ=255)
                                                                                   aaaaaTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACAta                                                                    >  1:709692/1‑149 (MQ=255)
                                                                                    aaaaTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCa                                                                   >  1:927348/1‑150 (MQ=255)
                                                                                    aaaaTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCa                                                                   <  2:19654/150‑1 (MQ=255)
                                                                                         cGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGtt                                                              >  2:1647713/1‑150 (MQ=255)
                                                                                          ggACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTa                                                             >  1:41763/1‑150 (MQ=255)
                                                                                          ggACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTa                                                             >  1:1676978/1‑150 (MQ=255)
                                                                                             cGCTCGGTGGTGCGCTACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGc                                                          <  2:1837244/150‑1 (MQ=255)
                                                                                             cGCTCGGTGGTGCGCATCTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCTTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGc                                                          <  2:1361941/150‑1 (MQ=255)
                                                                                              gCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCg                                                         >  1:1791481/1‑150 (MQ=255)
                                                                                                     ggtgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTc                                                  <  1:1642249/150‑1 (MQ=255)
                                                                                                      gtgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCa                                                 <  1:504295/150‑1 (MQ=255)
                                                                                                         cgcTACTCAACGGCGGCCTGTAACAGGAGGCGTTGGCGAATTTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAAGATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACCTCAATGTTAAGCGGCTGGTCAata                                              <  2:52427/150‑1 (MQ=255)
                                                                                                          gcAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAatat                                             >  2:1644727/1‑150 (MQ=255)
                                                                                                             aCTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCt                                          >  1:1871194/1‑150 (MQ=255)
                                                                                                             aCTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCt                                          <  1:624651/150‑1 (MQ=255)
                                                                                                             aCTCAACGGCGGCCAGTAACAGGCGGCGGTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCt                                          <  1:828868/150‑1 (MQ=255)
                                                                                                               tCAACGGCGGCCAGTAACCGGCGGCGTTGGCGAATGTAAGCGCCAAGCGTCTGATGCGGCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCg                                        <  2:1144522/150‑1 (MQ=255)
                                                                                                               tCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCg                                        >  1:1379999/1‑150 (MQ=255)
                                                                                                                cAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGt                                       >  1:1734098/1‑150 (MQ=255)
                                                                                                                 aaCGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTc                                      <  2:1562931/150‑1 (MQ=255)
                                                                                                                 aaCGGCGGCCAGTAACAGGCGGCGGTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCAGCGTGCGGAACATTCGTTGCAAGTCCCACTTAGAATAGCCTGATTTTTTTGCGACTACATCATTGTTAAGCGGCTGGTCAATATGCTCGTc                                      <  2:191793/150‑1 (MQ=255)
                                                                                                                       ggcCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCa                                <  1:1105720/150‑1 (MQ=255)
                                                                                                                         ccAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATg                              <  2:1744141/150‑1 (MQ=255)
                                                                                                                         ccAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATg                              >  2:1247477/1‑150 (MQ=255)
                                                                                                                            gTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGCACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGa                           >  1:1601997/1‑150 (MQ=255)
                                                                                                                            gTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGa                           >  2:1294528/1‑150 (MQ=255)
                                                                                                                             tAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGAt                          <  2:1302220/150‑1 (MQ=255)
                                                                                                                             tAACAGGCGGCATTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATGCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGAt                          <  2:103498/150‑1 (MQ=255)
                                                                                                                                cAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAg                       >  2:1311958/1‑150 (MQ=255)
                                                                                                                                  ggcggcGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGAt                     >  2:1617025/1‑150 (MQ=255)
                                                                                                                                  ggcggcGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGAt                     >  1:658346/1‑150 (MQ=255)
                                                                                                                                    cggcgTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATcc                   <  2:578162/150‑1 (MQ=255)
                                                                                                                                     ggcgTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCAACGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTAGTCAATCCATGCGATAAGATCCt                  >  1:142357/1‑150 (MQ=255)
                                                                                                                                       gcttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGa                >  2:416488/3‑150 (MQ=255)
                                                                                                                                       cgTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGa                >  2:497554/1‑150 (MQ=255)
                                                                                                                                        cttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaa               >  1:1534353/2‑150 (MQ=255)
                                                                                                                                        cttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaa               >  1:977047/2‑150 (MQ=255)
                                                                                                                                         ttGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaat              <  2:124616/150‑1 (MQ=255)
                                                                                                                                          tGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaata             >  2:739093/1‑150 (MQ=255)
                                                                                                                                           ggCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaataa            >  2:864943/1‑150 (MQ=255)
                                                                                                                                               aaTGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAAtttt        <  2:799244/150‑1 (MQ=255)
                                                                                                                                                aTGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTc       <  1:1111568/150‑1 (MQ=255)
                                                                                                                                                   tAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTGa    <  1:1898951/150‑1 (MQ=255)
                                                                                                                                                   tAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAAGCCATGCGATAAGATCCTGAATAATTTTCTGa    >  1:947869/1‑150 (MQ=255)
                                                                                                                                                     aTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTGATg  >  1:1235369/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTGATG  >  NC_000913/4277079‑4277377

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