Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 840,656 A→G 14.1% T292T (ACT→ACC)  fiu ← catecholate siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913840,6560AG14.1% 178.3 / 3.4 84T292T (ACT→ACC) fiucatecholate siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (40/29);  new base G (10/3);  total (52/32)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACATCTTTGGTATTCGCCGTGCGTGACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGACTTTGTCGCGTCCGGCATCGTGCGTTTTTTCGCCCATTACATTC  >  NC_000913/840421‑840893
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             
tCACATCTTTGGTATTCGCCGTGCGTGACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTa                                                                                                                                                                                                                                            >  1:323748/1‑239 (MQ=255)
    aTCTTTGGTATTCGCCGTGCGTGACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTTCTCAAACCGCATGGTGGCGGGGTCGGGGGTCGAATCGTCGGAATCGGAAACCGTGCCGGTAAATTTAT‑TAggcaca                                                                                                                                                                                                                                         >  1:1325582/1‑233 (MQ=255)
    aTCTTTGGTATTCGCCGTGCGTGACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGGGCCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAATAGTTAT‑GAGTAtg                                                                                                                                                                                                                                          >  1:717335/1‑236 (MQ=255)
          ggTATTCGCCGTGCGTGACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTc                                                                                                                                                                                                                                  <  1:850786/239‑1 (MQ=255)
              ttCGCCGTGCGTGACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTTT‑GAGGATAACCTTTTCCGGa                                                                                                                                                                                                                              >  1:703110/1‑239 (MQ=255)
              ttCGCCGTGCGTGACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGGGGTGTTATCGTTGGTTTCGTGCTCAAACCGCATGGTGGCGGTGTCGGGGGTCGGATGGTCGGACACGGAATCTGTGCCGGAAAAGCTAT‑GTGGAGCAAATTTTTCGGa                                                                                                                                                                                                                              >  2:1336027/1‑239 (MQ=255)
               tCGCCGGGGGGGCCCAGGTCCAGCTTTTCCCATCGCTGGGGGGCTGAGTAATTGTCGACGCCCCGCCCATAATCGG‑CGTCATCAGGTAATCCTGCTTTACGGGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGaa                                                                                                                                                                                                                             <  2:703110/239‑1 (MQ=255)
                cGCCGTGCGTGACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCCAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAATGTAT‑GAGTATCACCTTTTCCGGAAt                                                                                                                                                                                                                            >  1:1311193/1‑239 (MQ=255)
                cGCCGTGCGTGACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGCCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GGGTATCAACTTTTCCGGAAt                                                                                                                                                                                                                            >  1:1145144/1‑239 (MQ=255)
                     tgcgtgACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGGGTCGGTGGGCGAATCGTCGTAATCGGAATCCGGTCCGGAAAAGTTATGGGGAATCAATTTTTcggataaatat                                                                                                                                                                                                                        >  2:705915/1‑229 (MQ=255)
                      gcgtgACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGGGTTATCGTTGATGTCGGGCTCAAAACGCATGGTGGCGGGGTCGGGGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAATTTAT‑GGGTATAAACTTTTCCGGAAGAATTCa                                                                                                                                                                                                                      >  2:467933/1‑239 (MQ=255)
                                     gCTTTTCACAGGGTTGGGGGGTTGGGTAAGTTCCGACCCCCCGCCATAAATCGC‑CGGCATCAGGTAATCCTGCTTTACGGGCGGCGAACGGGTAGTATTGCGAATGGTGGTGTTAGCGTTGATGTCGTGCCCAAAACGCATGGTGGCGGTGTCGGGGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAG‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCg                                                                                                                                                                                                       <  2:66245/239‑1 (MQ=255)
                                         ttCCCATGGCTGGGTGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAAGGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGa                                                                                                                                                                                                   <  1:705915/239‑1 (MQ=255)
                                          tCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGGCGAATCTTCGTAATCGGAATCCGGGCCGGAAAAGTTAT‑GGGTTTAAACTTTTCCGGAATGATTCAGGGGCGCCGTTTCCCCCGAt                                                                                                                                                                                                  >  2:1019839/1‑239 (MQ=255)
                                           cacaTCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGGGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGCATGATTCAGGGCCGCCGTTCCCGCAGATg                                                                                                                                                                                                 >  1:1221967/1‑239 (MQ=255)
                                              ctccttggtggTCTGTGTTATTCTCGCCGCCCCCCCCAAAACCGC‑AGTCATCAGGTAATCCTGCTGTCCGCGCGACCAAGGGGTATTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAACCGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGGGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATgggg                                                                                                                                                                                              <  2:1221967/234‑1 (MQ=255)
                                              aTCGGTGGGGGGCTGAGTAATATTCGACGCCCCGGCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATgggg                                                                                                                                                                                              <  1:1165902/239‑1 (MQ=255)
                                                cGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGGGTCGGTGGGCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCACATTTTCCGGGATGATTCAGGGCCCCCGGTCCCGCAGATGGGGCa                                                                                                                                                                                            >  2:1030105/1‑239 (MQ=255)
                                                   tggtggGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCCCCGTTCCCGCAGGTGGGGCCGaa                                                                                                                                                                                         >  1:491053/1‑239 (MQ=255)
                                                            gAGTTATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCAGCCAACGGGTAGTATTGCGAATGGTTGTGTTATCGTTGATGTCGTGCTAAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAAGCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCgg                                                                                                                                                                                 <  1:113588/238‑1 (MQ=255)
                                                             aGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GGGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCa                                                                                                                                                                               >  1:223808/1‑239 (MQ=255)
                                                                     ttgcacccccccccATAATCGCACTCCAGTAGGTAATCGTGCTTTCGGGGGGAACAACGGTTAGTATTGCGAAGGGGGGGTTTATCGTTGAGGTCGTGCTCAAAACGCATGGGGGCGGTGGCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGAt                                                                                                                                                                        <  2:929286/233‑1 (MQ=255)
                                                                       gccgccccgcccATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGGGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGt                                                                                                                                                                     <  1:283166/237‑1 (MQ=255)
                                                                       gccgccccgcccATAATCGC‑CGGCATCAGGTAATCCTGATTTACGAGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGCGGCGGTGTCGGCGGTCGTATCGTCGTAATCGGACACCATGCCGCAAAAGTGCT‑TGGTATCAACTTTTCCGGAACGATTAAGAGCCGCCGTCCCCGCAACTGGGCAAGAATAAACCGGCaacacgcacgc                                                                                                                                                                     >  2:296301/3‑230 (MQ=255)
                                                                       gACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTGCGTGGTGGAATTGTCGTAAGCGGAATCGTGGGCGTCAAAGTTAT‑GAGCATCAACTTTTCCGGAAGGAGTACGGGGCGCCGGTCCCGCAACTGGTGAGGAAAAGCCCGGCAAACCGATGGt                                                                                                                                                                     >  2:280691/1‑239 (MQ=255)
                                                                       gACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGTGGGTCCGCTGCTCGTAATCGGAATCACGGCCGTAAAAGTTAT‑GGGTAACAAAGTTTCCGGACTGATTCAGGGCCCACGGTCCCCCAGGTGGCGCCGAAGAGcagcacaaaacgacgag                                                                                                                                                                     >  2:897069/1‑223 (MQ=255)
                                                                       gACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GGGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCGGATGGGGCAGAATAACCCGGCAAACCGGTggc                                                                                                                                                                     >  1:381054/1‑238 (MQ=255)
                                                                       gACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAATTAT‑GGGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCCCAGATGGGGCAGAATAGCCCGGCAAACCGATgg                                                                                                                                                                      >  1:901203/1‑238 (MQ=255)
                                                                       gACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGGGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATAATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAAAAACCCGGCAAAcgagaggg                                                                                                                                                                     >  2:674663/1‑232 (MQ=255)
                                                                       gACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGCGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGGCGGTGGTCGAATCGTCGTAATTGGTATCCGTGCCGTAAAAGTTAT‑GCGTATCAACTTTTCCGGAATTATTCGAGGCCTCCGTTCCCGCAGATGGGGGCGAAAAAACAGGAAAACAGATgga                                                                                                                                                                     >  2:205270/1‑238 (MQ=255)
                                                                         cgccccgcccATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCCGATGGTGCAGAATAGCCCGGCAAACCGATGGTCg                                                                                                                                                                   >  1:771281/1‑239 (MQ=255)
                                                                                cccATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGGATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCCCAGATGGGGCAGAATAGCCCGGCAAACCGATGGCCGGAATgc                                                                                                                                                             >  1:409842/1‑238 (MQ=255)
                                                                                  gaggggtataagagcaagggTAATCCGGCTTTACGCGCTACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATgcc                                                                                                                                                            <  1:850879/220‑1 (MQ=255)
                                                                                        cGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTg                                                                                                                                                    <  1:1075613/239‑1 (MQ=255)
                                                                                               catcatGTAATCCCGGTTGACGCTGAACCAACGGGTAGTATTGGGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCgtg                                                                                                                                               <  1:1283069/238‑1 (MQ=255)
                                                                                               catcaGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCgtg                                                                                                                                               <  1:1019839/238‑1 (MQ=255)
                                                                                                   aGGTAATCCTTCTTTACGGGCGACCAACGGGTAGTATTGGGAATGGTGGTGTTATCGTTGATGTCGTGCTAAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTtgtg                                                                                                                                          <  1:467933/239‑1 (MQ=255)
                                                                                                      tAATCCTGCTTTACGCGCGCCCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTtgtg                                                                                                                                          <  1:1024082/236‑1 (MQ=255)
                                                                                                          ccTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGAAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTggggg                                                                                                                                   >  1:609214/1‑237 (MQ=255)
                                                                                                               tttACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGTTGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGCATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GGGGATCAACTTTTCCGGCAGGATTCAGGCACGCCGTTCCCGCAAATGGGGCAGAAGAGCCGGGCAAACCGATGGACGGGAGCGCCCCGTCTggacctgtccttttaaccaa                                                                                                                                 >  2:850879/1‑218 (MQ=255)
                                                                                                               tttACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAAAAACCCGGCAACCCGATGGTCGGAATGCCGCCCTCTGGCGGTTTGGGCTGGGGGACAt                                                                                                                               >  1:1334516/1‑238 (MQ=255)
                                                                                                                         aCCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCCGATGGGGCCGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGATGGGCTGGGTGACATGCAGATAATTa                                                                                                                    >  1:789290/1‑239 (MQ=255)
                                                                                                                         aCCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGGGGCGGGGTCGGTGGTCGAATCGGCGTAATCGGAATCCGGTCCGTAAAAGTTAT‑GAGGATCAACTTTTCCGGAATGCTTCAGGGCCGCCGTTCCCCCAGATGGGGCAGAAAAGCCCGGCAAACCGCTAGTCGGAAGCCCGCCGGCTGGGGGGTTGTTCGGGGGGAccgcgcgacaatta                                                                                                                    >  2:1256574/1‑226 (MQ=255)
                                                                                                                                     taTTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGGGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCAGCAAAACGATGGACGGAATGCCGCCGGCTGGGGTGTTGGGTGGGGGGACATGGAGAAAAATAAAAAACAAAAGa                                                                                                        >  2:1211285/1‑239 (MQ=255)
                                                                                                                                      tttGCGGATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGa                                                                                                        <  1:1256574/237‑1 (MQ=255)
                                                                                                                                         gggAATGGTGGTGTTATCGTTTATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGCAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCg                                                                                                    <  1:1030105/237‑1 (MQ=255)
                                                                                                                                          cGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGGGACATGCAGATAATTAAGATACAAACGATTc                                                                                                     >  1:290065/1‑237 (MQ=255)
                                                                                                                                           gAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGACCGGCGTGTTGTTCTGGGTGACATGCAGATAATTAAGATCCAAACGATTc                                                                                                     >  1:944844/1‑236 (MQ=255)
                                                                                                                                               ggtggtGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCACCTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCCGATGGGGCAGAAAAGCCCGCACAACCGCTGGTCGGAATGCCGCCGTCTGCCGTGTTGTGCTGGGGGACATGCAGTTAATAAAAAAAAAAACAATTGGCTGAAc                                                                                              >  2:15476/1‑239 (MQ=255)
                                                                                                                                               ggtggtGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGGAAATTAAGATACAAACAATTCGCTGTAc                                                                                              >  1:839447/1‑239 (MQ=255)
                                                                                                                                               ggtggtGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTACCAACTTTTCCGGAATGATTCAAGGCCCCCGTTCCCGCAGATGGGGCGGAATAGCCCGGCAAACCGAGGGATGGAATGACGCCGGCTGGTGTTGTTCGCTGGGGGACACGCAGCTAAtataaaaaaaaaagattgctctca                                                                                               >  2:407350/1‑215 (MQ=255)
                                                                                                                                                 tggtgTTATCGTTGATGTCGTGCTCAAAACGCGTGGTGGCGGTGTCGGTGGTCAAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTAcc                                                                                             <  1:205270/238‑1 (MQ=255)
                                                                                                                                                 tggtgTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGGGACATGCCGATAATTAAGATACAAACg                                                                                                         >  1:1148487/1‑226 (MQ=255)
                                                                                                                                                 tggtgTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGGTAT‑GAGTATCCACTTTTCCGGAATGATTCAGGGCCGCAGTTCCCGCAGATGGGGCAGAAAAGCCCGGCAAACCGATGGTCGGAATGCCGCCGCGGTGCGGGGTTGGTCGGGTGACATGAAGCAAATTAAAATACAAAcagttcacgggaaca                                                                                            >  2:326061/1‑225 (MQ=255)
                                                                                                                                                       tttgttgggttggttGCTCAAAAGGAGTGGGGGGGGTGGCGGTGGTGGAATGGTCGTAATGGGAATCCGTGCCGGAAAAGTTAT‑GAGTATGAACTTTTCCGGAATGATCCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCa                                                                                      <  2:1227116/225‑1 (MQ=255)
                                                                                                                                                              gATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGGCGGAATGCCGCCGTCTGGGGGGTTGTGCTGGGTGACATGCAAGTAATTAAAAAACAAACGATTCCCTGTACCAAGGCCAAAAGGGa                                                                               >  2:1130368/1‑239 (MQ=255)
                                                                                                                                                               aTGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGa                                                                                                                                                                                                   <  2:1148345/122‑1 (MQ=255)
                                                                                                                                                               aTGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGa                                                                                                                                                                                                   >  1:1148345/1‑122 (MQ=255)
                                                                                                                                                                    gTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAg                                                                         >  1:17411/1‑239 (MQ=255)
                                                                                                                                                                                     ttttaaagggaagggtCGAATCGTCGTAAGCGGAATCCGGGCCGTAAAAGTTAT‑GGGTATCAACTTTTCCGGAATGATGCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGc                                                        <  2:992236/226‑1 (MQ=255)
                                                                                                                                                                                      gggtggctgggtggtgacagcgtggtAAGCGGAGTGCGTGCCCGAAAAGGTGT‑GAGTATCACAGTTCCCGGATGTATTAGGGGCCGCCGTGCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGAAATGACGCCGTCTGGCGTGTTGTGCTTGGTGACATGCAGATAAGTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGAGGTAAAGc                                                        <  2:410063/218‑1 (MQ=255)
                                                                                                                                                                                      ggggttttCGGTGGTTCATTCGTCGTAATCGGGACCCGGGCCGTAAAATTTAG‑GAGTATCAATTTTTCCGGAATGATTCGGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGc                                                        <  2:223808/231‑1 (MQ=255)
                                                                                                                                                                                      gggggTGGCGGTGGTCAAATTGTCGTAAGGGAAATCCGGGCCGTAAAAGTTAT‑GAGTATCAACTTTGCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGc                                                        <  2:17411/235‑1 (MQ=255)
                                                                                                                                                                                      ggCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGAAAACCGATGGTCGGGATTCCGCCGACGCGCGGGTTGGTCTGGGGGACTGGCCGATAATTAGGATACAAAACATTGGCGGTACAAAGCCCAAAAGCGACAAAAAGGGCGAcccccaaaagc                                                        >  2:70792/1‑231 (MQ=255)
                                                                                                                                                                                            gtCGGTGGTCGAATCGGCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAATACCAACTTTTCCGGAATGGTTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCAGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCAtttt                                                 <  1:1211285/239‑1 (MQ=255)
                                                                                                                                                                                             tCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGaa                                                                                                                                                                                                                             >  1:471561/1‑66 (MQ=255)
                                                                                                                                                                                             tCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGaa                                                                                                                                                                                                                             <  2:471561/66‑1 (MQ=255)
                                                                                                                                                                                               ggtggtCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCGGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCGGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCCCTTTACCAAGGCCCAAAGGGCAAAAAGGGGGGACGCCGCAACACCCATTTTTg                                               >  2:537818/1‑238 (MQ=255)
                                                                                                                                                                                               ggtggtCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGGGCCATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGGGACGCCGTAAACCCCATTTTTGa                                              >  2:411957/1‑239 (MQ=255)
                                                                                                                                                                                                 ggtggcaaccgtcgaAATCGGAATCCGTGCCGGAAAATTTAT‑GAGTTTCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGAc                                             <  1:280691/229‑1 (MQ=255)
                                                                                                                                                                                                      ctttcgtgttatgcgtaaTCCGGTCCGTAAAATTTAT‑GGGGATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGGTGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGACTTTGTc                                       <  2:1269669/223‑1 (MQ=255)
                                                                                                                                                                                                              gtAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGGGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGCCCGAAGGGGCGACGCCGTAACGCTCATTTTTGa                                              >  1:992236/1‑224 (MQ=255)
                                                                                                                                                                                                                      aaTCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGCAGCCGCAACGCTCATTTTTGACTTTTTGCCGTCCGGCATCGTGc                       >  1:972931/1‑239 (MQ=255)
                                                                                                                                                                                                                       aTCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGACTTTGTCGCGTCCGGCATCGTGCg                      <  1:1130368/239‑1 (MQ=255)
                                                                                                                                                                                                                            tGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCCGGTAATTAAGATACAAACGATTCGCTGGACCAAGGCCAAAAGCGAAACAAGGGGCGACGCCGTAAAGCCCATTTTTGAATTTGACACGTCCGGAATCGGGGGgtttt                 >  2:291725/1‑239 (MQ=255)
                                                                                                                                                                                                                            tGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGCTACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGACTTTGTCGCGTCCGGCATCGTGCGttttt                 >  1:1001352/1‑239 (MQ=255)
                                                                                                                                                                                                                            tGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAACCGATTCGCGGTACCAAGGCCAAAAGCGACGGAAGGGGCACCGCCGAAACGCTCATTTTTGACTTTTTCGCGTCCGGCACCGTGCGttttt                 >  1:351072/1‑239 (MQ=255)
                                                                                                                                                                                                                                gTAAAAGGTCT‑GAGTATCAACTTGTCCGGAATGAGTCGGGGCCGCCGTTCCCGCAAATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGACTTTGTCGCGTCCGGCATCGTGCGTTTTTTCGc             <  1:291725/239‑1 (MQ=255)
                                                                                                                                                                                                                                ctAAAGGTTTT‑GGGTTTCAATTTTTCCGGAATGATCCGGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGACTTTGTCGCGTCCGGCATCGTGCGTTTTTTCGc             <  2:290065/238‑1 (MQ=255)
                                                                                                                                                                                                                                       ttAT‑GAGTATCAACTTTTCCGGAAGGATTAGGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGAtt                                                                                                      <  2:1316407/143‑1 (MQ=255)
                                                                                                                                                                                                                                         aT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCCTTTTTGACTTTGTCGCGTCCGGCATCGTGCGTTTTTTAGCCCATTAAAt    >  1:648567/1‑239 (MQ=255)
                                                                                                                                                                                                                                         aT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGACTTTGTCGCGTCCGGAATCGTGCGTTTTTTCGCCCATTACa     >  1:660175/1‑238 (MQ=255)
                                                                                                                                                                                                                                         aT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGAGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCCAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGACTTTGTCCCGCCCGGCCTCGTGCGTTTTTTCGCCCATTACa     >  2:351343/1‑238 (MQ=255)
                                                                                                                                                                                                                                            gAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGACTTTTCCGCGTCCGGCATCGTGCGTTTTTTCGCCCATTACCTTc  >  1:838849/1‑239 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                             
TCACATCTTTGGTATTCGCCGTGCGTGACCAGGTCCAGCTATTCACATCGCTGGTGGGCTGAGTAATATTCGACGCCCCGCCCATAATCGC‑CGTCATCAGGTAATCCTGCTTTACGCGCGACCAACGGGTAGTATTGCGAATGGTGGTGTTATCGTTGATGTCGTGCTCAAAACGCATGGTGGCGGTGTCGGTGGTCGAATCGTCGTAATCGGAATCCGTGCCGTAAAAGTTAT‑GAGTATCAACTTTTCCGGAATGATTCAGGGCCGCCGTTCCCGCAGATGGGGCAGAATAGCCCGGCAAACCGATGGTCGGAATGCCGCCGTCTGGCGTGTTGTGCTGGGTGACATGCAGATAATTAAGATACAAACGATTCGCTGTACCAAGGCCAAAAGCGACAGAAGGGGCGACGCCGTAACGCTCATTTTTGACTTTGTCGCGTCCGGCATCGTGCGTTTTTTCGCCCATTACATTC  >  NC_000913/840421‑840893

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