Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,473,900 C→T G52R (GGA→AGA)  bdcA ← c‑di‑GMP‑binding biofilm dispersal mediator protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,473,9000CT98.0% 159.7 / NA 50G52R (GGA→AGA) bdcAc‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/25);  minor base G (1/0);  total (25/25)
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.87e-01

ACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCT  >  NC_000913/4473656‑4474144
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
aCGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTAAGCACTATCTGTGAATACTGCTGTCGCTCTAGtctc                                                                                                                                                                                                                                                >  1:267472/1‑251 (MQ=255)
   gAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTtg                                                                                                                                                                                                                                             <  1:269525/251‑1 (MQ=255)
           gatTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGc                                                                                                                                                                                                                                     <  2:137412/251‑1 (MQ=255)
                    gcgcCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGCACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCgg                                                                                                                                                                                                                            >  2:88928/1‑251 (MQ=255)
                    gcgcCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCgg                                                                                                                                                                                                                            <  2:368903/251‑1 (MQ=255)
                         cgccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCt                                                                                                                                                                                                                       <  1:83370/251‑1 (MQ=255)
                               cGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGAc                                                                                                                                                                                                                 >  1:459410/1‑251 (MQ=255)
                                 ggCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTTTGAATACTGCTGTCGCTCTAGTCGCTTGTGCCAGGCGTTTAGCGGCATCTTTCGAccc                                                                                                                                                                                                               <  2:465469/251‑1 (MQ=255)
                                  gCATCTGCAGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGAcccc                                                                                                                                                                                                              <  1:264825/251‑1 (MQ=255)
                                        tccGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGAGCAAGATCGTCGGCATTTAAGTCCAGGGCCTCGCCAGAGAAGCCAATACCTGCATTGACCGCCAGGATATCCAAGGCGCCTCTCTTACGAACGACATTAATGACCGCGTCTCGGTCAGCACTAGCGGTGAATACTGCTGTCGCGCTAGGCTGTTGTGCCAGGCGTTTATCGTCATCGTTCGACCCAGCATa                                                                                                                                                                                                         <  2:451378/249‑1 (MQ=255)
                                        gCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATAGTGATTGTGAAAAGGCGATCAATATCGTCGGCAGTTAATTCCAGGGCCTCGCCAAAGACGCCAATATCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCAt                                                                                                                                                                                                          <  1:505923/249‑1 (MQ=255)
                                            ggcggcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTg                                                                                                                                                                                                     >  2:125681/1‑250 (MQ=255)
                                                   tCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGCTATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                                                                                                                                             >  2:438432/1‑251 (MQ=255)
                                                   tCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTa                                                                                                                                                                                              <  2:274723/250‑1 (MQ=255)
                                                              aTGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAt                                                                                                                                                                                  <  1:435263/251‑1 (MQ=255)
                                                                    aGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTg                                                                                                                                                                             >  1:235851/1‑250 (MQ=255)
                                                                            gAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATGTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCAGTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGACCCATCTGGGACGAAacg                                                                                                                                                                     <  2:364006/250‑1 (MQ=255)
                                                                                      tttGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATa                                                                                                                                                           <  1:407280/250‑1 (MQ=255)
                                                                                       ttGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAg                                                                                                                                                          <  1:313790/250‑1 (MQ=255)
                                                                                         gAAAAGGCGCTCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATAACTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAACCGACGTACGATAGCgg                                                                                                                                                       >  2:334796/1‑251 (MQ=255)
                                                                                         gAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCg                                                                                                                                                        <  1:356121/250‑1 (MQ=255)
                                                                                           aaaGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGc                                                                                                                                                      >  2:269927/1‑250 (MQ=255)
                                                                                            aaGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTAGACCCCGCATAGGTGAATCGTACATTGGCCCCCTCGGTGACGAAACGACGTACGATAGCGGccc                                                                                                                                                    >  2:467130/1‑249 (MQ=255)
                                                                                             aGGCGATCAATATCTTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAAGGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCAcc                                                                                                                                                   <  2:135865/251‑1 (MQ=255)
                                                                                               gCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCTGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCg                                                                                                                                                  >  1:363955/1‑250 (MQ=255)
                                                                                                 ggacAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCTAAGACTTTCAGTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCATAGACAGCGTCTCTTTCAGCCCTATCTGTGAATACTGCTGTCGCTCTAGGCTCGTGTGCCAGGCGTTTAGCGGCTTCTTTCGGCCGCGAATAGGTGAATCGGACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGTACCGATa                                                                                                                                               <  2:40484/248‑1 (MQ=255)
                                                                                                    cAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGTCCCCATCGGTGACGAAACGGCGTACGATAGCGGCACCGATAcc                                                                                                                                             <  2:315984/250‑1 (MQ=255)
                                                                                                        atCGTAGGCATTTAATTCCAGGGCCTCGCAAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGCACGACATCAATGAAAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGGGAATCGTACATTGGCCCCATCGGTGCCGTAACGCCGTACGATAGCGGCACCGATACCACGAc                                                                                                                                        >  2:347483/1‑251 (MQ=255)
                                                                                                           gtcgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGc                                                                                                                                     >  1:90521/1‑251 (MQ=255)
                                                                                                            tcgGCATTTAATTACAGAGCCTCGCCGAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATTCGCCGCTCTTACGAACGACATCAATGACCGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGGCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCAATCGGTGACGGAACGACGTACGATAGCGGCACCGAGACCACGACTGcc                                                                                                                                    >  2:1502/1‑251 (MQ=255)
                                                                                                                            ggtccTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACt                                                                                                                    <  1:189959/248‑1 (MQ=255)
                                                                                                                                   gCCAAAGACGCCAATACATGCATTTCCCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAc                                                                                                             >  2:176857/1‑251 (MQ=255)
                                                                                                                                         gaccctaacacttGCATTTACCACCAGGATATCCAAGGCGCACCTCTTACGAACGACATCAGTGACAGCGTCTCTGTCAGCACTATCGGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGTCGTTTAGCGGCATCTTTCGACCCCGCATAGGGGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCTCCGATACCACGTCTGCCACCGAGTATGAGATTTGGTTTACCTGTaa                                                                                                       <  2:337534/239‑1 (MQ=255)
                                                                                                                                            gCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGTCCCCGCATAGGTGAATCGTCCATTGGCCCCATCGGTGACGAAACGACTTACGATAGCGGCACCGATACCACGACGGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAg                                                                                                    <  2:182794/251‑1 (MQ=255)
                                                                                                                                                                 caGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCAGAAGGACTCCTTGAt                                                                               <  1:387102/251‑1 (MQ=255)
                                                                                                                                                                     atatCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATATGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGACTCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGCGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATttat                                                                           >  2:326630/1‑251 (MQ=255)
                                                                                                                                                                          cAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGt                                                                      <  2:381333/251‑1 (MQ=255)
                                                                                                                                                                                 ccgcTCTTACGAACGACATCAATGACAGCGTCTCTGTAAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGc                                                               >  2:83064/1‑251 (MQ=255)
                                                                                                                                                                                        tACGAACGACATCAATTATAGCGTCTCTGTCAGCACTATCTGTGAATACTTCTGTCGCTCTATTCTCTTGTGCCAGGCGTTTAGCGGCATTTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCTAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATGATGTAACATGCATTACa                                                         <  2:461507/250‑1 (MQ=255)
                                                                                                                                                                                        tACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACa                                                         <  2:144305/250‑1 (MQ=255)
                                                                                                                                                                                                       tGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGACTCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCCCCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACttt                                         >  1:54643/1‑251 (MQ=255)
                                                                                                                                                                                                                    gTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGCCGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGTACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGtt                            >  2:235725/1‑251 (MQ=255)
                                                                                                                                                                                                                       agctcTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCAGCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGGAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTg                         <  2:83069/247‑1 (MQ=255)
                                                                                                                                                                                                                           cTATCTGTGAATACTGCTGTCGCTCTAGTCTATTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGACACCGAGGATGAGAACTGTCTTAACTGTAAAAGCGACCATAAGGCCTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAAt                     >  2:333538/1‑251 (MQ=255)
                                                                                                                                                                                                                                tgtgAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTc                 >  1:100450/1‑250 (MQ=255)
                                                                                                                                                                                                                                  tgAATACTGCTGTCGCTCTAGTCTCTTGTGCCCGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACt              >  2:323543/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ctgatgTCGCTCGGGTCTCTTGTGCCAGGCGTTTAGCGACATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTTCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTCTGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAAAAGGGTTTGTAATgttacgtaaaaa         >  2:39212/5‑239 (MQ=255)
                                                                                                                                                                                                                                           ctgTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGAACCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGCACCATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACa      >  1:70507/1‑250 (MQ=255)
                                                                                                                                                                                                                                             gTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATTTTTCGACCCCGAATAGGTGAATCGTACCTTGGCCCCATCGGTGACGAAACGACGTACGTTAGCGGCACCGATACCACGAATGCCACCGAGGATGAGAACTGTCTTACATGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAAAATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAt  >  2:451932/1‑250 (MQ=255)
                                                                                                                                                                                                                                              tCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAtct  >  1:504827/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
ACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCT  >  NC_000913/4473656‑4474144

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