Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,223 A→T 5.1% I54N (ATT→AAT)  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,2230AT 5.1% 400.0 / 7.2 136I54N (ATT→AAT) soxSsuperoxide response regulon transcriptional activator; autoregulator
Reads supporting (aligned to +/- strand):  ref base A (67/62);  new base T (5/2);  total (72/64)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 139 total aligned reads displayed.

ATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTC  >  NC_000913/4277075‑4277372
                                                                                                                                                    |                                                                                                                                                     
aTAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAAt                                                                                                                                                      <  2:1844798/150‑1 (MQ=255)
 tAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATg                                                                                                                                                     <  1:1135169/150‑1 (MQ=255)
  aaTCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGt                                                                                                                                                    >  1:832033/1‑150 (MQ=255)
   aTCGCTGGGAGTGCGATCAAACTGCCGACGGAACACGCGGGAGAAGGTCTGCTGCGAGCCCTAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCACCGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTa                                                                                                                                                   >  1:147687/1‑150 (MQ=255)
    tCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCACTAACAGGCGGCGTTGGCGAATGTaa                                                                                                                                                  >  1:915235/1‑150 (MQ=255)
       cTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCg                                                                                                                                               <  2:1254151/150‑1 (MQ=255)
       cTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCg                                                                                                                                               <  2:1235666/150‑1 (MQ=255)
         gggAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGcc                                                                                                                                             <  2:31703/150‑1 (MQ=255)
            cgTGCGATCAAACTGCCGACGGAAAACGCGTGAGAAGGTCTGCTGCGATACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCACCGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAg                                                                                                                                          <  1:322740/149‑1 (MQ=255)
             gTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGATTACATTCGCCaacg                                                                                                                                         <  1:1690005/150‑3 (MQ=255)
              tGCGAGCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGATTACATTCGCCaacgc                                                                                                                                        <  1:819352/150‑4 (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)
                          gCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTAAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCgacc                                                                                                                            >  1:459424/1‑147 (MQ=255)
                           ccGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCAc                                                                                                                           >  1:124616/1‑150 (MQ=255)
                             gACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCg                                                                                                                         >  1:1868219/1‑150 (MQ=255)
                                ggAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGc                                                                                                                      >  1:1744141/1‑150 (MQ=255)
                                ggAAAACGCGGGAGAAGGGCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGc                                                                                                                      <  2:1851032/150‑1 (MQ=255)
                                 gAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGAGATATCAAAAATCGTACGCTCGGTGGTGCGCAACTCAACGCCGGCCCGTAACAGGCGGCGTTCGCAAATGTAAACGCCAAGCGACTGATGCGTCCACGTGCg                                                                                                                     >  2:966858/1‑150 (MQ=255)
                                  aaaaCGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCgg                                                                                                                    >  1:1288516/1‑150 (MQ=255)
                                  aaaaCGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCgg                                                                                                                    <  2:58764/150‑1 (MQ=255)
                                       gcgGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACAt                                                                                                               <  2:1630933/150‑1 (MQ=255)
                                           gagaAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGt                                                                                                           <  1:505706/150‑1 (MQ=255)
                                                   ctgctgCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTa                                                                                                   >  1:799244/1‑150 (MQ=255)
                                                   ctgctgCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGCCGGCCAGTTACAGGCGGCGTTGGCGAATGTAATAGCCAAGCGTCTGATGCGTCAAAGTGCGGAACATTAGTTGCAAGTa                                                                                                   >  2:169354/1‑150 (MQ=255)
                                                    tgctgcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCGACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTAc                                                                                                  >  2:379726/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)
                                                        gcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACt                                                                                              >  1:1031890/1‑150 (MQ=255)
                                                         cGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACtt                                                                                             >  2:1381941/1‑150 (MQ=255)
                                                            gaCATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTAcc                                                                                                 >  2:356844/1‑143 (MQ=255)
                                                            gaCATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGa                                                                                          >  1:1753019/1‑150 (MQ=255)
                                                              cATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAAt                                                                                        <  2:1516056/150‑1 (MQ=255)
                                                              cATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAAt                                                                                        >  2:1898951/1‑150 (MQ=255)
                                                                 aaCCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGc                                                                                     <  1:130416/150‑1 (MQ=255)
                                                                   cccAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCt                                                                                   <  2:432192/150‑1 (MQ=255)
                                                                   cccAGGTCCATTGCGATATCAAAAATCGGACCCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAAAAGGCGGCGTTCGCGAATGTACTCGCCAAGCGTCTGAGGCGTCACCGAGCGGAACAATCGTTGCAAGTACCACTTTGAATAGCCt                                                                                   >  1:1547551/1‑150 (MQ=255)
                                                                       ggTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAtt                                                                               >  2:1732777/1‑150 (MQ=255)
                                                                        gTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGCTGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAttt                                                                              >  1:306000/1‑150 (MQ=255)
                                                                        gTCCAATGCGATCTCAAACATCGGACGCGCGGGGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGGCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAttt                                                                              <  2:909697/150‑1 (MQ=255)
                                                                             ttGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTg                                                                         <  1:418175/150‑1 (MQ=255)
                                                                             ttGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCGGATGCGTCACCGTGCGGACCATTCGTTGCACGTACCACTTTGAATAGCCTGATTTTTTTg                                                                         >  1:163900/1‑150 (MQ=255)
                                                                              tGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGc                                                                        <  1:973789/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)
                                                                                      cAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACAt                                                                >  2:1316973/1‑150 (MQ=255)
                                                                                       aaaaaTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGTCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGATTCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATc                                                               <  1:1414278/150‑1 (MQ=255)
                                                                                       aaaaaTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACAta                                                               >  1:709692/1‑149 (MQ=255)
                                                                                             cGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGtt                                                         >  2:1647713/1‑150 (MQ=255)
                                                                                              ggACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAGCATTCGTTGCAAGTACCACTTTGAATAGCCTGATGTTTTTGCGAGTACATCATGGTTa                                                        >  2:1156963/1‑150 (MQ=255)
                                                                                              ggACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTa                                                        >  1:41763/1‑150 (MQ=255)
                                                                                                 cGCTCGGTGGTGCGCATCTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCTTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGc                                                     <  2:1361941/150‑1 (MQ=255)
                                                                                                  gCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCg                                                    >  1:1791481/1‑150 (MQ=255)
                                                                                                   cTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCgg                                                   >  2:673200/1‑150 (MQ=255)
                                                                                                         ggtgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTc                                             <  1:1642249/150‑1 (MQ=255)
                                                                                                          gtgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCa                                            <  1:504295/150‑1 (MQ=255)
                                                                                                             ctctaCTCAACGGCGGCCAGTTAGAGGCTGCGTTGGCGAATGTATTCGCCAAGCGTCTGATGCGTCACCGTTCTGAACATTCGTTGCAAGTACCACTTTGAATATCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAata                                         <  1:276738/146‑1 (MQ=255)
                                                                                                             cgcTACTCAACGGCGGCCTGTAACAGGAGGCGTTGGCGAATTTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAAGATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACCTCAATGTTAAGCGGCTGGTCAata                                         <  2:52427/150‑1 (MQ=255)
                                                                                                              gcAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAatat                                        <  1:71153/150‑1 (MQ=255)
                                                                                                                 aCTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCt                                     <  2:1760039/150‑1 (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)
                                                                                                                   tCAACGGCGGCCAGTAACAGGCGGCGTTGGAGAATGTAATCGCCAAGAGTATGATGAGTCACCGTGCGGAACATTAGTTGCAAGTTCGACTTTGAATATCCTTATTTTTTTGCGACTTCATCAATGTTAAGCGGCTGGGCAATATGCTCg                                   >  2:270086/1‑150 (MQ=255)
                                                                                                                    cAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGt                                  >  1:1378558/1‑150 (MQ=255)
                                                                                                                     aaCGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTc                                 <  2:486612/150‑1 (MQ=255)
                                                                                                                     aaCGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTc                                 <  2:1562931/150‑1 (MQ=255)
                                                                                                                     aaCGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTc                                 <  1:332862/150‑1 (MQ=255)
                                                                                                                     aaCGGCGGCCAGTAACAGGCGGCGGTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCAGCGTGCGGAACATTCGTTGCAAGTCCCACTTAGAATAGCCTGATTTTTTTGCGACTACATCATTGTTAAGCGGCTGGTCAATATGCTCGTc                                 <  2:191793/150‑1 (MQ=255)
                                                                                                                       cggcggCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCaa                               <  2:1906736/150‑1 (MQ=255)
                                                                                                                           ggcCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCa                           <  1:1105720/150‑1 (MQ=255)
                                                                                                                            gcCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTAGTCAATCCAt                          >  1:1174176/1‑150 (MQ=255)
                                                                                                                             ccAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATg                         >  2:1247477/1‑150 (MQ=255)
                                                                                                                             ccAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATg                         <  2:1744141/150‑1 (MQ=255)
                                                                                                                               aGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCg                       <  2:1379896/150‑1 (MQ=255)
                                                                                                                                 tAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGAt                     <  2:1302220/150‑1 (MQ=255)
                                                                                                                                  aaCAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATa                    >  2:54448/1‑150 (MQ=255)
                                                                                                                                    cAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAg                  >  2:1311958/1‑150 (MQ=255)
                                                                                                                                      ggcggcGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGAt                >  1:658346/1‑150 (MQ=255)
                                                                                                                                      ggcggcGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGAt                >  2:1617025/1‑150 (MQ=255)
                                                                                                                                         ggcgTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCAACGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTAGTCAATCCATGCGATAAGATCCt             >  1:142357/1‑150 (MQ=255)
                                                                                                                                         acgcttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCt             >  1:1716493/5‑150 (MQ=255)
                                                                                                                                            cttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaa          >  1:977047/2‑150 (MQ=255)
                                                                                                                                            cttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaa          >  1:1534353/2‑150 (MQ=255)
                                                                                                                                              tGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaata        >  2:739093/1‑150 (MQ=255)
                                                                                                                                               ggCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaataa       >  2:864943/1‑150 (MQ=255)
                                                                                                                                                 cGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGCCTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAttt     >  1:1662056/1‑150 (MQ=255)
                                                                                                                                                   aaTGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAAtttt   <  2:799244/150‑1 (MQ=255)
                                                                                                                                                    aTGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTc  <  1:1111568/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTC  >  NC_000913/4277075‑4277372

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