Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,533 A→T 6.3% intergenic (+32/+22) cspE → / ← flc constitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein/fluoride efflux channel, dual topology membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,5330AT 6.3% 497.3 / 15.2 142intergenic (+32/+22)cspE/flcconstitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein/fluoride efflux channel, dual topology membrane protein
Reads supporting (aligned to +/- strand):  ref base A (68/65);  new base T (7/2);  total (75/67)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 143 total aligned reads displayed.

ACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCAC  >  NC_000913/657386‑657680
                                                                                                                                                    |                                                                                                                                                   
acTTCTCTGCAATCCAGACTAATGGTTTCAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAAt                                                                                                                                                    >  1:838250/1‑149 (MQ=255)
  ttCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCa                                                                                                                                                  >  1:82061/1‑149 (MQ=255)
  ttCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCa                                                                                                                                                  >  2:265099/1‑149 (MQ=255)
  ttCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCa                                                                                                                                                  >  2:410611/1‑149 (MQ=255)
            tCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGttttt                                                                                                                                        >  2:259711/1‑149 (MQ=255)
            tCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGttttt                                                                                                                                        >  1:47640/1‑149 (MQ=255)
             ccAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAAACCCGCTGATTAAGCGGGttttg                                                                                                                                        <  1:39395/149‑2 (MQ=255)
               aGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGtttttttt                                                                                                                                     <  2:416127/149‑1 (MQ=255)
                 aCTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTgg                                                                                                                                   >  3:101778/1‑149 (MQ=255)
                    aaTGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCt                                                                                                                                >  1:777519/1‑149 (MQ=255)
                     aTGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCtt                                                                                                                               >  1:241601/1‑149 (MQ=255)
                      tGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCttt                                                                                                                              >  3:810570/1‑149 (MQ=255)
                      tGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCttt                                                                                                                              >  1:152579/1‑149 (MQ=255)
                       ggTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTa                                                                                                                             <  4:10943/149‑1 (MQ=255)
                          tttAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGG‑TTTTTTTGGTCTTTAgaga                                                                                                                         >  3:522293/1‑146 (MQ=255)
                             aaaaCTCTTGGTGAAGGTCAGCGCGTAGAGTTCAAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCACACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGc                                                                                                                       >  4:557868/1‑149 (MQ=255)
                             aaaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAAt                                                                                                                                                    <  4:694615/120‑1 (MQ=255)
                              aaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCg                                                                                                                      <  1:314508/149‑1 (MQ=255)
                               aaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCgg                                                                                                                     <  1:745364/149‑1 (MQ=255)
                                 ctctTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGtt                                                                                                                   <  2:599521/149‑1 (MQ=255)
                                  tctTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                                                                  >  2:779726/1‑149 (MQ=255)
                                      gCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                                                                              <  1:57500/149‑1 (MQ=255)
                                      gCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                                                                              <  1:648682/149‑1 (MQ=255)
                                          aaGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaa                                                                                                          >  2:815329/1‑149 (MQ=255)
                                             gTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCCCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCg                                                                                                                      <  1:74169/134‑1 (MQ=255)
                                             gTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCCCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCg                                                                                                                      >  2:74169/1‑134 (MQ=255)
                                              tCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                                                                      <  3:780962/149‑1 (MQ=255)
                                                aGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGc                                                                                                    <  3:669414/149‑1 (MQ=255)
                                                aGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGc                                                                                                    >  4:699933/1‑149 (MQ=255)
                                                 gcgcGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGcc                                                                                                   <  4:771074/149‑1 (MQ=255)
                                                   gcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                                                                                 >  1:435366/1‑149 (MQ=255)
                                                     gTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGaa                                                                                               >  3:545175/1‑149 (MQ=255)
                                                           ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                                                                         >  2:527994/1‑149 (MQ=255)
                                                            tCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGt                                                                                        <  3:223159/149‑1 (MQ=255)
                                                                    aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCAt                                                                                <  2:196500/149‑1 (MQ=255)
                                                                     cTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                                                                               <  1:707191/149‑1 (MQ=255)
                                                                     cTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                                                                               >  3:9573/1‑149 (MQ=255)
                                                                           gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaaa                                                                         <  1:674177/149‑1 (MQ=255)
                                                                              ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                                                                      >  1:232119/1‑149 (MQ=255)
                                                                              ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                                                                      <  4:896771/149‑1 (MQ=255)
                                                                                aaaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                                                                    <  2:621819/149‑1 (MQ=255)
                                                                                aaaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                                                                    >  3:17397/1‑149 (MQ=255)
                                                                                 aaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                                                                   >  3:828160/1‑149 (MQ=255)
                                                                                 aaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                                                                   <  4:101778/149‑1 (MQ=255)
                                                                                    gCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                                                                >  2:661161/1‑149 (MQ=255)
                                                                                    gCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                                                                >  3:636989/1‑149 (MQ=255)
                                                                                    gCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGa                                                                >  3:658249/1‑149 (MQ=255)
                                                                                      ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAg                                                              >  1:293306/1‑149 (MQ=255)
                                                                                       cTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAgg                                                             <  4:828160/149‑1 (MQ=255)
                                                                                         tCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGtt                                                           <  2:742849/149‑1 (MQ=255)
                                                                                           tgctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa                                                         <  4:545175/149‑1 (MQ=255)
                                                                                            gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc                                                        <  2:82061/149‑1 (MQ=255)
                                                                                            gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc                                                        <  3:706561/149‑1 (MQ=255)
                                                                                                 aaaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaaaa                                                   <  4:658249/149‑1 (MQ=255)
                                                                                                 aaaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaaaa                                                   >  4:855821/1‑149 (MQ=255)
                                                                                                  aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTa                                                                                                                             >  2:495603/1‑74 (MQ=255)
                                                                                                  aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTa                                                                                                                             <  1:495603/74‑1 (MQ=255)
                                                                                                  aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCATAAGACCCCAGAAGGTTGACGAAAAc                                                  >  3:604142/1‑149 (MQ=255)
                                                                                                  aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc                                                  <  1:661161/149‑1 (MQ=255)
                                                                                                  aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc                                                  >  4:397562/1‑149 (MQ=255)
                                                                                                   aCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACg                                                 >  1:178928/1‑149 (MQ=255)
                                                                                                     gTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCa                                                                                 <  2:708577/115‑1 (MQ=255)
                                                                                                     gTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                               <  4:17397/149‑1 (MQ=255)
                                                                                                       aaTCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCa                                             <  2:877684/149‑1 (MQ=255)
                                                                                                       aaTCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCa                                             <  4:636989/149‑1 (MQ=255)
                                                                                                        aTCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAg                                            <  2:47640/149‑1 (MQ=255)
                                                                                                           gCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTaa                                         >  1:906291/1‑149 (MQ=255)
                                                                                                              ctGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGc                                      <  3:141742/149‑1 (MQ=255)
                                                                                                              ctGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGc                                      <  1:12526/149‑1 (MQ=255)
                                                                                                                gTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGccc                                    <  2:777519/149‑1 (MQ=255)
                                                                                                                gTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGccc                                    <  2:241601/149‑1 (MQ=255)
                                                                                                                   aGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGc                                 >  2:96412/1‑149 (MQ=255)
                                                                                                                    gATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTc                                                                                  <  3:380811/99‑1 (MQ=255)
                                                                                                                    gATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTc                                                                                  >  4:380811/1‑99 (MQ=255)
                                                                                                                         gTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAACAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGc                           >  2:672969/1‑149 (MQ=255)
                                                                                                                            aGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGc                        >  3:391529/1‑149 (MQ=255)
                                                                                                                             gCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGcc                       >  2:574927/1‑149 (MQ=255)
                                                                                                                              cAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGccc                      >  3:592094/1‑149 (MQ=255)
                                                                                                                               aaGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCct                     >  4:911487/1‑149 (MQ=255)
                                                                                                                               aaGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCct                     >  4:745835/1‑149 (MQ=255)
                                                                                                                                 gAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCctct                   >  1:892177/1‑149 (MQ=255)
                                                                                                                                   aTTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTg                 <  1:728181/149‑1 (MQ=255)
                                                                                                                                    aaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCg                          <  4:171837/136‑1 (MQ=255)
                                                                                                                                    aaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCg                          >  3:171837/5‑140 (MQ=255)
                                                                                                                                      c‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTaa              <  2:178928/149‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTa               >  2:570023/3‑149 (MQ=255)
                                                                                                                                         aaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                               >  3:618748/1‑114 (MQ=255)
                                                                                                                                         aaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                               <  4:618748/114‑1 (MQ=255)
                                                                                                                                         aaaCCCGCTGATTAAGCGGGTTTTTTTTGGGCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACa            >  3:891721/1‑149 (MQ=255)
                                                                                                                                          aaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGACTGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaa           >  1:832993/1‑149 (MQ=255)
                                                                                                                                           aCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaa          >  3:287080/1‑149 (MQ=255)
                                                                                                                                           aCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaa          >  3:125216/1‑149 (MQ=255)
                                                                                                                                           aCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaa          >  4:129288/1‑149 (MQ=255)
                                                                                                                                           aCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaa          >  3:506260/1‑149 (MQ=255)
                                                                                                                                             ccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaaaa        >  3:255300/1‑149 (MQ=255)
                                                                                                                                             ccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaaaa        >  4:297887/1‑149 (MQ=255)
                                                                                                                                             ccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaaaa        <  3:297887/149‑1 (MQ=255)
                                                                                                                                                cTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAacac     <  1:672969/149‑1 (MQ=255)
                                                                                                                                                cTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAacac     >  3:691727/1‑149 (MQ=255)
                                                                                                                                                   aaTCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAcaccac  >  3:297400/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTC‑AAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCAC  >  NC_000913/657386‑657680

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑