Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,516,959 C→T P253P (CCG→CCA)  fecR ← ferric citrate regulator FecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,516,9590CT94.6% 246.0 / ‑5.8 95P253P (CCG→CCA) fecRferric citrate regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (1/1);  major base T (44/43);  minor base G (5/0);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 5.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑CGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGT  >  NC_000913/4516717‑4517204
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
gcGAACGAGAAGGGTAAAAAGGACAAGCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAGGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCGCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATc                                                                                                                                                                                                                                                     <  1:403821/250‑1 (MQ=255)
  gAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAA‑CGGGTTAGCGGGGCTATCACCTC‑ACCC‑AG‑TGGGTTATC‑Gc                                                                                                                                                                                                                                                  >  1:285162/1‑251 (MQ=255)
     cGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGa                                                                                                                                                                                                                                               <  1:664282/251‑1 (MQ=255)
      gagaAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGa                                                                                                                                                                                                                                               >  1:293106/1‑250 (MQ=255)
               aaaaaGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGAATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGTTCGCAGAGAAGCACACCGTTTCCGTGA‑CGGGCTGGCGAGGCTAACAACCC‑ACCC‑AA‑TGGGTTATT‑GCTGCAGCTCAggc                                                                                                                                                                                                                                      >  1:642241/1‑249 (MQ=255)
                 aaaGGACAATCGAAATAATAATTATTTTCCTTATAGTGGTGAAATGTTTATCCAGTACCGCGTAATAAACTGAATTTTAACGTGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGTTAACTTCCCTCTCAGGCTCAGCCCGGCAACGGCGGGATCGCATCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑TCTGAAGCTCAGGATGt                                                                                                                                                                                                                                   <  2:523342/251‑1 (MQ=255)
                   aGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAAGTTTAACGGGAAGCGTTTGCGCTATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCGCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTcc                                                                                                                                                                                                                                 <  1:248499/251‑1 (MQ=255)
                     gACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCtt                                                                                                                                                                                                                               <  1:500435/251‑1 (MQ=255)
                     gACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCtt                                                                                                                                                                                                                               <  1:205807/251‑1 (MQ=255)
                         aTCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTc                                                                                                                                                                                                                           <  1:174508/251‑1 (MQ=255)
                            gAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGATTGCGCGATAACGTTCAGGCTCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCCGCCGTTGCGGTAA‑CGGGGTAGCGTGGCTATCCCCTC‑CCCC‑AG‑TGGGTTATC‑GCCTCAGCTCTGGATGTCCTTCGTcccg                                                                                                                                                                                                                        >  1:86908/1‑249 (MQ=255)
                                    aaTTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAA‑CGGGTTAGCGTGGCTATCACCTC‑CCCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTAct                                                                                                                                                                                                                >  1:447496/1‑251 (MQ=255)
                                              cTTACAGTTGTTAAATGTTTACCCAGTACCGCGTAATAGACTGAATTTTAACGCGAAGCGTTGGCGCGATAACGTTCAGGAGCGCATCG‑GTATTTTTCAGCGGGAACGGCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCTGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑TCCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGt                                                                                                                                                                                                                            <  2:568806/229‑1 (MQ=255)
                                              cTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGAGGGAA‑CGGGTTAGCGTGGCTATCACCTC‑AACC‑AG‑TGGTTTATC‑GATGAAGCTCAGGATGTCCTTCGt                                                                                                                                                                                                                            >  1:568806/1‑229 (MQ=255)
                                                   aGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑CG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTCATCCAG‑Cgg                                                                                                                                                                                                >  1:668006/1‑251 (MQ=255)
                                                      ggtgAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑CCCC‑CG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTCATCCAG‑Cgggtt                                                                                                                                                                                             >  1:642709/1‑248 (MQ=255)
                                                                  aTCCAGTACCGCTGAATTGACTGAATGTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCAGCG‑GTATTTTGCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGGCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCaa‑a                                                                                                                                                                             <  2:351799/251‑1 (MQ=255)
                                                                         aCCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCaga                                                                                                                                                                       <  1:623328/250‑1 (MQ=255)
                                                                             cgTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGgcgc                                                                                                                                                                  <  1:572905/251‑1 (MQ=255)
                                                                               tAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGGTTGCCGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGGGTGGCTGTCCCCCC‑ACCC‑AG‑TGGGTTTTC‑GCTGAAGCTCAGGGTGTCCTTCGTCCAGCTTGTCCTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACCCAGAGgcgcgg                                                                                                                                                                >  2:385606/1‑249 (MQ=255)
                                                                                aaTAGGCTTAATTTTAACGGGAAGTGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTg                                                                                                                                                                <  1:124028/250‑1 (MQ=255)
                                                                                aaTAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑AGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCCGAGGCGCTg                                                                                                                                                                >  1:351246/1‑250 (MQ=255)
                                                                                aaTAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACCCCGGTGCGGTAA‑AGGGATAGCGGGGCCATCAACCC‑ACCC‑CG‑GGGTTTTTC‑GGCGAAGCTCAGGATGTCCCTTGTCCAGCTTGTACTCTCGTCATCCAC‑AGGTTTTCA‑T‑G‑CGCCAA‑AACCCGAGGCGCTg                                                                                                                                                                >  2:660739/1‑250 (MQ=255)
                                                                                   aGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCCGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGGTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACt                                                                                                                                                            >  1:367587/1‑251 (MQ=255)
                                                                                      cTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGAACGCAGCGCAGCACGCCGTTGCGGGAA‑CGGGGTAGCGTGGCTATCACCCC‑AACC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACA‑G‑CGCCCA‑ACTCAGAGGCCCTGAACttc                                                                                                                                                          >  2:284191/1‑248 (MQ=255)
                                                                                       tGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCCGCCGTTGCGGTAA‑CGGGGTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCCTCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGc                                                                                                                                                          >  1:284687/1‑249 (MQ=255)
                                                                                         aaTTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACATCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑CG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCCG‑CGGTTT‑CACT‑G‑CGCCCA‑ACTCAGAGGCGCTGCACTGCAGGc                                                                                                                                                      >  1:369704/1‑251 (MQ=255)
                                                                                             ttAACGGGAAGCGTTGGCTCGATAACTTTCAGGATCGCATCG‑GTATTTTGCAGCGGTAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATTTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCttt                                                                                                                                                   <  2:658978/250‑1 (MQ=255)
                                                                                              tAACGGTAAGCGTTGGCGCGATAACGTGCAGGATAGCGTCG‑GTATTTTTCAGCGGGAACGTCCCGCTCATGCGCAGCCCGGCACCGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGTCTATCACCTC‑ACCC‑AG‑TGGTTTAGC‑GCTGAAGCGCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTc                                                                                                                                                  <  2:588174/250‑1 (MQ=255)
                                                                                               aaCGGGAAGCGTTTGCGCCATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCGCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCa                                                                                                                                                 <  2:367587/250‑1 (MQ=255)
                                                                                                aCGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCCCGCCGTTGCGGTAA‑CGGGTTAGCGGAGCGATCCCCTA‑ACCC‑AGATGGTTT‑TC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGACATCCAG‑CGGGTT‑CACT‑G‑CGCCAA‑ACTCCGAGGCGCTGAAATGCAGGCTTTCAcc                                                                                                                                               >  1:25665/1‑251 (MQ=255)
                                                                                                 cGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCAcc                                                                                                                                               >  1:162396/1‑250 (MQ=255)
                                                                                                   ggAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATGGCAACGCAGCAAGCCGTTGCGGAAA‑CGGGTTAACGTGGCTATCAACTC‑ACCC‑AG‑TGGTTTATC‑GGTGACACTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CAGGTT‑CAAT‑G‑CGACAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTTACCAGc                                                                                                                                            >  1:107958/1‑251 (MQ=255)
                                                                                                   ggAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAACACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑CCCC‑CG‑GGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCCGAGGCGCTGAACTTCAGGCTTTCAcccgc                                                                                                                                            >  1:616647/1‑248 (MQ=255)
                                                                                                   ggAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCCCAGAGCCACACGCCGTTGCGGTAAACGGGT‑AGCGTGGCCATCACCTC‑ACCC‑AG‑GGGGTTATC‑GCTGAAGCTCCGGCTGACCTTCCTCCACCTTGTACTCTCGTCATCCCGTCGGTT‑‑CACT‑GGCCCCAA‑ACT‑AGAGGAGCTGAAcctcaggatttccccacc                                                                                                                                            >  2:662575/1‑234 (MQ=255)
                                                                                                     aaGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGt                                                                                                                                          >  1:387515/1‑251 (MQ=255)
                                                                                                      aGCGTTTGCGCGATAACGTTAAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGGTAGCGTGGCTATCACCCC‑ACCC‑AG‑TGGTTTTAC‑GCTGAAGCCAAGGAGGTCCTTCGCCCAGCTTGGACTCTCTCCATCCAG‑CGGTTT‑CACT‑G‑CCCCCA‑ACTCAGAGGAGCGAAACAGCAGGCTTTTACCAGCGtt                                                                                                                                         >  2:467375/1‑251 (MQ=255)
                                                                                                          gttGCGCGATATCGTTCAGGATCGCATCG‑GGATGTTGCAGCGGGTACGTTCCTCTCAGGCGCAGCCCGGCAACGGCGGGACCGCAGGTCAGCACGCCGTTGCGGGAA‑CGGGTTAGCGTGGCTAGCACCGC‑ACCC‑AG‑TGGGGTAGC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTGGTACTCGCTTCATCCAG‑CGGTTT‑CACT‑G‑CGCTAA‑ACTCAGAGGCGCTGAACGGCAGGCTGTCACCAGCGGTCACg                                                                                                                                     <  2:86908/250‑1 (MQ=255)
                                                                                                            ggcgcGATAACGTGTAGGATCCCA‑CGTGTGTTTTTCAGCGGAAACGTCCCGCGCAGGAGCAGCCCGGCAACGTCGGGATCGCAGCGCTGCACTCCGTTGCGGTAG‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCGGTAGCGCAGGATGTCCTTCGTCTAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCTGCGTTCACGAt                                                                                                                                   <  2:636325/250‑1 (MQ=255)
                                                                                                             gcgcGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATg                                                                                                                                  <  1:514150/251‑1 (MQ=255)
                                                                                                               gcgATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGTAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGGCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAAAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGCTGCg                                                                                                                                >  1:175727/1‑251 (MQ=255)
                                                                                                                cgATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCCGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCCGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGGACTGCAGGCTTTCACCAGCGTTCAAGATGCGt                                                                                                                               >  1:445614/1‑251 (MQ=255)
                                                                                                                cgATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGt                                                                                                                               <  1:170390/251‑1 (MQ=255)
                                                                                                                    aaCGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCCGCGCAGCACGCCGGTGCGGTAA‑CGGGATAGCGGGGCCATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGCTGTCCTTCGTCCAGCTTGTAATCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCCA‑ACTCAGAGGCGCTGAACTACAGGCTTTCACCAGCGTTCACGATGCGttttt                                                                                                                           >  1:74631/1‑251 (MQ=255)
                                                                                                                    aaCGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGttttt                                                                                                                           <  1:108412/251‑1 (MQ=255)
                                                                                                                    aaCGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGttttt                                                                                                                           <  1:284191/251‑1 (MQ=255)
                                                                                                                       gTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGGTAGCGGGGCCATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCCGGGTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCAACAGCGGTCACGATGCGGTTTTgcg                                                                                                                        >  1:132116/1‑251 (MQ=255)
                                                                                                                         tcatgatagCATCG‑GTATTTTTCAGCAGGAACGCCCCGCTCAGGCGCGGACCTGCAACTGCGGGAGCTCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGGACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCTCGATTCGTTTTTgcgc                                                                                                                       <  2:462605/242‑1 (MQ=255)
                                                                                                                             gATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGGTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCggggg                                                                                                                  >  1:560022/1‑251 (MQ=255)
                                                                                                                                cGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCa                                                                                                                <  1:511348/250‑1 (MQ=255)
                                                                                                                                cGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGAATGCAGAGCAGAAAGAGGCTGCGGTAA‑CGGGGTAGGGCGGCCATCACCTC‑GACC‑AG‑TGGGTTATCAGC‑GGCGCTAAGGGTGTCCTTCGCCAAGCGTGTATTCTCGCCATCCAG‑CGGGTT‑CACTCG‑CACCAACAC‑CAG‑CGCGCTACACCGCAGCCTTTCAACAGCGTTCACGATGCGgtttttagccggggcg                                                                                                                >  1:666805/1‑239 (MQ=255)
                                                                                                                                 gCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACt                                                                                                              <  1:162289/251‑1 (MQ=255)
                                                                                                                                 gCATCG‑GTATTTTTCAGCGGGAACGTCCCGATAAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCACCACCCCACTTCGGTAA‑AGGGGTAGCGAGGGTATCACCTC‑ACCCCAG‑‑GGGTTATC‑GCCAAAACAAAGGAAGTACTTCCTCCAACTTGTAATACCTTCATCCAC‑AGGTTT‑CCCA‑A‑CTCCAA‑ACTCAGAGGAACTTAACTGCAGGCCTTCACCCCCGTTCACGactaggttttgagacggcggccc                                                                                                              >  2:594667/1‑229 (MQ=255)
                                                                                                                                   agcG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCAc                                                                                                                                      <  2:400871/223‑1 (MQ=255)
                                                                                                                                   aTCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGGTGCGGGAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCCGGCTTTCACCAGCGTTCAc                                                                                                                                      >  1:400871/1‑225 (MQ=255)
                                                                                                                                     cG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGc                                                                                                           <  1:88817/250‑1 (MQ=255)
                                                                                                                                     cG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCg                                                                                                          <  1:362099/251‑1 (MQ=255)
                                                                                                                                         tATTTTTCAGCGGGAACGTACCGCTCAGGCGCAGCCCGGCAACGGCGGGCTCGCAGCGCAGCCCGCCGTTGCGGTAA‑CGGGCTCGCGGGGCCATCACCTC‑ACCC‑AG‑TGGGTTATC‑GCTCAAGGACAGGGTTGCCTTCCTCCCGCTTGTACTCTCGTCATCCAG‑CGGTGT‑CACC‑G‑GGCCAA‑ACTCCGCGGAGCTGAACAGCAGGCTTTCACCAGCGTTACCCATACGTTTTTGCGCGGGGAGACTGGCgcg                                                                                                        >  2:435297/1‑248 (MQ=255)
                                                                                                                                                   cGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑CGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGTGAGAAGCACTTccc                                                                                             >  1:318910/1‑250 (MQ=255)
                                                                                                                                                   cGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGGGGCTATCACCTC‑CCCC‑AG‑TGGGTTATA‑GCCGAAGCTCAGGCTGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGGTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCTCCGAACTGCAGGCTTTCACCAGCGTTCACCATACGTTTTTGTGCGGGGCGACTGGCGCGAAggccttcac                                                                                             >  2:439379/1‑243 (MQ=255)
                                                                                                                                                     ggAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGGTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CCCT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCCGGCTTTCACCAGCGTTCACGATGCGTTTTTTCGCGGGGGCACTGGCGAGAAGCACTTCcacc                                                                                           >  1:435973/1‑250 (MQ=255)
                                                                                                                                                          gTCCCGCTCAGGCGCAGCCCGGCACCGGCGGGATCGCAGCGCAGCACGTCGTTGCGTTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACGCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtg                                                                                      <  2:598189/251‑1 (MQ=255)
                                                                                                                                                            cccGCTCATGCGCAGCCGTGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTTA‑CGGGTTAGCGTTGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTTCTCTCGTCATCCAG‑CGGTTT‑CTCT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGGTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgct                                                                                    <  2:447496/251‑1 (MQ=255)
                                                                                                                                                             ccGCCCAGGCGCAGCCCGGCAACGGCGTGATCGCTGCGCGTCACTCCGTTGCGGTAA‑CGTGTTAGCGTGGCTATCACCTC‑TCCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGAGGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CCCT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCTCTTCCCCAGCGtgctg                                                                                   <  2:449960/251‑1 (MQ=255)
                                                                                                                                                              cGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGGAA‑CGGGGTAGCGTGGCTATCACCTC‑ACCC‑AG‑GGGTTTATC‑GCTGAAACTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCCA‑ACTCAGAGGCGCTGAACTGCCGGCTTTCACCCGCGTTCCCGATGCGTTTTTGCGCGGGGGCCCTGGCGAGAAACACTTCCACAGCGggctgt                                                                                  >  2:172785/1‑250 (MQ=255)
                                                                                                                                                                   aGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgctgctgaac                                                                             <  1:179737/251‑4 (MQ=255)
                                                                                                                                                                      cgcAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTCGCGTGGCTATCACCTC‑ACAC‑AG‑TGGTTTTTC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CCCT‑G‑CGCCAA‑ACTCAGAGGCCCTGAACTGCAGGCTTTCACCAGCGTTCCCGATGCGTTTTTGCGCGGGGGCACTGGCGCGAAGGACTTCCACAGCGTGCTGCTGAACGTc                                                                          >  1:634668/1‑251 (MQ=255)
                                                                                                                                                                          gCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCCAGc                                                                      >  1:22214/1‑251 (MQ=255)
                                                                                                                                                                                cAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTg                                                                <  1:223241/251‑1 (MQ=255)
                                                                                                                                                                                cAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTg                                                                <  1:467011/251‑1 (MQ=255)
                                                                                                                                                                                    ggcggGATCGCAGTGCAGCACGCCGTTTCTTTAT‑TGGGCTAGCGTGCCTTTCCCCTC‑CCCC‑AG‑TCGTTTATC‑GCTGAAGGTCAGGATGTCCTTCGTGCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCTA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTGCACGATGCGTTTTTTCGCGGGGGCACTGGCGAGAAGCACTGCCACTGCTTGCTGCTGAACGTCAAGCTGCGTGAAAt                                                            <  2:410395/251‑1 (MQ=255)
                                                                                                                                                                                     gcggGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCTCAGCGTGCTGCTGAACGTCAAGCTGCGTGAAAtt                                                           <  1:176172/251‑1 (MQ=255)
                                                                                                                                                                                     gcggGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAAtt                                                           >  1:589717/1‑251 (MQ=255)
                                                                                                                                                                                        tgctcgccccgcagcACGCCGTTGCGGCCA‑CGCGATTGCGAAGCTATCCCCCC‑ACCC‑CG‑TGCTTTCTC‑GCTCCAGATCATGATGTCCTTCGACAATCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACC‑G‑CGCACA‑ACTCCGACGCCCTGAACCGCAGGCTTTCAACCTAGTTCACGCTACGTTTTTACGCGGGGGCCCTGGCAAGAAGCACTTCCACAGCGTGATGCTGAAAGTCAAGCTGCggtctattacc                                                        >  2:341489/10‑241 (MQ=255)
                                                                                                                                                                                            cgcagcgcagCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGt                                                                                       >  1:517100/1‑216 (MQ=255)
                                                                                                                                                                                            cgcagcgcagCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCctg                                                     <  1:45141/250‑1 (MQ=255)
                                                                                                                                                                                             gcagcgcagcACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCctgct                                                   >  1:48254/1‑251 (MQ=255)
                                                                                                                                                                                               agcgcagcACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTgg                                                 >  1:679571/1‑251 (MQ=255)
                                                                                                                                                                                                gcgcagcACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGc                                                <  1:42633/251‑1 (MQ=255)
                                                                                                                                                                                                    agcaCGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGAc                                            <  1:583683/251‑1 (MQ=255)
                                                                                                                                                                                                       ccGCCGTTGCGGTAA‑CGGGTTAGCGTGACTATCACCTCTACCC‑‑G‑GGGTTTATA‑GCTGAAGCTCAGGATGACCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCGGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAGGCTCTTCCCAAGCGTGCGGCTGATCGTCCAACTGCGTAAAATTCTCCTTCGGGGGGACggg                                         >  2:368197/2‑250 (MQ=255)
                                                                                                                                                                                                                             aGCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGCCGGTAAATTCTGTCCCTAAAGCAgg                    >  1:131453/1‑250 (MQ=255)
                                                                                                                                                                                                                              gCGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTg                   <  1:385606/251‑1 (MQ=255)
                                                                                                                                                                                                                               ggtggcttgcacctc‑tccC‑GG‑TGGTTGATC‑GCTGAAGCTCAGGAGTTCCTTCGTCCAGCTTGTACTCTCGGCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGa                  <  2:15954/242‑1 (MQ=255)
                                                                                                                                                                                                                               cGTGGCTCTCACCTC‑ACCC‑AG‑TGGTTTCTC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCGGGCGGACGGGAAATTCTGCCCCTAAAGCAGTgg                  >  2:350731/1‑250 (MQ=255)
                                                                                                                                                                                                                               cGTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCGCAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGa                  <  1:18559/251‑1 (MQ=255)
                                                                                                                                                                                                                                gTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAACGCAGTg                   >  1:448741/1‑249 (MQ=255)
                                                                                                                                                                                                                                gTGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTg                   >  2:322847/1‑249 (MQ=255)
                                                                                                                                                                                                                                 tGGCTATCACCTC‑ACCC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGc                >  1:518223/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggCTATCTCCTC‑GCCC‑AG‑CGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGGCATCCTG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTTCACGTCAAGCTGCGTGAAATTATCCTGCTGGTGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCt               <  2:333333/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aTCACCTC‑ACCC‑AG‑TTGTTTTAC‑GCTGGAACTCCGGATGTCCTTTGTCCTGCTTGTACTCTCGTCAACCAG‑CCGTTT‑CACT‑G‑CGCCAA‑ACTCAAAGGTGCTGGACTGCCGGCTTTCAACAGCCTTCCCGATGCGTTTTTGCTCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGcc           <  1:100337/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tcacctc‑accC‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑AACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGccc          >  1:21824/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ccc‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGtgct                                                                                    <  1:372135/169‑1 (MQ=255)
                                                                                                                                                                                                                                                ccc‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGAACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTACATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGt  >  2:107852/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ccc‑AG‑TGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGAGtgct                                                                                    >  2:372135/1‑169 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           
TCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTCCTTACAGTGGTGAAATGTTTATCCAGTACCGCGTAATAGACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCG‑GTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAA‑CGGGTTAGCGTGGCTATCACCTC‑ACCC‑AG‑CGGTTTATC‑GCTGAAGCTCAGGATGTCCTTCGTCCAGCTTGTACTCTCGTCATCCAG‑CGGTTT‑CACT‑G‑CGCCAA‑ACTCAGAGGCGCTGAACTGCAGGCTTTCACCAGCGTTCACGATGCGTTTTTGCGCGGGGGCACTGGCGAGAAGCACTTCCACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGT  >  NC_000913/4516717‑4517204

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