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,6220TC98.9% 276.7 / ‑4.8 91D496D (GAT→GACrdcAregulator of diguanylate cyclase RdcA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (38/51);  total (39/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

CCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTC  >  NC_000913/4328375‑4328857
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
ccTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGAAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGAcgcg                                                                                                                                                                                                                                          >  1:591020/1‑251 (MQ=255)
         aTGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGt                                                                                                                                                                                                                                  <  1:166648/250‑1 (MQ=255)
         aTGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGt                                                                                                                                                                                                                                  <  2:213885/250‑1 (MQ=255)
            ctctGCGATCTGACGCGCATAACCTGATGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCATGAACTCCGCCAACCACTTTCTGCGTCAACAGCAAGCCGCGCTGAAGGTGc                                                                                                                                                                                                                                >  2:437005/1‑249 (MQ=255)
             tctGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGt                                                                                                                                                                                                                             >  1:126374/1‑251 (MQ=255)
               tGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGtt                                                                                                                                                                                                                            >  2:302515/1‑250 (MQ=255)
                 cGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCCGGATTTTCGTGCGCAAGGGTTAAACGTCGCTTATGAAAAATTGCGACAAAATATCCCACAGATCGAAAAAAGTTTGCAGCTATTGCAACTCCATCAGGCGCAGGTGAGCGGAGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCACCACCCACTTTCTGCTTCAACAGCAAGACGCACTGAAGGTGCAGTTAgc                                                                                                                                                                                                                         >  1:557859/1‑251 (MQ=255)
                   aTCTGCCGCGCATAAACTGTTGCATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAgccg                                                                                                                                                                                                                       >  2:74013/1‑251 (MQ=255)
                     cTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGTGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCCACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTCgccgcc                                                                                                                                                                                                                     >  1:502765/1‑251 (MQ=255)
                       gCCGCGCTTAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGTAAAATATCCCCCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCGCCGCCAACCACTTTCCGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCtt                                                                                                                                                                                                                   <  2:522045/251‑1 (MQ=255)
                       gCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAgccgc                                                                                                                                                                                                                      <  2:561546/248‑1 (MQ=255)
                       gCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAgc                                                                                                                                                                                                                         <  1:82174/245‑1 (MQ=255)
                           cgcATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGttt                                                                                                                                                                                                               <  2:82540/251‑1 (MQ=255)
                              aTAAACTGTTGAATTACGCGCCGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAACGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTc                                                                                                                                                                                                              >  2:564627/1‑249 (MQ=255)
                              aTAAACTGTTGAATTACGCGCAGCAGGCGCGGGACTACCTGGATGTTCGTGCGCCCGGGTTAAACGTCGCTTGTGAACAATTGCGGTAAAATATCCTCCAGATCGAAGAAAGTGTGCAGCTCTTGCAACTAAATCCGGCGCAGGTGAGCGGCGAGATTATTCATGAAATCGAGCAGGCCCTGACCTCGGCCAACCACTTTCGGCGTAAAGAGCAAGACGCGGTAAATGTGCAGTTAGCCGCCTTGTTTCAg                                                                                                                                                                                                            <  2:607840/251‑1 (MQ=255)
                              aTAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTCTCGAGCGCACGGATTAAACGTCGCTTGTGAACAATTGCGGCAACATATACACCAGCTCGAAGCAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACCGCAAGACGCGCTGAAGGTGCAGTTAGCCGCATTGTTTc                                                                                                                                                                                                              >  2:589236/1‑249 (MQ=255)
                              aTAAACGGTTGAAGTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACAGGTTAAACGTCTCTTGTGAACAAGTGCGGCAAAATATCCACCAGAGCGAAGAAAGGTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGGCCTCCGCCAACCACTTTCTTCGTCAACATCAAGACGCGCTGAAGGTGCAGTTAGCCTCCTTGTTTCAg                                                                                                                                                                                                            <  2:79513/251‑1 (MQ=255)
                                  aCTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATGTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGATCTGGCCCGGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGgatg                                                                                                                                                                                                        <  2:105459/251‑1 (MQ=255)
                                     gTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGgatga                                                                                                                                                                                                       <  1:68769/249‑1 (MQ=255)
                                        gAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCgg                                                                                                                                                                                                  >  1:395579/1‑251 (MQ=255)
                                               gcgcAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGtt                                                                                                                                                                                           <  1:30148/251‑1 (MQ=255)
                                                   agcagGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTGAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGc                                                                                                                                                                                       <  1:286685/251‑1 (MQ=255)
                                                         cgcgGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGAtt                                                                                                                                                                                 <  2:459288/251‑1 (MQ=255)
                                                         cgcgGGAAGCCCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCATCCTCTTTCTGCGTCAACAGCAAGTCGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGAtt                                                                                                                                                                                 <  2:428003/251‑1 (MQ=255)
                                                          gcgGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTc                                                                                                                                                                                <  1:170705/251‑1 (MQ=255)
                                                              gAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTAc                                                                                                                                                                            >  1:548073/1‑251 (MQ=255)
                                                                  acatGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCg                                                                                                                                                                         <  2:280728/247‑1 (MQ=255)
                                                                           ttCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGac                                                                                                                                                                <  1:96837/250‑1 (MQ=255)
                                                                            tCGTGCGCACGGGGTATACAACGCTTGTGCACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAACTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGATGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTACGCGAGATTCGTACCTGCTGTGAGacac                                                                                                                                                              >  2:313805/1‑251 (MQ=255)
                                                                             cGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACt                                                                                                                                                             <  1:42677/251‑1 (MQ=255)
                                                                               tGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCAAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTCAGCGAGATTCGTACCCGCTGTGAGACACTGt                                                                                                                                                           >  1:309384/1‑251 (MQ=255)
                                                                                  gcACGGGTTAAACGTCGCTTGTGATCAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTGTGCAGCTATTTCAACGCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCGCTTTCTGCGTCAACAGCAAGACTCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGAGGATTCGGAGCCGTTAAGCGAGAGTCGGACCCGCTGTGAGACACTGTTAc                                                                                                                                                        <  1:435553/251‑1 (MQ=255)
                                                                                    aCGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAg                                                                                                                                                      <  1:218545/251‑1 (MQ=255)
                                                                                             aCGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGGGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAg                                                                                                                                             <  1:589236/251‑1 (MQ=255)
                                                                                                 cGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGCGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAac                                                                                                                                          >  1:113868/1‑250 (MQ=255)
                                                                                                   cTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACAcc                                                                                                                                       >  2:563766/1‑251 (MQ=255)
                                                                                                    ttGTGATCAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAAAAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGCCACTGTTACAGACGGCGCAGAACAcc                                                                                                                                       >  2:159536/1‑250 (MQ=255)
                                                                                                    ttGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACAcc                                                                                                                                       >  2:203527/1‑250 (MQ=255)
                                                                                                      gtgAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCAt                                                                                                                                     <  1:312681/250‑1 (MQ=255)
                                                                                                      gtgAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATc                                                                                                                                    >  1:242285/1‑251 (MQ=255)
                                                                                                        gAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGAAACTGTTACAGACGGCGCAGAACAccctc                                                                                                                                    >  2:550613/1‑246 (MQ=255)
                                                                                                            aaTTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTcgc                                                                                                                              <  1:275212/251‑1 (MQ=255)
                                                                                                             atggCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGTCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTcgcg                                                                                                                             <  1:438795/248‑1 (MQ=255)
                                                                                                             aTTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTcgcg                                                                                                                             >  2:132333/1‑251 (MQ=255)
                                                                                                               tGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGAAACTCAATCAGGCGCAGGTGAGAGGCGAGATTAAACATGCCATCGAGCTGGCCCTGACCTCCGCCAACCACATTCTGCGTCAACCGCAAGCCGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTAAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGCACACCATCAGTCGCGAc                                                                                                                           >  1:368481/1‑251 (MQ=255)
                                                                                                                       aaTATCCACCGGATCGAAGACAGTTTGTAGCTATTGCACCTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGAAAGACGCGCTGAAGGTGCAGTTCGCCGCCTTGTTTCAGGATGATTCGGCGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACCACATCAGTCGCGACTTTACGCt                                                                                                                   >  2:537160/1‑251 (MQ=255)
                                                                                                                        atatCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTCCGCTg                                                                                                                  >  1:270077/1‑251 (MQ=255)
                                                                                                                        atatCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTg                                                                                                                  <  1:462378/251‑1 (MQ=255)
                                                                                                                        atatCCACCAGAGCGAAGAAATGTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGGCAACAGCAAGACGCGCTGAAGGTGCAGGTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGTCGCAGAACACCATCAGTCGCGATTTTACGCTg                                                                                                                  <  2:319423/251‑1 (MQ=255)
                                                                                                                          atCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGc                                                                                                                 <  2:454289/250‑1 (MQ=255)
                                                                                                                               ccaGAACGAAGAAAGGTTGCAGATATTGCAACTCAATCAGGCGCAGGAGCGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGAAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGTTGATTCGGAGCCGTTAAGCGAGATTCGTAACAGCTGTGAGAAAATGTTACAGACGGCGCAGACCACCATCAGTAGGGACTTTAACCTGCGTTTTg                                                                                                           >  2:421879/1‑251 (MQ=255)
                                                                                                                                   aTCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCg                                                                                                        >  2:526790/1‑250 (MQ=255)
                                                                                                                                      gaagaaAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGCGCt                                                                                                    >  1:51270/1‑251 (MQ=255)
                                                                                                                                           aaGTCTGCAGCTATTGCAACGCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAAt                                                                                               <  2:86675/251‑1 (MQ=255)
                                                                                                                                             gttggCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATAAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGTTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGATATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACATCATCAGTCGCGACTTTACGCTGCGTGTTGCCGAGCTTGAATcc                                                                                             <  1:13883/247‑1 (MQ=255)
                                                                                                                                                 gCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCAccc                                                                                         >  1:192343/1‑251 (MQ=255)
                                                                                                                                                  cAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCt                                                                                        <  2:611038/251‑1 (MQ=255)
                                                                                                                                                   cgctagtgcatcgcAATCAGGCGCAGTTGAGCGGCGAGAGTAAACATGAAATCGAACTGGCCCTCACCTCCTACAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCATTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTCCAGTCGGCGCAGAACACCATCATTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCtt                                                                                       <  2:392647/238‑1 (MQ=255)
                                                                                                                                                       aTTGCAACTAAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTATAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTg                                                                                     <  1:419884/249‑1 (MQ=255)
                                                                                                                                                         tGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGCAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCgg                                                                                 >  2:613451/1‑251 (MQ=255)
                                                                                                                                                           caaatCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTTGCCCTGACCTCCGCCAACCACGGTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTGGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGGGAGACACTGTGACAGTCGGCGCAGAACACCATCAGTTGCTACTTTACGCTGCGTTTTGCTGAGCGTGAATCCACCTTTTTCCGGgt                                                                               <  2:592828/247‑1 (MQ=255)
                                                                                                                                                           cAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGgt                                                                               >  1:8340/1‑251 (MQ=255)
                                                                                                                                                           cAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGgt                                                                               <  2:333299/251‑1 (MQ=255)
                                                                                                                                                                cAATAAGGCGCAGGTGATCGGCGAGATTAAACAAGAAATCGAGCTGGCCCTGCCCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTGTCCGCTGTGAGACACTGTTACAGACGGCGCAGAGCACCATCAGTCGCGACGGTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTaa                                                                          <  1:250281/251‑1 (MQ=255)
                                                                                                                                                                  agcAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAAcc                                                                        <  2:95454/249‑1 (MQ=255)
                                                                                                                                                                     aGGCGCCGGTCAACGGCGCGATTAAGCATGAATTCGAGCTGGCCCTGACATCCGCCAACCACTTTCTGCGTCAACAGCAAGAAGCGCTGAAGGGGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGCGCCGTTAAGCGAGATTCGTACCCGCTGTGAGCCACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGAt                                                                     >  2:597665/1‑251 (MQ=255)
                                                                                                                                                                        cgccggTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGtt                                                                  <  2:369249/247‑1 (MQ=255)
                                                                                                                                                                         gcAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGTGTGTTAACCGATGTTa                                                                 <  1:550880/251‑1 (MQ=255)
                                                                                                                                                                          cAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTGTGCCGGGTGTTAACCGATGTTAt                                                                <  1:148916/251‑1 (MQ=255)
                                                                                                                                                                           aGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTAtt                                                               <  2:502145/251‑1 (MQ=255)
                                                                                                                                                                            ggagAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCCGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGGGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGTCTTTACGCTGCGTTTTGCCGAGCTTGACTCCACCCTTTGCCGGGTGTTAACCGATGTTATTc                                                              >  2:3135/4‑251 (MQ=255)
                                                                                                                                                                             gTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCg                                                             <  2:126374/251‑1 (MQ=255)
                                                                                                                                                                                   agcgAGATTAAATATGAAATCGTGCTGGCCCCGACCTCCGCCAACCACTTTATGCGTCAACAGCAAGCCGCGCTGCAGGTGCAGCTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGAACCCGCTGTGAGACACTGTTCCAGACTGCGCAGAAAACCTTAAGCAGCGCCTTTACGCTGCGTTTTGCCGAACTTGAATCCACCCTTTGCCGGGGGTTATCCGGTGTTATCCGCCCCAt                                                       >  2:347454/2‑251 (MQ=255)
                                                                                                                                                                                     cgAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTg                                                     <  1:489764/251‑1 (MQ=255)
                                                                                                                                                                                       agaTTAAACATGAAATCGAGCTGGCCCTGAACTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGCGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTg                                                     >  1:537572/1‑249 (MQ=255)
                                                                                                                                                                                               aaTGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGt                                           <  1:277803/250‑1 (MQ=255)
                                                                                                                                                                                                        gAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGGGCAACAAGTCAAAATgg                                   >  2:601957/1‑250 (MQ=255)
                                                                                                                                                                                                             ggCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTg                              <  1:293248/250‑1 (MQ=255)
                                                                                                                                                                                                              gCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTg                              >  1:265841/1‑249 (MQ=255)
                                                                                                                                                                                                              gCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAg                            <  2:520818/251‑1 (MQ=255)
                                                                                                                                                                                                               cccggACCTCCGCCAACCCCTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCTTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGc                           <  1:277144/247‑1 (MQ=255)
                                                                                                                                                                                                                ccTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCg                          <  1:214020/251‑1 (MQ=255)
                                                                                                                                                                                                                   gACCTCCGCCAACCGCTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGAGTCGTACCCGCTGTGATACACTGTGACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTAGTCGCCCCATGGAGCAACAAGTCAAAATGGAATTGAGCGa                         <  2:569716/249‑1 (MQ=255)
                                                                                                                                                                                                                         cGCCAACCACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGTTGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTAAGGGt                 >  1:34065/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccAACCTCTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGttt               <  2:450951/251‑1 (MQ=255)
                                                                                                                                                                                                                               ccACTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAAGTGAGCGAGTCAGGGTTTCGTc           <  1:441263/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGCGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACACGTCAAAATGGAATTGAGCGAGTCAGGGttt               >  2:554445/1‑244 (MQ=255)
                                                                                                                                                                                                                                  cTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGttt               <  1:554445/244‑1 (MQ=255)
                                                                                                                                                                                                                                  cTTTCTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTATTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAACATGGAATTGAGCGAGTCAGGGTTTCGTCCTg        >  2:150452/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cTGCGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGtt    >  1:438317/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gCGTCAACAGCAAGACACGCTGAAGGTGCAGGTAGCCGCCTTGTTTCAGTAGGATTCGGGGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCTCAGAACACCATCAGTCGCGACTTGACTCAGCGTTTTGCTTAGCTCGAATCCACCCTTTGCCGGGTGTGAACTGATGGTATTCGCCCCTGTGATCAACAAGGCAAAATGGAAATGAGCGAGTCAGGGTTTCGTTCTTGGTTTc  <  2:269776/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cGTCAACAGCAAGACGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTc  >  1:151074/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
CCTCGTTGTATGCTCTGCGATCTGCCGCGCATAAACTGTTGAATTACGCGCAGCAGGCGCGGGAATACCTGGATTTTCGTGCGCACGGGTTAAACGTCGCTTGTGAACAATTGCGGCAAAATATCCACCAGATCGAAGAAAGTTTGCAGCTATTGCAACTCAATCAGGCGCAGGTGAGCGGCGAGATTAAACATGAAATCGAGCTGGCCCTGACCTCCGCCAACCACTTTCTGCGTCAACAGCAAGATGCGCTGAAGGTGCAGTTAGCCGCCTTGTTTCAGGATGATTCGGAGCCGTTAAGCGAGATTCGTACCCGCTGTGAGACACTGTTACAGACGGCGCAGAACACCATCAGTCGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTC  >  NC_000913/4328375‑4328857

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