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,9000CT95.0% 235.4 / 3.1 80G52R (GGA→AGA) bdcAc‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (39/37);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

GATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCG  >  NC_000913/4473667‑4474146
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
gatTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGGGCCAGGc                                                                                                                                                                                                                                       >  1:222746/1‑251 (MQ=255)
gactcaGATGCACCCGCCTTCGGGCATCTGCCGGGCGGCTTCTACAGAGGCATGCTAAGGCGCATGAATATTGATTTTGTAAAGGCGATCAATATCGTCGGCCTTTAATTACAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACAAGGATAT‑CCAATGCTCCGCCCTTACTACAGACATCAATGACCGCGTCCCAGTCAGCCCTATCTGTGAATACTGCTGTCTCACTAGTCTCATGTGCCAGGc                                                                                                                                                                                                                                       >  2:434214/6‑251 (MQ=255)
         gcgcCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCg                                                                                                                                                                                                                               <  1:502558/250‑1 (MQ=255)
         gcgcCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGGCGCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCCCTTGTGCCAGGCGTTTAGCg                                                                                                                                                                                                                               <  2:632205/250‑1 (MQ=255)
               gccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGACCGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATCCTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCtt                                                                                                                                                                                                                        >  2:691981/1‑251 (MQ=255)
                  ttCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCg                                                                                                                                                                                                                     >  1:235262/1‑251 (MQ=255)
                      ggCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGCGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGCCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCACTTGAGCCAGGCGGTTAGTGGCATCTTTCGAccc                                                                                                                                                                                                                 >  2:284310/1‑251 (MQ=255)
                           cTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCAt                                                                                                                                                                                                            >  2:302042/1‑251 (MQ=255)
                                 ggcggcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAAGCGATCAATATAGTCGGCATTTCATTCACGGGCCTCGCCACAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGACGCTATTAAGAACCACATCAATGACAGCGTCTCTGTCAGCACTCTCTGTGAATACTGCCGTCGCTCTCGTCTCTTTTGCCCGGCGTTTAGAGGCATCTTTCGAACCCGCATAGGTgg                                                                                                                                                                                                      >  2:410435/1‑250 (MQ=255)
                                     gcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCg                                                                                                                                                                                                  >  1:547210/1‑251 (MQ=255)
                                         cAACAGAGGCATGATAAGGAGAATTAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                                                                                                                                               <  2:193948/250‑1 (MQ=255)
                                                   aTGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAt                                                                                                                                                                                    >  1:598810/1‑251 (MQ=255)
                                                         aGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGa                                                                                                                                                                              <  1:519847/251‑1 (MQ=255)
                                                          gtaGCATGAATATTGATTATGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGAc                                                                                                                                                                             <  1:410435/249‑1 (MQ=255)
                                                             gCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCCCTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGACAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCgggcacca                                                                                                                                                                           >  2:460876/1‑244 (MQ=255)
                                                             gCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCa                                                                                                                                                                                                                                          <  2:170366/187‑1 (MQ=255)
                                                             gCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCa                                                                                                                                                                                                                                          >  1:170366/1‑187 (MQ=255)
                                                             gCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGa                                                                                                                                                                           >  1:541409/1‑250 (MQ=255)
                                                                tGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCACTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGCAacg                                                                                                                                                                       >  1:481670/1‑251 (MQ=255)
                                                                                aaaGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATCGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGCCGTACGATAGCGGc                                                                                                                                                        >  2:546763/1‑250 (MQ=255)
                                                                                aaaGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCa                                                                                                                                                       <  1:558809/251‑1 (MQ=255)
                                                                                  aGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAAAGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGAACATTGGCCCCATCGGTGACGAAACGACGTCCGATAGCGGCAcc                                                                                                                                                     >  2:366518/1‑251 (MQ=255)
                                                                                         aaaTATCGTAGGCATTTAAATCCAGGGCCTCGCCAAAGACGCCAATACCCGCATTTACCTCCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCa                                                                                                                                              <  2:598810/250‑1 (MQ=255)
                                                                                            tatCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGACTACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCCTCGGTGACGAAACGGCGGACGATAGCGGccacgctaccacc                                                                                                                                            >  2:129344/1‑238 (MQ=255)
                                                                                                gtcgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGCCTGc                                                                                                                                       >  2:132689/1‑251 (MQ=255)
                                                                                                  cgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGcc                                                                                                                                      <  1:698206/250‑1 (MQ=255)
                                                                                                  cgGCATTTAATTCCAGGCCCTCGCAAAAGACGTCAATACCTGCATTTACCACTCGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGCCAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGAATGcc                                                                                                                                      >  2:333393/1‑250 (MQ=255)
                                                                                                         tAATTCCAGGGCCCCGACAATGACGCCAATACCTGCATTTACCACCATG‑TGTACCTATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATATCGGCACCGATACCACGACTGCCACCGAgg                                                                                                                               <  2:524588/250‑1 (MQ=255)
                                                                                                           aTTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                                                                                            <  1:101866/251‑1 (MQ=255)
                                                                                                              ccAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                                                                                            >  1:152238/1‑248 (MQ=255)
                                                                                                               cAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATgaga                                                                                                                         >  1:481792/1‑250 (MQ=255)
                                                                                                                     cTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCCTAGGTGAATCGTACATTGGCCCCATCGGGGCCGAAACGACGTACGCTAGCGGCCCCGATACCACGACTGCCACCGAGGATGCGAACTGTCt                                                                                                                  >  2:200977/1‑251 (MQ=255)
                                                                                                                         ccAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATTCGCCGCTCTTACGAACCACATCAAGGACAGCGTCTCTGTCCGCACTAGCTGTGAATACTGCTGTCGCTCTAGTCCCTTGTGCCAGGCGTTTAGCTGCATCTTTCGACCCCGCATAGGTTAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAACACGACTGCCACCGAGGATGAGAACTGTCttgcc                                                                                                              <  2:237924/251‑4 (MQ=255)
                                                                                                                          cAAAGACGCCAATACCTGCATTTACCTCCAGGATAT‑CCAATTCGCCGCTATTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGGTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGTGGACGAAACGACGTACGATAGCTGCACCGATACCTCGACTGCCACCGAGGATGAGAACTGTCTTAcgt                                                                                                             <  2:622892/251‑3 (MQ=255)
                                                                                                                          cAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAcc                                                                                                              >  1:535129/1‑250 (MQ=255)
                                                                                                                          cAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAcc                                                                                                              <  2:535129/250‑1 (MQ=255)
                                                                                                                          cAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                                                                                             >  1:121825/1‑251 (MQ=255)
                                                                                                                             aGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACCGTCTTACCTGt                                                                                                           >  2:360703/1‑250 (MQ=255)
                                                                                                                                cGCCAATACCTGCATTTCCCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGGCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTaaaa                                                                                                       <  2:541409/251‑1 (MQ=255)
                                                                                                                                           gCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATaa                                                                                             <  1:132689/250‑1 (MQ=255)
                                                                                                                                           gCATTTACCACCAGGATAT‑ACAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATa                                                                                              <  1:502765/249‑1 (MQ=255)
                                                                                                                                            cATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGCACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTCAAAGCGCCCATAAgg                                                                                           >  2:708649/1‑251 (MQ=255)
                                                                                                                                                tACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAc                                                                                                                                                <  1:331097/194‑1 (MQ=255)
                                                                                                                                                tACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAc                                                                                                                                                >  2:331097/1‑194 (MQ=255)
                                                                                                                                                tACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTc                                                                                       >  1:677881/1‑251 (MQ=255)
                                                                                                                                                  ccaccaGGATTT‑CCAATGCGCCGCTCTAACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCt                                                                                     <  2:481792/251‑1 (MQ=255)
                                                                                                                                                   ccccaGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACCGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTg                                                                                                            <  2:377044/225‑1 (MQ=255)
                                                                                                                                                   caccaGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTg                                                                                                            >  1:377044/1‑227 (MQ=255)
                                                                                                                                                   caccaGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCtt                                                                                    <  1:298989/251‑1 (MQ=255)
                                                                                                                                                      caGGATAG‑CCATTGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCTGCACGTTCTGTGAATACTGCTGTCGCTATAGTCGCTTGTGACAGGCGTTTAGCGGCATCTTCAGACCCCGCATAGGTGAATCGTAAATTGGCCCCATCGGTGACGAAACGAAGTACGATAGCGGCACCGATACCACGACTGCCTCCTAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACACCTTGAt                                                                                 <  2:17544/251‑1 (MQ=255)
                                                                                                                                                        ggATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAttt                                                                               <  1:708649/251‑1 (MQ=255)
                                                                                                                                                             acccAATGCGCCGCCCTGACGAGAGACATCAATGTCAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCATGCGTTTATCGGCATCTTTCGACCCCTCATAGGGGAATTGTACATTGGCCCCAGCGCTGACGAAACGACGTACGATAGCGGGCCCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATttatg                                                                            <  2:558277/248‑2 (MQ=255)
                                                                                                                                                                       ccgcTCTTACGAACGAGATCAATGACATAGTCGCTGTCAGCACTAGCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCTGGCTTTTAGCTGCAGCTGGCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGTTGACGAAACGACGTACGATAGCGGCACCGATACCTCGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGTCTCCTTGATTTATTATGGAACATGc                                                                 <  2:97594/251‑1 (MQ=255)
                                                                                                                                                                       ccgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGc                                                                 >  1:115183/1‑251 (MQ=255)
                                                                                                                                                                        cgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCa                                                                <  1:208969/251‑1 (MQ=255)
                                                                                                                                                                         gcTCTTACGACCGACATCAATGCCAGCGTCGCTGTCAGCACTATCTGTGAATACTGCTGCCGCTCTAGTCTCTTGTGCCAGGGGTTTAGCGGCATATTTCGCCCCCCCATAGGTGAATCGTACATTGGCCCCATCGGTGACGACACGACGTACGCAAGCGGCACAGATACCCCGACTGCCCCCGAGGATGAGAACAGTATTACCTGTTAACGCGCCCATAAGGACTCCTTGGTTTATTATGTAacatacat                                                               >  2:97048/1‑251 (MQ=255)
                                                                                                                                                                         gcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTGCGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                                                                                            >  2:316811/1‑190 (MQ=255)
                                                                                                                                                                         gcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTGCGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                                                                                            <  1:316811/190‑1 (MQ=255)
                                                                                                                                                                              tACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTa                                                             <  1:384731/248‑1 (MQ=255)
                                                                                                                                                                              tACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACa                                                           <  2:262543/250‑1 (MQ=255)
                                                                                                                                                                               acgaacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACaaa                                                         <  1:160046/251‑1 (MQ=255)
                                                                                                                                                                               acgaacgaCAGCAATGACCGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTTTGCCAGGCGTTTAGCGGCAGCTTTCAACCCCGCATAGGTGAATCGTACATTGGTCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACaaa                                                         <  2:671885/251‑1 (MQ=255)
                                                                                                                                                                                cgaacgaCATCAATGACAGCGTCTATGTCAGCACTCTCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGAATGAGAACTGTCTTACCTGTAAAAGCGCCCAGCAGGACTCCTTGATTTATTTTGTAACTTGCATTACaaa                                                         >  2:629055/1‑250 (MQ=255)
                                                                                                                                                                                        aTCAATGACAGCGGCTCTGTCAGCTCTAACTTGGAATACTGCTGTCGCCCTACTTGCTTGTGCCAGGCGATTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACTTTGTCCCCATCGGTGACGAAACGACGTACTATAGAGTCACCGAGACCACGACTGTCACCGAGGATGAGAACTTTCTTACCTGTAAAAGCGCTCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTCTTa                                                <  2:173632/251‑1 (MQ=255)
                                                                                                                                                                                              gACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTc                                          >  1:35420/1‑251 (MQ=255)
                                                                                                                                                                                               aCAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCt                                         >  1:223024/1‑251 (MQ=255)
                                                                                                                                                                                                  gCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACTAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTc                                      >  2:413863/1‑251 (MQ=255)
                                                                                                                                                                                                   cGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGCTACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGAACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTc                                      >  1:687998/1‑250 (MQ=255)
                                                                                                                                                                                                   cGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCa                                     <  1:333393/251‑1 (MQ=255)
                                                                                                                                                                                                        ctGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAAc                                   <  1:478240/248‑1 (MQ=255)
                                                                                                                                                                                                          gTCAGCGCTATCTGTGAATCCTGCTGTCGCTCCAGTCTCTTGGGCCAGGCGTTTTGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCACCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAATGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGtt                              <  2:274444/251‑1 (MQ=255)
                                                                                                                                                                                                             aGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACCTGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTg                           >  1:135279/1‑251 (MQ=255)
                                                                                                                                                                                                              gCACTATCTGTGAATACTGCTGTCGCTCTAGTATCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTATCATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGt                          <  2:426696/251‑1 (MQ=255)
                                                                                                                                                                                                                       gtgAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATgg                     <  1:302042/247‑1 (MQ=255)
                                                                                                                                                                                                                         gAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTa               <  1:66086/251‑1 (MQ=255)
                                                                                                                                                                                                                             aCTGCTGTCGCTCTAGTCTCTTGTGCCCGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGAGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTAATTTATTCTGTACCATGCCTTCCAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGAAATGTACACTaaaa            >  2:525934/1‑250 (MQ=255)
                                                                                                                                                                                                                             aCTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGAt                                                                                                                             <  1:449150/137‑1 (MQ=255)
                                                                                                                                                                                                                             aCTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGAt                                                                                                                             >  2:449150/1‑137 (MQ=255)
                                                                                                                                                                                                                                     cGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAtctc   <  2:223024/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTcg  >  1:50223/1‑251 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
GATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATAT‑CCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCG  >  NC_000913/4473667‑4474146

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