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

GATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCG  >  NC_000913/4473667‑4474146
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
gatTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGGGCCAGGc                                                                                                                                                                                                                                       >  1:222746/1‑251 (MQ=255)
gactcaGATGCACCCGCCTTCGGGCATCTGCCGGGCGGCTTCTACAGAGGCATGCTAAGGCGCATGAATATTGATTTTGTAAAGGCGATCAATATCGTCGGCCTTTAATTACAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACAAGGATATCCAATGCTCCGCCCTTACTACAGACATCAATGACCGCGTCCCAGTCAGCCCTATCTGTGAATACTGCTGTCTCACTAGTCTCATGTGCCAGGc                                                                                                                                                                                                                                       >  2:434214/6‑251 (MQ=255)
         gcgcCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCg                                                                                                                                                                                                                               <  1:502558/250‑1 (MQ=255)
         gcgcCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGGCGCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCCCTTGTGCCAGGCGTTTAGCg                                                                                                                                                                                                                               <  2:632205/250‑1 (MQ=255)
               gccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGACCGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATCCTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCtt                                                                                                                                                                                                                        >  2:691981/1‑251 (MQ=255)
                  ttCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCg                                                                                                                                                                                                                     >  1:235262/1‑251 (MQ=255)
                      ggCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGCGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGCCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCACTTGAGCCAGGCGGTTAGTGGCATCTTTCGAccc                                                                                                                                                                                                                 >  2:284310/1‑251 (MQ=255)
                           cTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCAt                                                                                                                                                                                                            >  2:302042/1‑251 (MQ=255)
                                 ggcggcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAAGCGATCAATATAGTCGGCATTTCATTCACGGGCCTCGCCACAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGACGCTATTAAGAACCACATCAATGACAGCGTCTCTGTCAGCACTCTCTGTGAATACTGCCGTCGCTCTCGTCTCTTTTGCCCGGCGTTTAGAGGCATCTTTCGAACCCGCATAGGTgg                                                                                                                                                                                                      >  2:410435/1‑250 (MQ=255)
                                     gcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCg                                                                                                                                                                                                  >  1:547210/1‑251 (MQ=255)
                                         cAACAGAGGCATGATAAGGAGAATTAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                                                                                                                                               <  2:193948/250‑1 (MQ=255)
                                                   aTGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAt                                                                                                                                                                                    >  1:598810/1‑251 (MQ=255)
                                                         aGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGa                                                                                                                                                                              <  1:519847/251‑1 (MQ=255)
                                                          gtaGCATGAATATTGATTATGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGAc                                                                                                                                                                             <  1:410435/249‑1 (MQ=255)
                                                             gCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCCCTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGACAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCgggcacca                                                                                                                                                                           >  2:460876/1‑244 (MQ=255)
                                                             gCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCa                                                                                                                                                                                                                                          <  2:170366/187‑1 (MQ=255)
                                                             gCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCa                                                                                                                                                                                                                                          >  1:170366/1‑187 (MQ=255)
                                                             gCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGa                                                                                                                                                                           >  1:541409/1‑250 (MQ=255)
                                                                tGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCACTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGCAacg                                                                                                                                                                       >  1:481670/1‑251 (MQ=255)
                                                                                aaaGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATCGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGCCGTACGATAGCGGc                                                                                                                                                        >  2:546763/1‑250 (MQ=255)
                                                                                aaaGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCa                                                                                                                                                       <  1:558809/251‑1 (MQ=255)
                                                                                  aGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAAAGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGAACATTGGCCCCATCGGTGACGAAACGACGTCCGATAGCGGCAcc                                                                                                                                                     >  2:366518/1‑251 (MQ=255)
                                                                                         aaaTATCGTAGGCATTTAAATCCAGGGCCTCGCCAAAGACGCCAATACCCGCATTTACCTCCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCa                                                                                                                                              <  2:598810/250‑1 (MQ=255)
                                                                                            tatCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGACTACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCCTCGGTGACGAAACGGCGGACGATAGCGGccacgctaccacc                                                                                                                                            >  2:129344/1‑238 (MQ=255)
                                                                                                gtcgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGCCTGc                                                                                                                                       >  2:132689/1‑251 (MQ=255)
                                                                                                  cgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGcc                                                                                                                                      <  1:698206/250‑1 (MQ=255)
                                                                                                  cgGCATTTAATTCCAGGCCCTCGCAAAAGACGTCAATACCTGCATTTACCACTCGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGCCAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGAATGcc                                                                                                                                      >  2:333393/1‑250 (MQ=255)
                                                                                                         tAATTCCAGGGCCCCGACAATGACGCCAATACCTGCATTTACCACCATGTGTACCTATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATATCGGCACCGATACCACGACTGCCACCGAgg                                                                                                                               <  2:524588/250‑1 (MQ=255)
                                                                                                           aTTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                                                                                            <  1:101866/251‑1 (MQ=255)
                                                                                                              ccAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                                                                                            >  1:152238/1‑248 (MQ=255)
                                                                                                               cAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATgaga                                                                                                                         >  1:481792/1‑250 (MQ=255)
                                                                                                                     cTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCCTAGGTGAATCGTACATTGGCCCCATCGGGGCCGAAACGACGTACGCTAGCGGCCCCGATACCACGACTGCCACCGAGGATGCGAACTGTCt                                                                                                                  >  2:200977/1‑251 (MQ=255)
                                                                                                                         ccAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATTCGCCGCTCTTACGAACCACATCAAGGACAGCGTCTCTGTCCGCACTAGCTGTGAATACTGCTGTCGCTCTAGTCCCTTGTGCCAGGCGTTTAGCTGCATCTTTCGACCCCGCATAGGTTAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAACACGACTGCCACCGAGGATGAGAACTGTCttgcc                                                                                                              <  2:237924/251‑4 (MQ=255)
                                                                                                                          cAAAGACGCCAATACCTGCATTTACCTCCAGGATATCCAATTCGCCGCTATTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGGTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGTGGACGAAACGACGTACGATAGCTGCACCGATACCTCGACTGCCACCGAGGATGAGAACTGTCTTAcgt                                                                                                             <  2:622892/251‑3 (MQ=255)
                                                                                                                          cAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAcc                                                                                                              >  1:535129/1‑250 (MQ=255)
                                                                                                                          cAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAcc                                                                                                              <  2:535129/250‑1 (MQ=255)
                                                                                                                          cAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                                                                                             >  1:121825/1‑251 (MQ=255)
                                                                                                                             aGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACCGTCTTACCTGt                                                                                                           >  2:360703/1‑250 (MQ=255)
                                                                                                                                cGCCAATACCTGCATTTCCCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGGCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTaaaa                                                                                                       <  2:541409/251‑1 (MQ=255)
                                                                                                                                           gCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATaa                                                                                             <  1:132689/250‑1 (MQ=255)
                                                                                                                                           gCATTTACCACCAGGATATACAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATa                                                                                              <  1:502765/249‑1 (MQ=255)
                                                                                                                                            cATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGCACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTCAAAGCGCCCATAAgg                                                                                           >  2:708649/1‑251 (MQ=255)
                                                                                                                                                tACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAc                                                                                                                                                <  1:331097/194‑1 (MQ=255)
                                                                                                                                                tACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAc                                                                                                                                                >  2:331097/1‑194 (MQ=255)
                                                                                                                                                tACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTc                                                                                       >  1:677881/1‑251 (MQ=255)
                                                                                                                                                  ccaccaGGATTTCCAATGCGCCGCTCTAACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCt                                                                                     <  2:481792/251‑1 (MQ=255)
                                                                                                                                                   ccccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACCGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTg                                                                                                            <  2:377044/225‑1 (MQ=255)
                                                                                                                                                   caccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTg                                                                                                            >  1:377044/1‑227 (MQ=255)
                                                                                                                                                   caccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCtt                                                                                    <  1:298989/251‑1 (MQ=255)
                                                                                                                                                      caGGATAGCCATTGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCTGCACGTTCTGTGAATACTGCTGTCGCTATAGTCGCTTGTGACAGGCGTTTAGCGGCATCTTCAGACCCCGCATAGGTGAATCGTAAATTGGCCCCATCGGTGACGAAACGAAGTACGATAGCGGCACCGATACCACGACTGCCTCCTAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACACCTTGAt                                                                                 <  2:17544/251‑1 (MQ=255)
                                                                                                                                                        ggATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAttt                                                                               <  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)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      
GATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCTCG  >  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: