Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,535 C→A 9.3% intergenic (+34/+20) 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,5350CA 9.3% 547.3 / 33.3 161intergenic (+34/+20)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 C (75/71);  new base A (10/5);  total (85/76)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 163 total aligned reads displayed.

TTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTT  >  NC_000913/657388‑657682
                                                                                                                                                   |                                                                                                                                                   
ttCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCa                                                                                                                                                    >  4:860912/1‑149 (MQ=255)
  ctctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATCCAAAACCCGCTTAATCAGc                                                                                                                                                  >  3:258371/1‑149 (MQ=255)
   tctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                                                                 <  1:1009267/149‑1 (MQ=255)
   tctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                                                                 >  3:573935/1‑149 (MQ=255)
    ctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCgg                                                                                                                                                <  4:918646/149‑1 (MQ=255)
        aaTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttt                                                                                                                                            <  1:1241708/149‑1 (MQ=255)
        aaTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttt                                                                                                                                            <  1:175276/149‑1 (MQ=255)
          tCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttt                                                                                                                                          >  4:783180/1‑149 (MQ=255)
             aGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                                                                       <  4:261500/149‑1 (MQ=255)
              gACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                                                                      >  2:1129988/1‑149 (MQ=255)
                 tAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                                                                       >  2:1245755/1‑145 (MQ=255)
                  aaTGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCt                                                                                                                                  >  3:151964/1‑149 (MQ=255)
                       ttttAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAgt                                                                                                                             >  4:1041490/1‑149 (MQ=255)
                         ttAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAgtgt                                                                                                                           <  4:295271/149‑1 (MQ=255)
                           aaaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGc                                                                                                                         <  2:1192138/149‑1 (MQ=255)
                            aaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCg                                                                                                                        <  3:860912/149‑1 (MQ=255)
                             aaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                                                                       <  3:957910/133‑1 (MQ=255)
                             aaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                                                                       >  4:957910/1‑133 (MQ=255)
                                tctTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                                                                    >  1:934230/1‑149 (MQ=255)
                                tctTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                                                                    >  2:1627/1‑149 (MQ=255)
                                tctTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTg                                                                                                                    >  2:666023/1‑149 (MQ=255)
                                   tGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAgg                                                                                                                 <  1:320748/149‑1 (MQ=255)
                                    gCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAATAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                                                                                <  2:1407405/149‑1 (MQ=255)
                                    gCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                                                                                <  1:637007/149‑1 (MQ=255)
                                    gCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                                                                                <  4:1102749/149‑1 (MQ=255)
                                       gAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGa                                                                                                             >  2:398731/1‑149 (MQ=255)
                                          ggTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaaa                                                                                                          <  3:117301/149‑1 (MQ=255)
                                            tCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                                                                        >  4:376560/1‑149 (MQ=255)
                                              aGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGc                                                                                                      <  4:214321/149‑1 (MQ=255)
                                              aGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGc                                                                                                      <  1:666023/149‑1 (MQ=255)
                                               gcgcGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGcc                                                                                                     >  1:704875/1‑149 (MQ=255)
                                                cgcgCAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTACTGCTGCAATCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCa                                                                                                    <  4:1229468/149‑1 (MQ=255)
                                                    tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAt                                                                                                <  4:421392/149‑1 (MQ=255)
                                                        gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCa                                                                                            <  3:837900/149‑1 (MQ=255)
                                                         ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                                                                        >  1:569516/1‑136 (MQ=255)
                                                         ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                                                                           >  2:827323/1‑149 (MQ=255)
                                                               aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTc                                                                                                                                   <  1:243689/103‑1 (MQ=255)
                                                                tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTc                                                                                    <  2:1400363/149‑1 (MQ=255)
                                                                 cACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCa                                                                                   <  1:231013/149‑1 (MQ=255)
                                                                  aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAAt                                                                                                >  1:1390991/1‑135 (MQ=255)
                                                                  aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCAt                                                                                  <  4:890628/149‑1 (MQ=255)
                                                                      aCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCa                                                                              >  2:972627/1‑149 (MQ=255)
                                                                         gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaaa                                                                           <  3:376560/149‑1 (MQ=255)
                                                                            ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                                                                        >  1:1247367/1‑149 (MQ=255)
                                                                              aaaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCTAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                                                                      <  2:868569/149‑1 (MQ=255)
                                                                              aaaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                                                                      <  4:645445/149‑1 (MQ=255)
                                                                              aaaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                                                                      >  1:1452353/1‑149 (MQ=255)
                                                                              aaaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                                                                      <  1:761632/149‑1 (MQ=255)
                                                                               aaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                                                                     >  1:102755/1‑149 (MQ=255)
                                                                                 ggCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGCCAGCAAGAAGTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAg                                                                   <  4:912090/149‑1 (MQ=255)
                                                                                   cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGaa                                                                 <  1:1129988/149‑1 (MQ=255)
                                                                                    ccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCAt                                                                                  <  4:187601/131‑1 (MQ=255)
                                                                                     cTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAgg                                                               >  1:272935/1‑149 (MQ=255)
                                                                                       tCTGCTGCAAACGTAATCGCTCTGTAAGATATGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGtt                                                             <  4:57502/149‑1 (MQ=255)
                                                                                       tCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGtt                                                             <  3:1191384/149‑1 (MQ=255)
                                                                                        ctgctgCAAACGTAATCGCTCTCTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg                                                            <  3:894370/149‑1 (MQ=255)
                                                                                          gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc                                                          <  4:573935/149‑1 (MQ=255)
                                                                                            tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGa                                                        >  4:1336792/1‑149 (MQ=255)
                                                                                              cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaaa                                                      <  2:291517/149‑1 (MQ=255)
                                                                                              cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATCCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaaa                                                      <  4:258371/149‑1 (MQ=255)
                                                                                               aaaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa                                                           >  4:465019/1‑143 (MQ=255)
                                                                                               aaaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa                                                           <  3:465019/143‑1 (MQ=255)
                                                                                               aaaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaaaa                                                     <  3:1427432/149‑1 (MQ=255)
                                                                                               aaaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaaaa                                                     >  3:58330/1‑149 (MQ=255)
                                                                                                aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc                                                    <  2:791581/149‑1 (MQ=255)
                                                                                                 aCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGTCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACg                                                   <  4:1387905/149‑1 (MQ=255)
                                                                                                  cGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGt                                                  <  4:58330/149‑1 (MQ=255)
                                                                                                   gTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                                 <  4:1183254/149‑1 (MQ=255)
                                                                                                   gTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                                 <  2:1452353/149‑1 (MQ=255)
                                                                                                   gTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                                 >  1:1374449/1‑149 (MQ=255)
                                                                                                   gTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                                 >  1:1203948/1‑149 (MQ=255)
                                                                                                     aaTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCa                                               <  1:1106514/149‑1 (MQ=255)
                                                                                                     aaTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCa                                               <  1:861153/149‑1 (MQ=255)
                                                                                                       tCGCTCTGTAAGATACGTCAGCAAGAATTCAATACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGt                                             >  4:1249903/1‑149 (MQ=255)
                                                                                                           tctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGcc                                                                                                     >  1:1459837/1‑89 (MQ=255)
                                                                                                             tGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                                 <  2:1148314/139‑1 (MQ=255)
                                                                                                             tGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGcc                                       >  1:445737/1‑149 (MQ=255)
                                                                                                              gTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg                                                            >  2:560781/1‑127 (MQ=255)
                                                                                                              gTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGccc                                      >  2:1067995/1‑149 (MQ=255)
                                                                                                               tAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCa                                     >  3:599505/1‑149 (MQ=255)
                                                                                                                 aGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                                 >  3:1447885/1‑135 (MQ=255)
                                                                                                                 aGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                                 <  4:1447885/135‑1 (MQ=255)
                                                                                                                 aGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGc                                   >  3:698261/1‑149 (MQ=255)
                                                                                                                    tACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCaa                                >  2:177863/1‑149 (MQ=255)
                                                                                                                           gCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGc                                   <  1:503006/139‑1 (MQ=255)
                                                                                                                           gCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGc                                   >  2:503006/1‑139 (MQ=255)
                                                                                                                             aaGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCct                       >  2:571542/1‑149 (MQ=255)
                                                                                                                                caaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTTTGCGGTTGGGGCCGAAAACAGCCAAAATGCAAGTGCGGTCATGGCAAAGGCCCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCtt                    >  3:248300/6‑149 (MQ=255)
                                                                                                                                  ttCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGt                  <  3:1416451/149‑1 (MQ=255)
                                                                                                                                    cAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATCGCAAAAGACCCCAGAAGGTTGACGTAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTaa                <  4:841913/149‑1 (MQ=255)
                                                                                                                                     aaaaCCCGCTGATTAAGCGGGTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaaaa                                                                           <  3:1282424/88‑1 (MQ=38)
                                                                                                                                      aaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACa              >  1:868473/1‑149 (MQ=255)
                                                                                                                                       aaCCCGCTTAATCAGCGGGTTTTTATTGGTCTTTAGTGTGCGGTTGAGGCCGAAACCAGCCAGAATGCTAGTGCGGTCAAGGCAATAGATCCCAGAAGGTTGACGAAAACGTGCAGTAAGGCGTAGCCAAAGCGGCCCTCTTGTAACaa             >  4:1223761/1‑149 (MQ=255)
                                                                                                                                        aCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaa            >  1:450042/1‑149 (MQ=255)
                                                                                                                                          ccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaaaa          >  1:1182931/1‑149 (MQ=255)
                                                                                                                                          ccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaaaa          <  1:398731/149‑1 (MQ=255)
                                                                                                                                          ccGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaaaa          >  3:198565/1‑149 (MQ=255)
                                                                                                                                             cTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAacac       <  3:1033938/149‑1 (MQ=255)
                                                                                                                                               tAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAcacca     <  2:934230/149‑1 (MQ=255)
                                                                                                                                                   cAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                                 >  3:945650/1‑101 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTT  >  NC_000913/657388‑657682

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