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,9000CT97.6% 287.4 / ‑3.5 84G52R (GGA→AGA) bdcAc‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (35/47);  total (36/48)
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.79e-01

TTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCC‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCG  >  NC_000913/4473653‑4474146
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
ttCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGt                                                                                                                                                                                                                                                     >  1:35148/1‑251 (MQ=255)
       aGCCGATGATTACGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCCGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCCCTATCTGTGACAACTGCTGTCGCTCT‑AGTCTCTtgt                                                                                                                                                                                                                                              >  2:464665/1‑251 (MQ=255)
        gCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTtgtg                                                                                                                                                                                                                                             <  1:464665/251‑1 (MQ=255)
          cGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGcc                                                                                                                                                                                                                                           >  2:85015/1‑251 (MQ=255)
           gatgatTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCa                                                                                                                                                                                                                                          <  1:215995/251‑1 (MQ=255)
              gatTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGc                                                                                                                                                                                                                                       <  1:604986/251‑1 (MQ=255)
                  aaGATGCGCCCGCCTTCGGGCAGCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGtt                                                                                                                                                                                                                                    <  1:436482/250‑1 (MQ=255)
                   aGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGttt                                                                                                                                                                                                                                   >  1:706964/1‑250 (MQ=255)
                    gATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAg                                                                                                                                                                                                                                 <  1:148147/251‑1 (MQ=255)
                       gcgcCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCgg                                                                                                                                                                                                                              <  2:724740/251‑1 (MQ=255)
                         gcccgccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCa                                                                                                                                                                                                                            <  1:85015/251‑1 (MQ=255)
                           ccgccTTCGGGCTTCTGCCGGGCGTCGTCAACAGAGGCATTATAAGGAGCATGAATATTGATTTTGAAAATGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCACAGACGCCAATACCTGCATTTACCACCTGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATc                                                                                                                                                                                                                          <  2:787022/251‑1 (MQ=255)
                           ccgccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATc                                                                                                                                                                                                                          <  1:471568/251‑1 (MQ=255)
                                 tCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTc                                                                                                                                                                                                                      <  1:801776/249‑1 (MQ=255)
                                 gcGGGCATCTTCCGGGCGGCTTAAACAGAGGCATGATAAGGAGCATGAATATTGAGTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGGCGCGGTCTGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCAGCTTTc                                                                                                                                                                                                                      <  2:404554/248‑1 (MQ=255)
                                           gCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCAt                                                                                                                                                                                                            <  2:257451/249‑1 (MQ=255)
                                                  ggcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATGCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATc                                                                                                                                                                                                   <  1:459426/251‑1 (MQ=255)
                                                  ggcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATAGCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATc                                                                                                                                                                                                   <  2:4148/251‑1 (MQ=255)
                                                   gcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCg                                                                                                                                                                                                  >  2:617482/1‑251 (MQ=255)
                                                        aaCAGAGGCAGGATAAGGAGCATGAATCTTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGATTCGTACAt                                                                                                                                                                                             <  2:561653/251‑1 (MQ=255)
                                                         aCAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATCGGTGAATCGTACAtt                                                                                                                                                                                            >  2:29894/1‑251 (MQ=255)
                                                               gCATGATAAGGAGCATGAATATTGATTTTGCAACGGCGATCCATCTCGTCGGCATTTAATTCCAGATCCTCGCCAAAAACGACAATACCTGCATTTACCTCCAGACTGTCCAATGCACCGCTATTACGAACGACATCAATGACAGCGTCCCTGTCAGACCTATCTGTGATTCCTGCGGT‑TCTCTTAGTCTCTTGTACCAGGCGTTTAGCGGCATCATTCGACCCCgtctcggacaatctttcattggcccc                                                                                                                                                                                      >  2:186834/1‑226 (MQ=255)
                                                                      aaGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAGCGGTg                                                                                                                                                                               <  2:35148/251‑1 (MQ=255)
                                                                      aaGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCGCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGt                                                                                                                                                                                <  1:749607/250‑1 (MQ=255)
                                                                       aGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTTGCCACATCGGTg                                                                                                                                                                               >  1:510227/1‑250 (MQ=255)
                                                                       aGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTg                                                                                                                                                                               >  1:639748/1‑250 (MQ=255)
                                                                        ggAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGAc                                                                                                                                                                             >  1:285739/1‑251 (MQ=255)
                                                                              tGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacg                                                                                                                                                                       >  1:502325/1‑251 (MQ=255)
                                                                                    ttGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTAc                                                                                                                                                                 >  1:771557/1‑251 (MQ=255)
                                                                                      gATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACg                                                                                                                                                                >  2:512034/1‑250 (MQ=255)
                                                                                        ttttGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGAt                                                                                                                                                              >  1:435602/1‑250 (MQ=255)
                                                                                        ttttGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATAGGTGACGAAACGACGTACGAt                                                                                                                                                              >  1:446358/1‑250 (MQ=255)
                                                                                         tttGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATa                                                                                                                                                             <  2:706964/250‑1 (MQ=255)
                                                                                           ggAAAAAGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACTCCAATACCTGCAGTTACCACCAGGAGATCCAATTCTCCGCTCGTACGAACGACATCAATTACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCTTCTTTCGACCCCGCATAGGTGAATCATACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCg                                                                                                                                                          <  2:251609/250‑1 (MQ=255)
                                                                                             aaaaGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGc                                                                                                                                                        <  1:175119/251‑1 (MQ=255)
                                                                                                aGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCAcc                                                                                                                                                     <  1:208743/251‑1 (MQ=255)
                                                                                                   cgaacAATATCGTCGGCATTTAATTCCAGTGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACTACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACTTACGATAGCGGCACCGAt                                                                                                                                                  <  1:712533/247‑1 (MQ=255)
                                                                                                      tCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAcc                                                                                                                                               >  1:603665/1‑251 (MQ=255)
                                                                                                       cAATATCGTCGGCATTTAATTCCCGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCC‑CGTCTCTTGTGCCAgg                                                                                                                                                                                                                                        >  2:59177/1‑161 (MQ=255)
                                                                                                       cAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCC‑AGTCTCTTGTGCCAgg                                                                                                                                                                                                                                        <  1:59177/161‑1 (MQ=255)
                                                                                                               tcgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCCTTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGAACCGATACCACGCCTGcc                                                                                                                                      >  2:592035/1‑251 (MQ=255)
                                                                                                                cgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGccc                                                                                                                                     >  1:384032/1‑250 (MQ=255)
                                                                                                                    aTTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCgg                                                                                                                                                                                 >  2:393752/1‑203 (MQ=255)
                                                                                                                    aTTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCgg                                                                                                                                                                                 <  1:393752/203‑1 (MQ=255)
                                                                                                                        aaTTCCAGGGCCTCGCCAAAGACGCCAATACCCGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGa                                                                                                                              <  2:771557/250‑1 (MQ=255)
                                                                                                                          ccccAGGGTCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCCCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTTAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGAGAGCGGCACCGATACCACGACTGCCACCGAGGATga                                                                                                                           <  1:737467/249‑1 (MQ=255)
                                                                                                                             cAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATgaga                                                                                                                         >  1:354327/1‑250 (MQ=255)
                                                                                                                             cAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATgaga                                                                                                                         >  1:433972/1‑250 (MQ=255)
                                                                                                                               gggCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACt                                                                                                                      >  1:483351/1‑251 (MQ=255)
                                                                                                                                   cTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTAT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCt                                                                                                                  >  2:154266/1‑251 (MQ=255)
                                                                                                                                        cAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                                                                                             >  1:263269/1‑251 (MQ=255)
                                                                                                                                              cGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGGCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGCCGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTaaaa                                                                                                       <  1:478472/251‑1 (MQ=255)
                                                                                                                                               gCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCtt                                                                                                                                                                                                                                                >  2:64340/1‑113 (MQ=255)
                                                                                                                                               gCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCtt                                                                                                                                                                                                                                                <  1:64340/113‑1 (MQ=255)
                                                                                                                                                           aTTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGCTACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAgg                                                                                           >  2:411275/1‑250 (MQ=255)
                                                                                                                                                              tACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACt                                                                                        <  2:771674/250‑1 (MQ=255)
                                                                                                                                                                 caccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCtt                                                                                    <  1:475333/251‑1 (MQ=255)
                                                                                                                                                                     aGGATATCCAATGCGCCGCGCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCTCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCACATAAGGACTCCTTGAtt                                                                                <  1:592035/251‑1 (MQ=255)
                                                                                                                                                                      ggATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAttt                                                                               >  2:530045/1‑251 (MQ=255)
                                                                                                                                                                        atatCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACa                                                                                                                                                                                              >  1:208580/1‑138 (MQ=255)
                                                                                                                                                                        atatCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACa                                                                                                                                                                                              <  2:208580/138‑1 (MQ=255)
                                                                                                                                                                        atatCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATttat                                                                             >  1:19909/1‑251 (MQ=255)
                                                                                                                                                                          atCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTAT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATttatta                                                                           <  1:154266/251‑1 (MQ=255)
                                                                                                                                                                            ccAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTtattat                                                                          <  1:314772/250‑1 (MQ=255)
                                                                                                                                                                              aaTGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCCCCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGt                                                                        >  2:73632/1‑250 (MQ=255)
                                                                                                                                                                              aaTGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTa                                                                       <  1:544402/251‑1 (MQ=255)
                                                                                                                                                                               aTGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTaa                                                                      <  1:697062/251‑1 (MQ=255)
                                                                                                                                                                                          ttACGAACGACATCAATGACAGCGTCTCTGTCAGCACTAGCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTAc                                                            <  2:263269/250‑1 (MQ=255)
                                                                                                                                                                                                acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTg                                                     <  2:603665/251‑1 (MQ=255)
                                                                                                                                                                                                           gACATCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTc                                          <  2:433972/251‑1 (MQ=255)
                                                                                                                                                                                                           gACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTc                                          <  1:230329/251‑1 (MQ=255)
                                                                                                                                                                                                                    tctGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTTAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAAc                                   >  1:99876/1‑249 (MQ=255)
                                                                                                                                                                                                                    tctGTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGAGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAAc                                   >  2:791773/1‑249 (MQ=255)
                                                                                                                                                                                                                       gTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAgg                                <  2:235232/249‑1 (MQ=255)
                                                                                                                                                                                                                       gTCAGCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAg                                 <  1:29894/248‑1 (MQ=255)
                                                                                                                                                                                                                         cAGCACTATCTTTGAATACTGCTGTCGCTCT‑AGTCGCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCACATCGGTGCCGAAACGACGTACGATAGCGGCACCGCTACCCAGACGGCCACCGGGGATGAGAACGGTCTTACCTGTAAACGCGCCCAACAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTAAACAGGtttt                            >  2:119393/1‑251 (MQ=255)
                                                                                                                                                                                                                         cAGCACGATCTGTGAATACTGCTGGCGCTCT‑AGTGTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCTCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGtttt                            <  2:167460/251‑1 (MQ=255)
                                                                                                                                                                                                                           gCACTATCTGTGAATACTGCTGTCGCTCT‑AGTCGCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGt                          <  1:410683/251‑1 (MQ=255)
                                                                                                                                                                                                                            cACTATCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGt                          >  2:333247/1‑250 (MQ=255)
                                                                                                                                                                                                                                aTCTGTGAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAAt                       <  1:617482/249‑1 (MQ=255)
                                                                                                                                                                                                                                    gtgAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGt                    <  1:160782/248‑1 (MQ=255)
                                                                                                                                                                                                                                     tgAATACTGCTGTCGCTCT‑AGTCTCTTGTGCCAGGCGTTTAGAGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACt                >  1:339644/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aatcctgctgTCGCTCT‑AGTCTCTTGTGCCAGTCGTTTAGCGGCATCTTTCGACCCCGCATAAGTGAATCGTACATTGGCCCCATCGGTGTCGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGTCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTaa              <  2:57067/247‑1 (MQ=255)
                                                                                                                                                                                                                                                   gCTCT‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTGCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTcg  <  1:523368/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
TTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCC‑AGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCG  >  NC_000913/4473653‑4474146

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