Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,223 A→T 6.3% 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 6.3% 319.5 / 8.9 111I54N (ATT→AAT) soxSsuperoxide response regulon transcriptional activator; autoregulator
Reads supporting (aligned to +/- strand):  ref base A (52/52);  new base T (2/5);  total (54/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 = 4.97e-01

Only 100 of 112 total aligned reads displayed.

GATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAA  >  NC_000913/4277074‑4277367
                                                                                                                                                     |                                                                                                                                                
  aaaTCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATg                                                                                                                                                <  1:732341/149‑1 (MQ=255)
    aTCGCTGGGAGTGCGATCAACCTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACAAAACCCAGGTCCATTGCGATATCAAAAATCGGCCGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACCGGCGGCGTTGGCGAATGTa                                                                                                                                              >  1:84785/1‑150 (MQ=255)
      cGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAAt                                                                                                                                            <  1:1279933/150‑1 (MQ=255)
       gCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATc                                                                                                                                           >  1:1675740/1‑150 (MQ=255)
            gAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCATTGGCGATTACATTCGCCaa                                                                                                                                      <  2:1479919/150‑1 (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)
                  gATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGC‑TTGGCGATTACATTCGCCaacgccgcc                                                                                                                               <  2:1600338/150‑8 (MQ=255)
                       aaCTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGc                                                                                                                           >  1:537554/1‑150 (MQ=255)
                       aaCTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGc                                                                                                                           <  2:1675740/150‑1 (MQ=255)
                         cTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGt                                                                                                                         >  1:433686/1‑150 (MQ=255)
                           gCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCa                                                                                                                       <  1:623031/150‑1 (MQ=255)
                           gCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCa                                                                                                                       >  2:778986/1‑150 (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)
                                                      gctgcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTAcc                                                                                            <  2:1580560/150‑1 (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)
                                                                 tAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCATTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAg                                                                                                <  1:139958/135‑1 (MQ=255)
                                                                   aCCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGcc                                                                               >  1:855757/1‑150 (MQ=255)
                                                                   aCCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGcc                                                                               <  2:662203/150‑1 (MQ=255)
                                                                    cccAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCt                                                                              <  1:543936/150‑1 (MQ=255)
                                                                     ccAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTACCAGGCGGCTTAGGCGAAAGTAAACGCCAAGCGACTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATACCCTg                                                                             >  2:319232/1‑150 (MQ=255)
                                                                     ccAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTg                                                                             <  2:218069/150‑1 (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)
                                                                               tGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGc                                                                   <  1:202041/150‑1 (MQ=255)
                                                                               tGCGATATCAAAAATCGGACGCGCGGCGGTGCGCAACTCAACCGCGGCCAGTAACAGGCGGCGGTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGc                                                                   <  2:1252773/150‑1 (MQ=255)
                                                                                  gATATCAAAAATCGGCCGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTCCCACGTTGAATAGCCTGCTTTTTTTGCGAc                                                                >  2:829205/1‑150 (MQ=255)
                                                                                  gATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGAc                                                                >  2:451415/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)
                                                                                             tCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGt                                                     <  1:762370/150‑1 (MQ=255)
                                                                                                  cGCTCGGTGGTGCGCAACTCAACGGCGTCCAGTAACAGGCTGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGA‑TTTTTTGCGACTACATCAATGTTAAGCg                                               <  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:87122/150‑1 (MQ=255)
                                                                                                                cAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATg                                  <  2:1644481/150‑1 (MQ=255)
                                                                                                                 aaCTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGc                                 >  2:485446/1‑150 (MQ=255)
                                                                                                                   cTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTc                               >  1:890650/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)
                                                                                                                                 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)
                                                                                                                                           gcgTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTg       <  2:1566728/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                
GATAATCGCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAA  >  NC_000913/4277074‑4277367

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: ‑