Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 587,207 G→A G525G (GGC→GGT)  unknown ← DNA gyrase subunit A (EC 5.99.1.3)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1587,2070GA98.7% 186.0 / NA 74G525G (GGC→GGT) unknownDNA gyrase subunit A (EC 5.99.1.3)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (41/32);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTGATATACTTTCAGCCAGTATAAACGTCCGCGGCTGGAGAAGCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCGCCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATGGCTTGTGCTTGCACTTCGGTTAAATAATAT  >  1/586965‑587444
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
aGTTGATATACTTTCAGCCAGCATAAACGTCCGCGGCTGGAGAAGCAGAGAATGGCATCGTGGGTGTAGGCAACCAGAAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGCAGTCGGTAATCGGTTGATATTTCACATAACTTCCGCGGGAAAGGAACACCACCACATCTTCTTGCGCAATCAAATCTTCCAGATTTATATCACCGGAAg                                                                                                                                                                                                                                        >  1:669905/1‑250 (MQ=255)
    gATATACTTTCAGCCAGTATAAACGTCCGCGGCTGGAGAAGCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAAACTTCCAGATTTATATCACCGGAAgccgc                                                                                                                                                                                                                                    >  1:696472/1‑250 (MQ=255)
                  cAGTATAAACGTCCGCGGCTGGAGAAGCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTc                                                                                                                                                                                                                      <  1:561598/250‑1 (MQ=255)
                    gTATAAACGTCCGCGGCTGGAGAAGCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTg                                                                                                                                                                                                                    >  1:645339/1‑250 (MQ=255)
                    gTATAAACGTCCGCGGCTGGAGAAGCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTg                                                                                                                                                                                                                    >  1:708253/1‑250 (MQ=255)
                      ataAACGTCCGCGGCTGGAGAAGCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTgcg                                                                                                                                                                                                                  >  1:870643/1‑250 (MQ=255)
                                   gCTGGAGAAGCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAAt                                                                                                                                                                                                     <  1:828172/250‑1 (MQ=255)
                                       gagaAGCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCt                                                                                                                                                                                                 >  1:722355/1‑250 (MQ=255)
                                            gCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGa                                                                                                                                                                                           >  1:487949/1‑251 (MQ=255)
                                            gCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGa                                                                                                                                                                                           >  1:267455/1‑251 (MQ=255)
                                                 gaATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTg                                                                                                                                                                                       >  1:843051/1‑250 (MQ=255)
                                                   aTGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTaa                                                                                                                                                                                                       <  1:549155/232‑1 (MQ=255)
                                                   aTGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTc                                                                                                                                                                                    <  1:464859/251‑1 (MQ=255)
                                                                         ccAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTc                                                                                                                                                               >  1:784339/1‑250 (MQ=255)
                                                                             tAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTa                                                                                                                                                          <  1:863296/251‑1 (MQ=255)
                                                                                       aTAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAAt                                                                                                                                                                                                     <  1:40975/198‑1 (MQ=255)
                                                                                         aaaaTCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGc                                                                                                                                              <  1:254548/251‑1 (MQ=255)
                                                                                         aaaaTCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGc                                                                                                                                              <  1:829388/251‑1 (MQ=255)
                                                                                         aaaaTCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTAGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGc                                                                                                                                              <  1:196470/251‑1 (MQ=255)
                                                                                                tcttcttTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCaa                                                                                                                                       <  1:673988/251‑1 (MQ=255)
                                                                                                    cttTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATg                                                                                                                                   >  1:315241/1‑251 (MQ=255)
                                                                                                        cATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACa                                                                                                                               <  1:258721/251‑1 (MQ=255)
                                                                                                        cATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACa                                                                                                                               >  1:305759/1‑251 (MQ=255)
                                                                                                            ttGGTTGCGGATTTTCCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCGCCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACaataa                                                                                                                           <  1:201919/251‑1 (MQ=255)
                                                                                                            ttGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGTATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACaataa                                                                                                                           <  1:37502/251‑1 (MQ=255)
                                                                                                            ttGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACaataa                                                                                                                           <  1:391994/251‑1 (MQ=255)
                                                                                                            ttGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACaataa                                                                                                                           <  1:520009/251‑1 (MQ=255)
                                                                                                              ggTTGTGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAAttc                                                                                                                                                                              >  1:423677/1‑198 (MQ=255)
                                                                                                                 tGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTAc                                                                                                                      <  1:418304/251‑1 (MQ=255)
                                                                                                                    ggATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGc                                                                                                                                         >  1:771128/1‑229 (MQ=255)
                                                                                                                     gATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCg                                                                                                                    <  1:605169/249‑1 (MQ=255)
                                                                                                                     gATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAACTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGAt                                                                                                                  >  1:802311/1‑251 (MQ=255)
                                                                                                                          acctttaccACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCGGGGCTACGCAAGATGTACaataa                                                                                                                           <  1:419512/237‑1 (MQ=255)
                                                                                                                          acctttaccACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCAt                                                                                                             >  1:330504/1‑251 (MQ=255)
                                                                                                                               taccACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCaat                                                                                                        >  1:35313/1‑251 (MQ=255)
                                                                                                                                 ccaccaCGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACg                                                                                                                                                      >  1:850047/1‑203 (MQ=255)
                                                                                                                                    ccacGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAAttt                                                                                                   >  1:892357/1‑251 (MQ=255)
                                                                                                                                      acgacgCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAAttttt                                                                                                 >  1:98270/1‑251 (MQ=255)
                                                                                                                                              gtgCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACaataa                                                                                                                           >  1:884036/1‑217 (MQ=255)
                                                                                                                                              gtgCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCg                                                                                                                    <  1:274329/224‑1 (MQ=255)
                                                                                                                                                gCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAACGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGCCATCCAATAATTTTTTGTATTCGt                                                                                        >  1:728871/1‑250 (MQ=255)
                                                                                                                                                gCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGAt                                                                                                                                    >  1:598016/1‑206 (MQ=255)
                                                                                                                                                gCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACaataat                                                                                                                          <  1:598811/216‑1 (MQ=255)
                                                                                                                                                gCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGt                                                                                        >  1:846751/1‑250 (MQ=255)
                                                                                                                                                        gTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAAttc                                                                                                                                                                              >  1:412322/1‑156 (MQ=255)
                                                                                                                                                        gTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCACCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACaat                                                                                                                             >  1:204713/1‑205 (MQ=255)
                                                                                                                                                                       atatTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAAt                                                                >  1:97889/1‑251 (MQ=255)
                                                                                                                                                                             cacaTAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTa                                                          >  1:368204/1‑251 (MQ=255)
                                                                                                                                                                              acaTAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACaataa                                                                                                                           <  1:364020/185‑1 (MQ=255)
                                                                                                                                                                              acaTAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCaa                                                                                                         >  1:211429/1‑203 (MQ=255)
                                                                                                                                                                                   aCCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAAttttt                                                    >  1:798548/1‑251 (MQ=255)
                                                                                                                                                                                    ccTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTg                                                   >  1:820551/1‑251 (MQ=255)
                                                                                                                                                                                            ggAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGt                                            <  1:727974/250‑1 (MQ=255)
                                                                                                                                                                                             gAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTAc                                                                                                                      >  1:348738/1‑175 (MQ=255)
                                                                                                                                                                                             gAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTaa                                          >  1:818608/1‑251 (MQ=255)
                                                                                                                                                                                                  ggTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTa                                     >  1:772947/1‑251 (MQ=255)
                                                                                                                                                                                                  ggTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTa                                     >  1:270327/1‑251 (MQ=255)
                                                                                                                                                                                                     caccaccacATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTaaaa                                  >  1:689472/1‑251 (MQ=255)
                                                                                                                                                                                                       ccaccacATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTgg                                                             <  1:774704/222‑1 (MQ=255)
                                                                                                                                                                                                       ccaccacATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATg                                >  1:488860/1‑251 (MQ=255)
                                                                                                                                                                                                          ccacATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATgg                               >  1:752976/1‑249 (MQ=255)
                                                                                                                                                                                                             caTCTTCTTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAAt                                 >  1:331340/1‑244 (MQ=255)
                                                                                                                                                                                                               tcttctTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAAt                                 >  1:877363/1‑242 (MQ=255)
                                                                                                                                                                                                               tcttctTGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATGGCTTGTGc                        <  1:696425/251‑1 (MQ=255)
                                                                                                                                                                                                                 ttcttGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGc                                                                                                                                         <  1:378452/136‑1 (MQ=255)
                                                                                                                                                                                                                 ttcttGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATGGCTTGTGCtt                      >  1:128578/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATGGCTTGTGCTTGCAc                  <  1:119070/251‑1 (MQ=255)
                                                                                                                                                                                                                     tGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATGGCTTGTGCTTGCAc                  <  1:360679/251‑1 (MQ=255)
                                                                                                                                                                                                                     tGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATGGCTTGTGCTTGCAc                  <  1:657748/251‑1 (MQ=255)
                                                                                                                                                                                                                     tGGGCAATCAAATCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATGGCTTGTGCTTGCAc                  <  1:368705/251‑1 (MQ=255)
                                                                                                                                                                                                                                aTCTTCCAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTc                                                                              <  1:884343/180‑1 (MQ=255)
                                                                                                                                                                                                                                     ccAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATGGCTTGTGCTTGCCCTTCGGTTAAATAatat  <  1:384573/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ccAGATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATGGCTTGTGCTTGCACTTCGGTTAAATAatat  <  1:346318/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gATTTATATCACCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACaataa                                                                                                                           <  1:858464/127‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
AGTTGATATACTTTCAGCCAGTATAAACGTCCGCGGCTGGAGAAGCAGAGAATGGTATCGTGGGTGTTGGCAACCAGTAATTTTTCGATAAAATCTTCTTCTTTCATTTTGGTTGCGGATTTACCTTTACCACCACGACGCTGTGCTTCGTAGTCGGTAATCGGTTGATATTTCACATAACCTTCGTGGGAAAGGGTCACCACCACATCTTCTTGGGCAATCAAATCTTCCAGATTTATATCGCCGGAAGCCGCCGTGATTTCGGTTCTGCGCGGATCGTTAAATTGCTCGCGGATTTGCTCCAATTCTTCACGAATCACTTCCATTAAACGCTCAGGGCTACGCAAGATGTACAATAATTTACCGATGGCATCCAATAATTTTTTGTATTCGTCCAGGATTTCATCGTGACCTAAATTGGTTAATTTTTGCAAACGTAAATCTAAAATGGCTTGTGCTTGCACTTCGGTTAAATAATAT  >  1/586965‑587444

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