Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,535 C→A 6.5% 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 6.5% 335.5 / 8.4 101intergenic (+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 (43/51);  new base A (4/3);  total (47/54)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 101 total aligned reads displayed.

CTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAC  >  NC_000913/657390‑657675
                                                                                                                                                 |                                                                                                                                             
ctctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGc                                                                                                                                            >  6:414964/1‑149 (MQ=255)
 tctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                                                           >  8:2448/1‑149 (MQ=255)
   tGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGACCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCggg                                                                                                                                         >  5:287619/1‑149 (MQ=255)
      aaTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttt                                                                                                                                     <  3:34866/149‑1 (MQ=255)
        tCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttt                                                                                                                                   >  1:428913/1‑149 (MQ=255)
          cAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttttt                                                                                                                                 >  1:113366/1‑149 (MQ=255)
          cAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttttt                                                                                                                                 <  6:48057/149‑1 (MQ=255)
           aGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                                                                >  8:248356/1‑149 (MQ=255)
           aGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                                                                <  5:26073/149‑1 (MQ=255)
           aGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                                                                <  2:225975/149‑1 (MQ=255)
               tAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTc                                                                                                                            >  6:370764/1‑149 (MQ=255)
                 aTGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTc                                                                                                                            <  4:423064/147‑1 (MQ=255)
                 aTGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCtt                                                                                                                          <  5:414964/149‑1 (MQ=255)
                    gTTTTAAAACTCTTGCTGAATGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTa                                                                                                                        <  8:397227/148‑1 (MQ=255)
                    gTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTa                                                                                                                        <  2:428913/148‑1 (MQ=255)
                     ttttAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAgt                                                                                                                      <  5:265261/149‑1 (MQ=255)
                     ttttAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAgt                                                                                                                      <  7:2448/149‑1 (MQ=255)
                       ttAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAgt                                                                                                                      >  3:179426/1‑148 (MQ=255)
                       ttAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAgt                                                                                                                      <  4:179426/148‑1 (MQ=255)
                       ttAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAgtgt                                                                                                                    >  2:259681/1‑149 (MQ=255)
                        tAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                                                                   <  4:238504/133‑2 (MQ=255)
                        tAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGttttg                                                                                                                                   >  3:238504/1‑132 (MQ=255)
                         aaaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                                                               >  4:51141/1‑136 (MQ=255)
                         aaaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                                                               <  3:51141/136‑1 (MQ=255)
                         aaaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGc                                                                                                                  >  4:291036/1‑149 (MQ=255)
                         aaaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGATCTGTAAGATACTTCAGCAAGAATTCAAAACCCGCTTAATAAGCGGGTTATTTTTGGGCTTTAGTGTGc                                                                                                                  >  6:200686/1‑149 (MQ=255)
                          aaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCg                                                                                                                 >  8:160840/1‑149 (MQ=255)
                           aaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCgg                                                                                                                <  5:447707/149‑1 (MQ=255)
                            aCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAACTCTCTATCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGTGTTTTTTTGGTCTTTAGTGTGCGGtt                                                                                                              <  8:178460/149‑1 (MQ=255)
                                  gCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGa                                                                                                            <  1:457215/146‑1 (MQ=255)
                                  gCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAg                                                                                                           <  7:261767/147‑1 (MQ=255)
                                   cTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGc                                                                                                         <  7:160840/148‑1 (MQ=255)
                                      aaGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAATCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGTGGcc                                                                                                        <  5:200686/146‑1 (MQ=255)
                                       aGGTCAGCGAGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGaaa                                                                                                    <  8:165204/149‑1 (MQ=255)
                                            aGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGc                                                                                               <  1:239883/149‑1 (MQ=255)
                                             gcgcGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGc                                                                                               >  1:221419/1‑148 (MQ=255)
                                              cgcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGc                                                                                               >  5:196314/1‑147 (MQ=255)
                                               gcgTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAg                                                                                            <  5:334702/149‑1 (MQ=255)
                                                    gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGc                                                                                       <  7:248356/149‑1 (MQ=255)
                                                    gagTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGc                                                                                       <  3:291036/149‑1 (MQ=255)
                                                      gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                                                               <  2:383139/107‑1 (MQ=255)
                                                      gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                                                               >  1:383139/1‑107 (MQ=255)
                                                      gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCa                                                                                     <  8:132321/149‑1 (MQ=255)
                                                       ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACCGCCAGAATGCCAg                                                                                    >  8:254801/1‑149 (MQ=255)
                                                       ttCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                                                                    >  2:418683/1‑149 (MQ=255)
                                                         cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTGTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGt                                                                                   <  1:344902/148‑1 (MQ=255)
                                                         cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTg                                                                                  >  1:362643/1‑149 (MQ=255)
                                                          gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGc                                                                                 >  3:44832/1‑149 (MQ=255)
                                                          gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGc                                                                                 >  3:396180/1‑149 (MQ=255)
                                                           aaaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCg                                                                                <  4:312833/149‑1 (MQ=255)
                                                            aaTGACTAACGGTGCCAAAGGCGCTTCTGCTGCAAACGTAATCGGTCTGTAAGATATGTCAGCAATAATTCATTACCCTCTTATTCAGGGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCg                                                                                <  6:287619/148‑1 (MQ=255)
                                                             aTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                                                              >  1:273735/1‑149 (MQ=255)
                                                              tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTc                                                                             <  6:196314/149‑1 (MQ=255)
                                                                  tAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCAAAAACAGCCAGAATGCCAGTGCGGTCATgg                                                                         >  2:467356/1‑149 (MQ=255)
                                                                  tAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTATTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                                                                         >  6:7648/1‑149 (MQ=255)
                                                                    aCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                                                                         <  5:307720/147‑1 (MQ=255)
                                                                    aCGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGGCGTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGc                                                                        >  5:156037/1‑148 (MQ=255)
                                                                      ggTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaa                                                                      >  8:4031/1‑148 (MQ=255)
                                                                          ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGa                                                                  <  2:378721/148‑1 (MQ=255)
                                                                          ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGa                                                                  >  1:378721/1‑148 (MQ=255)
                                                                          ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                                                                 >  1:327002/1‑149 (MQ=255)
                                                                              aGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCa                                                             <  1:259681/149‑1 (MQ=255)
                                                                                   cTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAgg                                                        <  2:362643/149‑1 (MQ=255)
                                                                                        gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc                                                   >  2:443047/1‑149 (MQ=255)
                                                                                         ctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc                                                   >  5:81854/1‑148 (MQ=255)
                                                                                          tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg                                                  >  1:312741/1‑148 (MQ=255)
                                                                                            cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaaa                                               <  2:221419/149‑1 (MQ=255)
                                                                                            cAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaa                                                >  8:138989/1‑148 (MQ=255)
                                                                                               aCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACg                                            <  5:7648/149‑1 (MQ=255)
                                                                                                 gTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                          <  8:392467/149‑1 (MQ=255)
                                                                                                  tAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTc                                         <  4:396180/149‑1 (MQ=255)
                                                                                                        ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAAt                                   <  4:44832/149‑1 (MQ=255)
                                                                                                        ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAAt                                   >  3:397308/1‑149 (MQ=255)
                                                                                                           tGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGcc                                <  6:156037/149‑1 (MQ=255)
                                                                                                            gTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGccc                               >  8:265977/1‑149 (MQ=255)
                                                                                                                gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc                                                   >  3:384046/1‑125 (MQ=255)
                                                                                                                gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc                                                   <  4:384046/125‑1 (MQ=255)
                                                                                                                gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGGTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGc                            >  5:456673/1‑148 (MQ=255)
                                                                                                                           aaGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                                                               <  3:16043/102‑1 (MQ=255)
                                                                                                                           aaGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAccc                                                               >  4:16043/1‑102 (MQ=255)
                                                                                                                           aaGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAACCGTTCAGTAATGCCCAGCCAAAGCGGCCct                >  3:70344/1‑149 (MQ=255)
                                                                                                                           aaGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCct                >  7:363702/1‑149 (MQ=255)
                                                                                                                           aaGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCct                <  7:265977/149‑1 (MQ=255)
                                                                                                                            aGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCctc               <  2:327002/149‑1 (MQ=255)
                                                                                                                             gAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCctct              >  8:376258/1‑149 (MQ=255)
                                                                                                                              caaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAg                       <  7:81848/134‑1 (MQ=255)
                                                                                                                              aaTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCtt             <  7:138989/149‑1 (MQ=255)
                                                                                                                              aaTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCAAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCtt             <  1:467356/149‑1 (MQ=255)
                                                                                                                                ttCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGt           <  4:397308/149‑1 (MQ=255)
                                                                                                                                aaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                          >  2:436525/4‑118 (MQ=255)
                                                                                                                                aaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                          <  1:436525/115‑1 (MQ=255)
                                                                                                                                  cAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCTGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTa          <  5:102743/148‑1 (MQ=255)
                                                                                                                                   aaaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGc                      <  8:64000/135‑1 (MQ=255)
                                                                                                                                   aaaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGc                      >  7:64000/1‑135 (MQ=255)
                                                                                                                                   aaaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGc                      <  2:203452/135‑1 (MQ=255)
                                                                                                                                   aaaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGc                      >  1:203452/1‑135 (MQ=255)
                                                                                                                                   aaaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTaa         <  2:113366/148‑1 (MQ=255)
                                                                                                                                   aaaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAAc        <  1:418683/148‑1 (MQ=255)
                                                                                                                                    aaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAAc        >  7:462629/1‑148 (MQ=255)
                                                                                                                                          tcTTAATCATCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACATCCAGAATGCCCGTGCGGGCAATGCAAACGACCACCGAATGTTTAAGATAACGTTCAGTTATGCCCAGGCAAAGCGGACCTCTTTTAAGAAacac  >  6:355581/2‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAC  >  NC_000913/657390‑657675

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