Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,228 A→T 5.8% 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.8% 319.6 / 4.5 112D52E (GAT→GAAsoxSsuperoxide response regulon transcriptional activator; autoregulator
Reads supporting (aligned to +/- strand):  ref base A (52/52);  new base T (2/5);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

CGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTG  >  NC_000913/4277080‑4277374
                                                                                                                                                    |                                                                                                                                                  
cGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAAt                                                                                                                                                   <  1:1279933/150‑1 (MQ=255)
 gCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATc                                                                                                                                                  >  1:1675740/1‑150 (MQ=255)
      gAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCATTGGCGATTACATTCGCCaa                                                                                                                                             <  2:1479919/150‑1 (MQ=255)
         tGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGATTACATTCGCCaacgc                                                                                                                                          <  1:1476165/150‑4 (MQ=255)
          gCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCCGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAGCAGGCGGCGTTGGCGAAAGTAATCGCCAAGCgg                                                                                                                                         >  1:826407/1‑149 (MQ=255)
           cGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAACTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTc                                                                                                                                        >  1:87122/1‑150 (MQ=255)
           cGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTc                                                                                                                                        >  1:977907/1‑150 (MQ=255)
           cGATCAAACTGCCGACGGAAAACGCGGGAGAAGGGCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGATTACATTCGCCaacgccg                                                                                                                                        <  2:582440/150‑6 (MQ=255)
            gATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGCCATAACCCAGGTCCATTGCGATATCAAAAATCGGACGATCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCt                                                                                                                                       >  2:429698/1‑150 (MQ=255)
            gATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGATTACATTCGCCaacgccgc                                                                                                                                       <  2:791757/150‑7 (MQ=255)
            gATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTca                                                                                                                                       >  1:851278/1‑149 (MQ=255)
            gATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCt                                                                                                                                       >  1:361271/1‑150 (MQ=255)
            gATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCTTGGCGATTACATTCGCCaacgccgcc                                                                                                                                      <  2:1600338/150‑8 (MQ=255)
                 aaCTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGc                                                                                                                                  >  1:537554/1‑150 (MQ=255)
                 aaCTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGc                                                                                                                                  <  2:1675740/150‑1 (MQ=255)
                     gCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCa                                                                                                                              >  2:778986/1‑150 (MQ=255)
                     gCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCa                                                                                                                              <  1:623031/150‑1 (MQ=255)
                     gCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGGCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCTACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTACTCGCCAAGCGTCTGATGCGTcc                                                                                                                              >  1:116833/1‑149 (MQ=255)
                         aCGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGt                                                                                                                          >  1:667449/1‑150 (MQ=255)
                          cGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTg                                                                                                                         >  2:1099973/1‑150 (MQ=255)
                          cGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTg                                                                                                                         >  2:499491/1‑150 (MQ=255)
                          cGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGCACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTg                                                                                                                         <  2:293263/150‑1 (MQ=255)
                           ggAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGc                                                                                                                        >  1:705283/1‑150 (MQ=255)
                             aaaaCGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCgg                                                                                                                      <  2:1298446/150‑1 (MQ=255)
                              aaaCGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGa                                                                                                                     <  1:1631788/150‑1 (MQ=255)
                                aCGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAAc                                                                                                                   <  2:1104448/150‑1 (MQ=255)
                                aCGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAAc                                                                                                                   <  2:406320/150‑1 (MQ=255)
                                aCGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGGGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAAc                                                                                                                   <  1:894019/150‑1 (MQ=255)
                                 cgcgGGAGAAGGGCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGCCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACa                                                                                                                  <  2:1686583/150‑1 (MQ=255)
                                   cgGGAGAAGGTCTGCTGCGTGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGACTGATGCGTCACCGTGCGGAACAtt                                                                                                                >  2:1382233/1‑150 (MQ=255)
                                      gagaAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGt                                                                                                             >  1:536311/1‑150 (MQ=255)
                                        gaAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTg                                                                                                           <  2:309720/150‑1 (MQ=255)
                                         aaGGTCTGCTGCGAGATATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGACTGATCCGTCACCGTGCGGAACATTCGTTGc                                                                                                          >  2:516864/1‑150 (MQ=255)
                                         aaGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGc                                                                                                          <  2:194489/150‑1 (MQ=255)
                                          aGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCa                                                                                                         <  1:1016307/150‑1 (MQ=255)
                                              ctgctgCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTa                                                                                                     >  1:1649719/1‑150 (MQ=255)
                                              ctgctgCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACAATCGTTGCAAGTa                                                                                                     >  1:876390/1‑150 (MQ=255)
                                                gctgcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTAcc                                                                                                   <  2:1580560/150‑1 (MQ=255)
                                                    cGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACtt                                                                                               >  1:738050/1‑150 (MQ=255)
                                                      agaCATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCGCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGGGCGGAACATTCGTTGCAAGTACCACTTTg                                                                                             <  2:1354435/150‑1 (MQ=255)
                                                       gaCATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGAAGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGa                                                                                            <  1:1674756/150‑1 (MQ=255)
                                                        aCATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGaa                                                                                           <  1:778986/150‑1 (MQ=255)
                                                         cATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAAt                                                                                          >  2:1314706/1‑150 (MQ=255)
                                                           tAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCATTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAg                                                                                                       >  2:139958/1‑135 (MQ=255)
                                                             aCCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGcc                                                                                      <  2:662203/150‑1 (MQ=255)
                                                             aCCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGcc                                                                                      >  1:855757/1‑150 (MQ=255)
                                                              cccAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCt                                                                                     <  1:543936/150‑1 (MQ=255)
                                                               ccAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTACCAGGCGGCTTAGGCGAAAGTAAACGCCAAGCGACTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATACCCTg                                                                                    >  2:319232/1‑150 (MQ=255)
                                                               ccAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTg                                                                                    <  2:218069/150‑1 (MQ=255)
                                                                   gTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAttt                                                                                >  2:1158992/1‑150 (MQ=255)
                                                                    tCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAtttt                                                                               <  1:1279998/150‑1 (MQ=255)
                                                                      cATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAtttttt                                                                             >  1:518777/1‑150 (MQ=255)
                                                                      cATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAtttttt                                                                             <  2:918452/150‑1 (MQ=255)
                                                                       aTTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAttttttt                                                                            <  1:1158992/150‑1 (MQ=255)
                                                                            gATATCAAAAATCGGCCGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTCCCACGTTGAATAGCCTGCTTTTTTTGCGAc                                                                       >  2:829205/1‑150 (MQ=255)
                                                                            gATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGAc                                                                       >  2:451415/1‑150 (MQ=255)
                                                                             atatCAAAAATCGGACGCTCGGTGGTGCGCACCTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACt                                                                      >  1:250674/1‑150 (MQ=255)
                                                                             atatCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACt                                                                      >  2:360682/1‑150 (MQ=255)
                                                                                 cAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTCCCCCTTTGAATAGCCTGATTTTTTTGCGCCTACAt                                                                  >  2:146849/1‑150 (MQ=255)
                                                                                  aaaaaTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATc                                                                 <  1:428774/150‑1 (MQ=255)
                                                                                      aTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACGTGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATg                                                             >  2:621714/1‑150 (MQ=255)
                                                                                            cGCTCGGTGGTGCGCAACTCAACGGCGTCCAGTAACAGGCTGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTGCGACTACATCAATGTTAAGCg                                                      <  2:826407/150‑1 (MQ=255)
                                                                                            cGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGc                                                       <  2:1649719/150‑1 (MQ=255)
                                                                                             gCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCg                                                      <  1:941105/150‑1 (MQ=255)
                                                                                                  gtggtgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTgg                                                 >  1:1641817/1‑150 (MQ=255)
                                                                                                    ggtgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTc                                               <  2:433686/150‑1 (MQ=255)
                                                                                                    ggggcgcAACTCAACGGCGGCCAGTAACAGGAGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCAACGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTc                                               <  1:1054232/147‑1 (MQ=255)
                                                                                                     gtgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCa                                              <  2:537554/150‑1 (MQ=255)
                                                                                                      tgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCaa                                             >  1:284856/1‑150 (MQ=255)
                                                                                                      tgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCaa                                             >  1:342370/1‑150 (MQ=255)
                                                                                                        cgcAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAata                                           <  2:536311/150‑1 (MQ=255)
                                                                                                         gcAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAatat                                          <  2:977907/150‑1 (MQ=255)
                                                                                                          caccTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTTATCGCCAAGCGTCTGATGCGACACCGTGCAGAACCTGCGTTGCAAGGACCAAATTGAATAGCCTTCATATTTTGCTACTACAAGAATATTTAGCGTATTGTCAATATg                                         >  2:56814/4‑150 (MQ=255)
                                                                                                          cAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATg                                         <  2:1644481/150‑1 (MQ=255)
                                                                                                          cAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATg                                         <  2:87122/150‑1 (MQ=255)
                                                                                                           aaCTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGc                                        >  2:485446/1‑150 (MQ=255)
                                                                                                             cTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTc                                      >  2:474468/1‑150 (MQ=255)
                                                                                                              tgaaCGGCGGCCAGTAACAGGCGGCGTTTGCGTATGTAATCGTCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCg                                     <  1:326418/148‑1 (MQ=255)
                                                                                                              tCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCg                                     <  1:1314706/150‑1 (MQ=255)
                                                                                                              acAACTGCGGCCAGTAACAGTCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCg                                     <  2:1642837/149‑1 (MQ=255)
                                                                                                                 aCGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCa                                  <  2:667449/150‑1 (MQ=255)
                                                                                                                        ccAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATg                           >  1:1200180/1‑150 (MQ=255)
                                                                                                                         cAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGc                          >  1:1651226/1‑150 (MQ=255)
                                                                                                                          aGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCg                         >  2:1488568/1‑150 (MQ=255)
                                                                                                                           gTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCCTGCGa                        >  2:1055804/1‑150 (MQ=255)
                                                                                                                           gTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGa                        >  1:668199/1‑150 (MQ=255)
                                                                                                                           gTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCAGGCGa                        >  2:245056/1‑150 (MQ=255)
                                                                                                                             aaCAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATa                      <  1:360682/150‑1 (MQ=255)
                                                                                                                              aCAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATaa                     <  1:1488568/150‑1 (MQ=255)
                                                                                                                                  gcggcgTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATc                 <  2:855757/150‑1 (MQ=255)
                                                                                                                                  gcggcgTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGATCGTCAATCCATGCGATAAGATc                 >  1:775653/1‑150 (MQ=255)
                                                                                                                                       gTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaa            <  2:1482440/150‑1 (MQ=255)
                                                                                                                                       cttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGATTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaa            >  2:339682/2‑150 (MQ=255)
                                                                                                                                        ttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaat           >  1:87776/1‑150 (MQ=255)
                                                                                                                                         tGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaata          <  1:405563/150‑1 (MQ=255)
                                                                                                                                          ggCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATACTCaataa         >  2:44821/1‑150 (MQ=255)
                                                                                                                                          ggCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTACGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaataa         <  2:851278/150‑1 (MQ=255)
                                                                                                                                                tGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCt   <  2:1651226/150‑1 (MQ=255)
                                                                                                                                                 gTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTTGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTg  >  2:1049092/1‑150 (MQ=255)
                                                                                                                                                 gTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTg  >  1:1233967/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTG  >  NC_000913/4277080‑4277374

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