Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 657,533 A→T 6.8% 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.8% 298.1 / 6.9 108intergenic (+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 (47/52);  new base T (3/5);  total (50/58)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

CTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACT  >  NC_000913/657390‑657681
                                                                                                                                               |                                                                                                                                                     
ctctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCg                                                                                                                                                 <  1:314318/150‑1 (MQ=255)
   tGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGt                                                                                                                                             >  1:1077767/1‑150 (MQ=255)
    gCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                                                                            >  2:642354/1‑150 (MQ=255)
    gCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtt                                                                                                                                            >  2:1778654/1‑150 (MQ=255)
      aaTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTTCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtttt                                                                                                                                          <  2:592309/150‑1 (MQ=255)
      aaTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtttt                                                                                                                                          <  1:1159790/150‑1 (MQ=255)
      aaTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGATTAAGCGGGtttt                                                                                                                                          <  2:220475/150‑1 (MQ=255)
       aTCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGttttt                                                                                                                                         >  1:1221463/1‑150 (MQ=255)
         ccAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGAAAGATACGTCAGCAAGAATTCAAAACCCCCTTAATCAGCGGGttttttt                                                                                                                                       >  1:1186645/1‑150 (MQ=255)
          cAGACTAATGGTTTTAAAACTCTTGCTGAAGGGCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                                                                      <  2:1401429/150‑1 (MQ=255)
          cAGACTAATGGTCTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGAAAGAATTCAAAACCCGCTTAATCAGCGGGtttttttt                                                                                                                                      >  2:1518148/1‑150 (MQ=255)
           ataCTAATTGTTTTATAACTCTTGCTGAATGTCAGCGCGTAGAGTTCGTAATCACTAACGGTGCCAAAGGGCCTTATGCTGCAAACGTAATCGCTCTTTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                                                                     <  1:675072/148‑1 (MQ=255)
           aGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                                                                     <  1:839823/150‑1 (MQ=255)
           aGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTg                                                                                                                                     <  2:189233/150‑1 (MQ=255)
            gACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTgg                                                                                                                                    <  1:483309/150‑1 (MQ=255)
            gACTAATGGTTTTAAAACTCTTGCTGAAGGGCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCTAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTgg                                                                                                                                    <  1:1312778/150‑1 (MQ=255)
                aaTGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAAATCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCtt                                                                                                                                <  1:131893/150‑1 (MQ=255)
                 aTGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCttt                                                                                                                               <  1:189799/150‑1 (MQ=255)
                  tGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTa                                                                                                                              >  1:1565408/1‑150 (MQ=255)
                    gTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAgt                                                                                                                            >  2:451577/1‑150 (MQ=255)
                      tttAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCGTAATCAGCGGGTTTTTTTTGGTCTTTAgtgt                                                                                                                          >  2:1619914/1‑150 (MQ=255)
                       ttAAAACTCTTGCTGAAGGTCAGCGCGTAGAGGTCGTAATCACAAACGGTGCCTAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTAAATCAGCGGGTTTTTTTTGGTCTTTAGtgtg                                                                                                                         <  1:1757450/150‑1 (MQ=255)
                        tAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTATtgtgg                                                                                                                        >  2:356790/1‑149 (MQ=255)
                        tAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGc                                                                                                                        >  1:465803/1‑150 (MQ=255)
                          aaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCgg                                                                                                                      <  1:811024/150‑1 (MQ=255)
                          aaaCTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCgg                                                                                                                      <  1:760602/150‑1 (MQ=255)
                              tctTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGa                                                                                                                  >  2:1363141/1‑150 (MQ=255)
                                  gCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCg                                                                                                             <  2:188532/150‑1 (MQ=255)
                                        ggTCAGCGCGTAGAGTTCGAAATCACTAACTGTGCCAAAGGCCCTTATTCTGCAAACGTAATCGCTCTGTAATATACGTCAGCAAGAATTCAAAACCCGCTTAATAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAc                                                                                                        <  2:263430/150‑1 (MQ=255)
                                         gTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTGAATAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACa                                                                                                       >  2:270075/1‑150 (MQ=255)
                                          tCAGCGCGTAGAGTTCGAAATCACTAACTGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTAAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAg                                                                                                      <  2:909243/150‑1 (MQ=255)
                                          tCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAg                                                                                                      >  1:475158/1‑150 (MQ=255)
                                          tCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAg                                                                                                      <  2:1369155/150‑1 (MQ=255)
                                            aGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGCGGCCGAAAACAGcc                                                                                                    >  2:1265752/1‑150 (MQ=255)
                                                  tAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATg                                                                                              <  2:959025/150‑1 (MQ=255)
                                                     agTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCa                                                                                           >  2:1107136/1‑150 (MQ=255)
                                                      gTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                                                                          <  2:1284859/150‑1 (MQ=255)
                                                      gTTCGAAAGCACTAACGGTGCCACAGGCCCTTCTGCATCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTATTCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAg                                                                                          <  2:82686/150‑1 (MQ=255)
                                                         cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTATTCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGc                                                                                       >  2:808607/1‑150 (MQ=255)
                                                         cGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGc                                                                                       <  2:1767105/150‑1 (MQ=255)
                                                          gAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCg                                                                                      >  2:1632426/1‑150 (MQ=255)
                                                            aaTCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGt                                                                                    <  2:172126/150‑1 (MQ=255)
                                                              tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCa                                                                                  >  1:994192/1‑150 (MQ=255)
                                                              tCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCa                                                                                  >  2:931712/1‑150 (MQ=255)
                                                               cACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGATCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGTGTTTTTTTGGTCTTTAGTGTGCGGTTGAGTCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                                                                                >  2:1329283/1‑150 (MQ=255)
                                                                aCTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATg                                                                                >  2:900530/1‑150 (MQ=255)
                                                                 cTAACGGTGCTAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATgg                                                                               >  1:939715/1‑150 (MQ=255)
                                                                  tAACTGTGCCAAAGGCCCTTATGCTGCAATCGTAATCTCTCTGTAATATACGTCAGCAAGTATTCAAAACCCGCTTAATAAGCGGGTTTTTTTTTGTCTTTAGTGTGCGGTTGAGGCCGAATACAGCCAGAATGCCAGTGCGGTCATGGc                                                                              <  2:672420/150‑1 (MQ=255)
                                                                    aCGGTGCCAAAGTCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGTTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCaa                                                                            <  1:987379/150‑1 (MQ=255)
                                                                       gTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAg                                                                         <  2:1514875/150‑1 (MQ=255)
                                                                          ccAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAc                                                                       <  2:1723311/150‑1 (MQ=255)
                                                                            aaaGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGAcccc                                                                    <  2:566461/150‑1 (MQ=255)
                                                                                 cccTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCACCGGGTTTTTTTTGCTCTTTTGTGTGCGGTTGAGGACGAAAACAGCCAGAATGCCTGTGCGGTCAAGGCAAAAGACCAcagcag                                                               >  2:1233498/1‑150 (MQ=255)
                                                                                     tCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTg                                                           >  2:97676/1‑150 (MQ=255)
                                                                                      ctgctgCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGa                                                          <  2:576386/150‑1 (MQ=255)
                                                                                       tgctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAc                                                         <  2:638175/150‑1 (MQ=255)
                                                                                       tgcgccAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg                                                        <  2:1186645/145‑1 (MQ=255)
                                                                                        gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg                                                        <  2:465803/150‑1 (MQ=255)
                                                                                        gctgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg                                                        <  2:804297/150‑1 (MQ=255)
                                                                                        gcggcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACg                                                        <  2:1226504/147‑1 (MQ=255)
                                                                                         ctgcAAACGTAATCGCTCTGTAAGATACGGCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGa                                                       <  1:1243516/150‑1 (MQ=255)
                                                                                          tgcAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGaa                                                      >  1:45834/1‑150 (MQ=255)
                                                                                             aaaCTTAATGGCTCTGTAAGATACCTCAGCAAGAATTCAACACCCGCTAAATCATCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAAc                                                   <  2:330433/150‑1 (MQ=255)
                                                                                             aaaCGTAATCGCTCTGTAAGATACGGCAGCAAGAATTCAAAACCCGCTTACTCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGGTGACGAAAAc                                                   <  2:916359/150‑1 (MQ=255)
                                                                                              aaCGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGAGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACg                                                  >  2:413724/1‑150 (MQ=255)
                                                                                                cGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGtt                                                >  1:1226536/1‑150 (MQ=255)
                                                                                                 gTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTc                                               <  1:1701615/150‑1 (MQ=255)
                                                                                                  tactCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCa                                              <  1:1499951/147‑1 (MQ=255)
                                                                                                  tAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCa                                              >  2:1540628/1‑150 (MQ=255)
                                                                                                   aaTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCTGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCACAGAAGGTTTACGAAAACGTTCAg                                             >  1:361149/1‑150 (MQ=255)
                                                                                                    aTCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGt                                            >  1:1692334/1‑150 (MQ=255)
                                                                                                     tCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGCCCCCAGAAGGTTGACGAAAACGTTCAGTa                                           >  2:1663522/1‑150 (MQ=255)
                                                                                                        ctctGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATg                                        <  1:1047946/150‑1 (MQ=255)
                                                                                                            gTAATATACGTCATCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCa                                    <  2:1622193/150‑1 (MQ=255)
                                                                                                            gTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGATAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGGTGACGAAAACGTTAAGTAATGgccc                                    >  1:352851/1‑149 (MQ=255)
                                                                                                            gTAAGATACGTCAGCAAGAATCCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCa                                    <  2:1100315/150‑1 (MQ=255)
                                                                                                             tAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAg                                   >  1:1823957/1‑150 (MQ=255)
                                                                                                             tAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAg                                   <  2:1701911/150‑1 (MQ=255)
                                                                                                                gATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCa                                >  2:454074/1‑150 (MQ=255)
                                                                                                                    cGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGc                            >  1:1444229/1‑150 (MQ=255)
                                                                                                                       cAGCAAGAATTCAAAAACCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAATGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGc                         <  1:1644008/150‑1 (MQ=255)
                                                                                                                         gCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACCGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGccc                       >  1:1035192/1‑150 (MQ=255)
                                                                                                                          cAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTAGAGGCCGAAAACAGCCAGAATGCCAGTCCGGACATGGAAAAAGACCCCAGAAGGTTGCCGAAAACGTTCAGTAATGCCCAGCCAAAGAGGcccg                      >  2:132957/1‑149 (MQ=255)
                                                                                                                             gAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCATAGCGGCCCTCtt                   >  1:1357113/1‑150 (MQ=255)
                                                                                                                               aaaaaaaaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGt                 >  1:778112/5‑150 (MQ=255)
                                                                                                                               aTTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGCGCGTTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCATCTTGt                 >  2:419747/1‑150 (MQ=255)
                                                                                                                                 tCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCAAAAGCGGCCCTCTTGTaa               >  2:382742/1‑150 (MQ=255)
                                                                                                                                 tCAAAACCAGCTTCATCAGCGGGTTTTTTTTTGTCTTTAGTGTGCGGTTGAGGCCTAAAACAGACAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTaa               <  2:1257609/150‑1 (MQ=255)
                                                                                                                                  cAAAAGCCGCTTAATCAGCGGGTTTTTTTTTGTCTTTAGTGGGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCATAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAAc              <  1:1167387/150‑1 (MQ=255)
                                                                                                                                  cAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCGTCTTGTatc              >  2:894748/1‑148 (MQ=255)
                                                                                                                                   aacaCCCGCTTAATCAGCGGGTTTTTTTTGGGCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAACGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACa             <  1:920961/147‑1 (MQ=255)
                                                                                                                                   aaaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACa             <  1:1518148/150‑1 (MQ=255)
                                                                                                                                   aaaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGATGAGGCCGAAAACAGACAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACa             <  1:270075/150‑1 (MQ=255)
                                                                                                                                    aaaCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaa            >  1:438395/1‑150 (MQ=255)
                                                                                                                                     aaCCCGCTGATTAAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAAAAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGAAGAAAACGTTCAGTAATTCCCAGCCAAAGCGGCCCTCTTGTTACaaa           >  1:1577212/1‑150 (MQ=255)
                                                                                                                                      aCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaaa          <  2:1796407/150‑1 (MQ=255)
                                                                                                                                      aCCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACaaaa          <  2:1221463/150‑1 (MQ=255)
                                                                                                                                            ttAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAcacca    <  1:1476707/150‑1 (MQ=255)
                                                                                                                                             tAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAAcaccac   <  1:1491843/150‑1 (MQ=255)
                                                                                                                                              aaTCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTGCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACt  >  2:257673/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAAACCCGCTTAATCAGCGGGTTTTTTTTGGTCTTTAGTGTGCGGTTGAGGCCGAAAACAGCCAGAATGCCAGTGCGGTCATGGCAAAAGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACT  >  NC_000913/657390‑657681

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