Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,328,622 T→C D496D (GAT→GAC)  rdcA → regulator of diguanylate cyclase RdcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,328,6220TC97.5% 222.8 / ‑2.4 82D496D (GAT→GAC) rdcAregulator of diguanylate cyclase RdcA
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (40/38);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

TTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTT  >  NC_000913/4328380‑4328868
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
ttGTATGCTCTGCGCTCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCACAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCACTGACCTCCGCCAACCACTTTCTGCGTCAACCGCAAGACGCGCTGa                                                                                                                                                                                                                                                   >  2:616674/1‑250 (MQ=255)
  gTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAg                                                                                                                                                                                                                                                 <  1:344407/250‑1 (MQ=255)
    attctctGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTg                                                                                                                                                                                                                                              >  1:585576/4‑251 (MQ=255)
        gcggCGATCTGCCGCGCATAAACTGTTGAATTACGCGAAGCAGGCGCGGGAAGTCCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTTCGCCAACCACTTTCGGCGTCAACAGCAAGACGCGCTGAAGGTGCAGt                                                                                                                                                                                                                                          <  2:197735/248‑1 (MQ=255)
                  gCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAgcc                                                                                                                                                                                                                                     <  1:267546/246‑1 (MQ=255)
                  gCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGAGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAgcc                                                                                                                                                                                                                                     >  2:267546/1‑246 (MQ=255)
                             acggTTGAATTCCGCGCAGCAGGCGCTGGAATTCCTGGATTTTCGTGCGCACGAGTTAAACGTCGCGTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGgata                                                                                                                                                                                                                     <  2:532599/248‑2 (MQ=255)
                               tGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGAGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGGGCAGTTAGCCGCCTTGTTTCAGgatgat                                                                                                                                                                                                                   >  1:52887/1‑251 (MQ=255)
                                    aaTTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGACGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCggg                                                                                                                                                                                                              >  2:194443/1‑250 (MQ=255)
                                         cgcgCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGt                                                                                                                                                                                                         >  1:664390/1‑251 (MQ=255)
                                             cagcagGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGCCCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAg                                                                                                                                                                                                     <  1:519813/251‑1 (MQ=255)
                                             cagcagGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAg                                                                                                                                                                                                     <  1:7579/251‑1 (MQ=255)
                                              agcagGCGCGGGACTACCTGGATTTTCGTGCGCACGGGATAACCGTCGCTTGTGAACAATTGAGGCAAAATATCCACCAGATCGACGAAAGTTTGCAGCTATTGCAACTCAATAAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTCCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCAGCCTTGTTTCAGGCTGATTCGGGGCCGTTAAGc                                                                                                                                                                                                    >  2:323036/1‑251 (MQ=255)
                                              agcagGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCTAAGAAAGTTTTCAGCTATTGCCACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCATGACCTCCGCCAACACCTTTCTGCGACAACAGCAAGCCGCGCAGAAGGGGCAGTTAGCCGCCTTGTTTCACGAAGCTTCGGAGCCGTTAAGc                                                                                                                                                                                                    >  2:79502/1‑251 (MQ=255)
                                               gcagGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCg                                                                                                                                                                                                   >  1:609584/1‑251 (MQ=255)
                                                      tgggAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGATAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGGGA‑TTCg                                                                                                                                                                                           >  2:68018/2‑251 (MQ=255)
                                                        ggAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATAAAATCGAGCTGGCACTGACCTCCGCCAACCACTTTCGGCGTCACCAGACAGACGCGCTGAAGGGGCAGTTAGCCGACTTGTTTCAGGATTATTCGCAGCCGTTACGCGAGC‑TTCGt                                                                                                                                                                                          >  2:28349/1‑250 (MQ=255)
                                                         gAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGCTGATTCGGAGCCGTTAAGCGAGA‑TTCGTAc                                                                                                                                                                                        >  1:636143/1‑251 (MQ=255)
                                                          aaTACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTAcc                                                                                                                                                                                       >  1:110805/1‑251 (MQ=255)
                                                            tACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGCGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTAAGCGGCGAAATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCCGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTAcacg                                                                                                                                                                                     >  1:248196/1‑248 (MQ=255)
                                                               cTGGATTTTCGTGCGCACGGGTTAAACGTCGATTGTGAACAATTGCGGAAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGAGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAAACACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGCGCCGTTAGCCGCCTTGTTTCAGGATGATTCGGATCCGTTAAGCGAGA‑TTCGTCACCGCtg                                                                                                                                                                                  >  2:403822/1‑251 (MQ=255)
                                                                    ttttCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGC‑TTCGTACCCGCtgtg                                                                                                                                                                                >  1:311217/1‑248 (MQ=255)
                                                                     tttCGTTCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTgaga                                                                                                                                                                             <  2:315852/250‑1 (MQ=255)
                                                                         gTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTACGCGAGA‑TTCGTACCCGCTGTGAGACACTg                                                                                                                                                                        >  1:102511/1‑251 (MQ=255)
                                                                         gTGCGCACGGGTTAAACGTCGCTTGTGAACAAGTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTTCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGTCCGCCGCCAACCACGTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTg                                                                                                                                                                        <  2:565161/251‑1 (MQ=255)
                                                                           ccgcACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGtt                                                                                                                                                                      <  2:329751/250‑1 (MQ=255)
                                                                            cgcACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGtt                                                                                                                                                                      >  1:329751/1‑250 (MQ=255)
                                                                            cgcACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGtt                                                                                                                                                                      >  1:63285/1‑250 (MQ=255)
                                                                              aaCGGGGTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCCCCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAATGTGCAGTTAGCCGCCTTGTTTCAGGATTATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCTCTGTGAGACACTGTTACa                                                                                                                                                                   <  2:631108/250‑1 (MQ=255)
                                                                               aCGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGAACTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCG‑GACTTCGTACCCGCTGTGAGCCACTGTTAc                                                                                                                                                                    >  2:169710/1‑249 (MQ=255)
                                                                                    ttAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGgc                                                                                                                                                             >  1:492829/1‑251 (MQ=255)
                                                                                          gTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACATCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGATGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTCCCCGCTGTGTGACACTGTTACAGACGgcgccga                                                                                                                                                        >  2:503308/1‑247 (MQ=255)
                                                                                                    aaCAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑Gt                                                                                                                                            >  1:60180/1‑251 (MQ=255)
                                                                                                     aCAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGACATCGAGCTGGACCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGCCGCTCTGAAGGGGCCGTTAGCCGCCTTGTTTCTGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGAAACTGTTACAGACGGCGCAGAAAACCATCA‑Gt                                                                                                                                            >  2:198954/1‑250 (MQ=255)
                                                                                                     aCAATTGCGGCAAAATATCCACCAGATAGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCATCAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGCGCCGTTACGCGAGA‑TTCGTACCAGCTGTGCGACACTGTTACAGACGGCGCAGCACACCATCA‑GTc                                                                                                                                           >  2:649979/1‑251 (MQ=255)
                                                                                                     aCAATTCCGGCAAAATATCCCCCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGTTGAGCGGCCAGATTAAACATTAAATCGAGCTGGCCCTGACCTCCGCCAACCCCTTTCTTCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCTCAGTACACCATCA‑GTc                                                                                                                                           <  2:63285/251‑1 (MQ=255)
                                                                                                        aTTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCACTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTcgcg                                                                                                                                        >  2:137939/1‑251 (MQ=255)
                                                                                                                            aGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCGGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGcc                                                                                                                    >  1:433486/1‑251 (MQ=255)
                                                                                                                                        gTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATcc                                                                                                        <  1:68018/251‑1 (MQ=255)
                                                                                                                                         tttGCAGCTATTGCAACTCAATCAGGCGCAGTTGAGCGGCGAGTTTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCCCTTTCTGCGTCAACAGAAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTAACCGCAGTGAGCAACTGTTACAGACGGCGCAGGACACCATCA‑GTCCGGGCTTTACGCTGCGtttggcgtagctttaatacc                                                                                                       >  2:250053/1‑234 (MQ=255)
                                                                                                                                            gCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCAccc                                                                                                    >  1:350200/1‑251 (MQ=255)
                                                                                                                                                  aTTGCAACGCAATCCGGTGCAGGTGAGCGGCGAGAGTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGcc                                                                                              <  2:664390/251‑1 (MQ=255)
                                                                                                                                                           cAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCGCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTaa                                                                                     <  1:610984/251‑1 (MQ=255)
                                                                                                                                                              tCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCGGAACACCATCA‑GTCGCGCCTTTACGCTGCGTTTTGCCGAGACTGACTCCCCCCTTTGCCGGGTGTTAACCg                                                                                  >  2:329880/1‑251 (MQ=255)
                                                                                                                                                               cAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCg                                                                                  <  1:544521/250‑1 (MQ=255)
                                                                                                                                                                aGGCGCAGGTTAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCTCTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACTCCATCA‑TTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGAt                                                                                <  1:323036/251‑1 (MQ=255)
                                                                                                                                                                 ggCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTGTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATg                                                                               <  1:391921/251‑1 (MQ=255)
                                                                                                                                                                 ggCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACATCCGCCAACCACATTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑GTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATg                                                                               <  2:636143/251‑1 (MQ=255)
                                                                                                                                                                  gcgcAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGt                                                                              <  1:169710/251‑1 (MQ=255)
                                                                                                                                                                    gcAGGTGAGCGTCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCTGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGATTGAATCCACCCTTTGCCGGGTGTTAAccactgttt                                                                            >  2:181241/1‑244 (MQ=255)
                                                                                                                                                                    gcAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGACTTGTTTCAGGATGATTCGGAGACGTTAAGCGAGA‑TTCATACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGATGCGATTTGACGAGCTTGAGTCCACCCTTTGCCGGGTGTTAACCgcagttt                                                                            >  2:123533/1‑245 (MQ=255)
                                                                                                                                                                           aGCGGCGAGATTAAACATGGAATAGAGCTGGCCATGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGCCTTTACGCTGCGTTTTTCCGGGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGccc                                                                     >  2:671070/1‑251 (MQ=255)
                                                                                                                                                                           aGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGccc                                                                     <  1:616674/251‑1 (MQ=255)
                                                                                                                                                                                 gaggttAACCATGAAATCGGGCTGGCCCTGACCTCCGCCACCCTCTTTCGTCTTCAAACGCAACACGCTCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGTATGATTCGGAGCCGTTAAGCGAGC‑TTCGTACCCGCTGTGAGACGCTGTTACAGACGGCGCAGAACCCTATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTTTTCACCGATGTCATTCGCCCCATTGa                                                               <  2:79245/247‑1 (MQ=255)
                                                                                                                                                                                  agaTTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTg                                                                >  1:66615/1‑249 (MQ=255)
                                                                                                                                                                                      tAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGcaac                                                          >  1:288360/1‑251 (MQ=255)
                                                                                                                                                                                       aaaCATGAAATCGAGCTGTCTCTGACCCCCGCCAACCACGTGTTGCGTCCACAGCAAGCAGCGGTGAAGGTGCCGTTAGCCGCCTTGTTGCAGGATGAGTCGGAGCCGGTAAGCGAGA‑TTCGGACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACTCCATCATGT‑GCGACTTTACGCTGCGTTTTGCCGATCGTGATTCCACCCTTTGCCGGTTTTTAACCGATGTTATTCGCCCCTTTGAGcaac                                                          <  2:277835/250‑1 (MQ=255)
                                                                                                                                                                                              aaaTCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCg                                                                                  <  1:27510/219‑1 (MQ=255)
                                                                                                                                                                                              aaaTCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCAGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCaaa                                                  <  1:142916/251‑1 (MQ=255)
                                                                                                                                                                                                  cgggCTGGCCCTGACCTCCGCCAACCCCTTTCGGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCTGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑GTCGTACCCGCTGTGAGACACTGTTACAGTCGGCGCAGAACACCATCA‑GTCGCGACGTTACGCTGCGTTTTGCCGATCTTGAATCCACCCTTTGCCGGGTGGTAACCGATGTTATTCGCCCCATCGAGCAACAAGTCAAAATgg                                              <  2:183468/248‑1 (MQ=255)
                                                                                                                                                                                                              gACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCTTTAAGCGATA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAAGGGAATTGAGCGAg                                   <  1:79502/250‑1 (MQ=255)
                                                                                                                                                                                                                 cTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTc                                 >  2:351993/1‑249 (MQ=255)
                                                                                                                                                                                                                 cTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGCTGTTATTCGCCCCATTGAGCCACAAGTCAAAATGGAATTGAGCGGGTc                                 >  2:290766/1‑249 (MQ=255)
                                                                                                                                                                                                                    cGCCAACCACTTGCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTGCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGGTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGt                            <  1:656814/251‑1 (MQ=255)
                                                                                                                                                                                                                       cAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTc                         >  1:238410/1‑251 (MQ=255)
                                                                                                                                                                                                                        acacACTTTCTGCGTCACAAGCAAGACCCCCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTTCCCGCTGTGAGACACTGTTACAGACGGCGCAGAACCCCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCg                        <  2:465335/248‑1 (MQ=255)
                                                                                                                                                                                                                        aaCCACTTTCTGCGTCAACAGCAAGACGCGCTGTAGGTGCAGTTAGCCGCATTGTTTCAGGATGATTCGGAGACGTTAAGCGATA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGCACACCATCA‑GTCGCGAATTTACGCTGCGATTTGCCAAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCTACAAGTCAATATGGAATTGAGCgggtatggttttcg                        >  2:243952/1‑238 (MQ=255)
                                                                                                                                                                                                                        aaCCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGGACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTGAACCGATGTTATTCGCCCCATTGAGCAGCAAGTCAAAATGGAATTGAGCGAGTCAGGATTTct                        >  1:201246/1‑250 (MQ=255)
                                                                                                                                                                                                                         aCCACTTTATGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGt                       <  2:337862/251‑1 (MQ=255)
                                                                                                                                                                                                                            ccgttttGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAATCGAGA‑TTCGTACCCTCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCt                    <  2:102511/245‑1 (MQ=255)
                                                                                                                                                                                                                            ccTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGAAGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCt                    <  1:14495/250‑1 (MQ=255)
                                                                                                                                                                                                                               ttCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTgg                  <  1:351993/250‑1 (MQ=255)
                                                                                                                                                                                                                               ttCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTgg                  <  1:324753/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTc             <  1:403822/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACTCTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCa            <  1:290766/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCa            <  1:291873/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cAACAGTAAGAAGCGCTGTAGCTGCCGTTAACCGCATTGTTTCAGGATGCTTCGCAGCCGTTAAGCGCGA‑TTCGTACCCGCTGTGAGACACTGTTACCGACGGCGCAGAACACCCTCA‑GTCGCGCCTTTAAGCTGCGTTTTGCCGAGCTTGAATCCACTCTTTGCCGGGTGTTAAACGCTGTTATTCGCCCCATTGCGCATCAACTCAACATGGAATTGAGCGAGTCAGGACTTCGTCCTGGGTTTCCttt         >  2:286313/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aCAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGcc                                                                      >  1:262117/1‑188 (MQ=255)
                                                                                                                                                                                                                                         aCAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGcc                                                                      <  2:262117/188‑1 (MQ=255)
                                                                                                                                                                                                                                         aCAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTc       <  1:319196/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGAGGTTTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTga   <  2:508581/251‑2 (MQ=255)
                                                                                                                                                                                                                                              aaGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGtt  <  1:201659/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
TTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGA‑TTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCA‑GTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTTTCCTGTT  >  NC_000913/4328380‑4328868

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