Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,313,027 T→G intergenic (‑36/‑323) alsR ← / → rpiB DNA‑binding transcriptional repressor AlsR/allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,313,0270TG98.3% 228.7 / NA 60intergenic (‑36/‑323)alsR/rpiBDNA‑binding transcriptional repressor AlsR/allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (18/41);  total (18/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAG  >  NC_000913/4312789‑4313261
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
gATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGc                                                                                                                                                                                                                                 <  1:36915/250‑1 (MQ=255)
   aCCTTACCTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCa                                                                                                                                                                                                                                >  1:831289/1‑248 (MQ=255)
    ccTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGaa                                                                                                                                                                                                                            <  2:659966/251‑1 (MQ=255)
    ccTTAACTATCACCGCTTCTGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGaa                                                                                                                                                                                                                            <  2:453258/251‑1 (MQ=255)
      ttAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAAc                                                                                                                                                                                                                           >  2:56480/1‑250 (MQ=255)
       tAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAAc                                                                                                                                                                                                                           <  1:820175/249‑1 (MQ=255)
       tAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACa                                                                                                                                                                                                                          <  2:808345/250‑1 (MQ=255)
        aaCTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAg                                                                                                                                                                                                                         <  1:646378/250‑1 (MQ=255)
            atcatcGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCTCCACAGCAAGAACAGAAAAg                                                                                                                                                                                                                    <  1:179415/251‑1 (MQ=255)
               atcGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATg                                                                                                                                                                                                                 <  2:716584/251‑1 (MQ=255)
                  gCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTg                                                                                                                                                                                                              <  1:655404/251‑1 (MQ=255)
                     tCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTgg                                                                                                                                                                                                             <  1:788759/249‑1 (MQ=255)
                     tCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCCCCACAGCAAGAACAGAAAAGATGCTGGCg                                                                                                                                                                                                           >  2:599291/1‑251 (MQ=255)
                     tCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGc                                                                                                                                                                                                            <  1:530282/250‑1 (MQ=255)
                       aGATACCGCCAGAGCTGCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCAGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGAc                                                                                                                                                                                                         <  2:244652/251‑1 (MQ=255)
                          tACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTc                                                                                                                                                                                                      >  2:359213/1‑251 (MQ=255)
                          tACCGCCAGAGCTGCTGCGACATCTTTAATTGCGGGGGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTACTTCCTGCTTCATTCGCAGGAAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTACTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTc                                                                                                                                                                                                      >  1:282152/1‑251 (MQ=255)
                           aCCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTc                                                                                                                                                                                                      <  1:103317/250‑1 (MQ=255)
                               ccAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACaaa                                                                                                                                                                                                 <  1:192554/251‑1 (MQ=255)
                                  gagCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACaaaaaa                                                                                                                                                                                              <  1:425116/251‑1 (MQ=255)
                                  gagCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACaaaaaa                                                                                                                                                                                              <  2:476036/251‑1 (MQ=255)
                                           cGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCt                                                                                                                                                                                                               <  1:404432/225‑1 (MQ=255)
                                           cGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCt                                                                                                                                                                                                               >  2:404432/1‑225 (MQ=255)
                                                        tGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGtttattt                                                                                                                                                                        >  1:13600/1‑251 (MQ=255)
                                                         gCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGtttatttt                                                                                                                                                                       <  1:807557/251‑1 (MQ=255)
                                                                          gTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGa                                                                                                                                                      <  1:783599/251‑1 (MQ=255)
                                                                             aCCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTg                                                                                                                                                   >  2:561324/1‑251 (MQ=255)
                                                                                     tGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTaaa                                                                                                                                            <  1:124021/250‑1 (MQ=255)
                                                                                        gTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCTTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTTAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTTGCTTAAAAttt                                                                                                                                        <  1:375376/251‑1 (MQ=255)
                                                                                            ccACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTAtt                                                                                                                                    <  1:316392/251‑1 (MQ=255)
                                                                                            ccACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTAtt                                                                                                                                    <  1:735829/251‑1 (MQ=255)
                                                                                             cACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTAttt                                                                                                                                   <  1:331416/251‑1 (MQ=255)
                                                                                             cACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTAttg                                                                                                                                   <  2:61174/251‑2 (MQ=255)
                                                                                                    cGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTaa                                                                                                                            <  1:221884/251‑1 (MQ=255)
                                                                                                     gATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCGAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAAt                                                                                                                           >  2:240284/1‑251 (MQ=255)
                                                                                                         cGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGt                                                                                                                       >  1:430282/1‑251 (MQ=255)
                                                                                                            cTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGACGACGCCACTAGCTCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGAGAGTGGCTTAAAATTTTATTTCATTTATTTCGttt                                                                                                                     >  2:448325/1‑250 (MQ=255)
                                                                                                                   tttCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACat                                                                                                             >  2:486090/1‑251 (MQ=255)
                                                                                                                        gTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACatgat                                                                                                          <  1:56480/249‑1 (MQ=255)
                                                                                                                          cATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGtt                                                                                                       >  1:681380/1‑250 (MQ=255)
                                                                                                                          cATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGCAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAATTTTTATTTCCTTTAATTAGTTTAAAAACATGATGttt                                                                                                      >  1:242385/1‑250 (MQ=255)
                                                                                                                              ccttcctGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGTCTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTAttt                                                                                                   <  1:799347/250‑1 (MQ=255)
                                                                                                                                     gCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGttttt                                                                                             <  1:853746/249‑1 (MQ=255)
                                                                                                                                        tCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTgg                                                                                         <  2:427991/250‑1 (MQ=255)
                                                                                                                                        tCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTAGTTGTTTTTCTg                                                                                          <  2:558573/249‑1 (MQ=255)
                                                                                                                                             cGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGttttt                                                                                    <  1:199723/250‑1 (MQ=255)
                                                                                                                                                 ggTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGtttttt                                                                                   >  2:830494/1‑247 (MQ=255)
                                                                                                                                                            aaCCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGtt                                                                     <  1:44457/250‑1 (MQ=255)
                                                                                                                                                              cccTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCt                                                                  <  1:352513/251‑1 (MQ=255)
                                                                                                                                                                                    gAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAAc                                              <  1:479050/249‑1 (MQ=255)
                                                                                                                                                                                                 aCTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCa                               <  2:362240/251‑1 (MQ=255)
                                                                                                                                                                                                  cTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGCTCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTc                                >  2:273098/1‑249 (MQ=255)
                                                                                                                                                                                                    ggCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGAt                            <  2:308392/251‑1 (MQ=255)
                                                                                                                                                                                                     gCTCATCCCGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTGTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATg                           <  2:392709/251‑1 (MQ=255)
                                                                                                                                                                                                         atcaCGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAAc                       <  2:391847/251‑1 (MQ=255)
                                                                                                                                                                                                            aCGTTCTCTTGTCCTTGACTTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAActct                    >  2:711614/1‑251 (MQ=255)
                                                                                                                                                                                                             cGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTg                   >  1:110110/1‑251 (MQ=255)
                                                                                                                                                                                                                            gaattgaaGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAgg     <  1:187315/250‑1 (MQ=255)
                                                                                                                                                                                                                            gaattgaaGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGa    <  1:599291/251‑1 (MQ=255)
                                                                                                                                                                                                                              attgaaGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAg  <  1:399851/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
GATACCTTAACTATCATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAACTCTGTAAAACAGATCAGGAAG  >  NC_000913/4312789‑4313261

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