Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,512,642 T→C V256V (GTA→GTG)  fecC ← ferric citrate ABC transporter membrane subunit FecC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,6420TC96.0% 220.2 / ‑3.2 76V256V (GTA→GTG) fecCferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (34/39);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATG  >  NC_000913/4512406‑4512881
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
aaTTTTCATCCTCGCCTCCTCACAAGCCAGACAAGGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGCT‑CT‑CCGGGGAAGTCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGCCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCGGCCCAGAAGCGCGCCCGATGTGGCACCAGCAGACCGCTCAAC‑‑CCACCCGGACCTGCCACGCtaacccacac                                                                                                                                                                                                                                   >  2:78870/1‑242 (MQ=255)
   tttCATCCGCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACgcacgcac                                                                                                                                                                                                                                 <  1:83857/250‑1 (MQ=255)
        cccccccctcccCACACGCCAGAAAAAGCAAGGGC‑GGCAATACAGCGCCAGCCCTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCaccaacca                                                                                                                                                                                                                           <  2:406016/245‑1 (MQ=255)
         ccTCGCCTCCTCCCAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCGGCGGGCAGAT‑CT‑CCGGGAAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCAGCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAg                                                                                                                                                                                                                          <  2:28929/251‑1 (MQ=255)
            cGCCTCCGCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCAGGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAACC‑GCCACC‑GGGCCTGCCACGCTGACGCACGCACCAACCAGAag                                                                                                                                                                                                                       <  1:20986/251‑1 (MQ=255)
              cctcctCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAagc                                                                                                                                                                                                                      >  1:127626/1‑250 (MQ=255)
                 cctcACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCCCGGC‑GCCGGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGCCCTGCCACGCTGACGCACGCACCAACCAGAAgcagca                                                                                                                                                                                                                  <  2:483227/251‑1 (MQ=255)
                 cccaacaAGCAAGAAAAAGCAAGGGCTGCCAAT‑CAGGGCCAGCACTGC‑GCCTGCGGGCAGAC‑CT‑CCGGGGAAGGCCAGTGCGCGGGCGAGTACATCTGCCAGCAGCATAAGCGTGGCCCCCAGCAGCATGCTCACTGGCCGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGAAGACCGATAAAC‑GCCACG‑GGACCAGCCACGCTGACGCACGCACCAACCAGAAgcagca                                                                                                                                                                                                                  <  2:224703/247‑1 (MQ=255)
                                aaaGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGACGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGCTAAAC‑GCCCCC‑GGACCTGCCCCGCTGACGCACGCACCACCCAGAAGCCACACTAACCTATtgatc                                                                                                                                                                                                    >  2:400198/1‑249 (MQ=255)
                                aaaGCAAGGGCTGCCAAC‑AAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGGGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgat                                                                                                                                                                                                     <  1:400198/249‑1 (MQ=255)
                                 aaGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGAc                                                                                                                                                                                                  <  1:288981/251‑1 (MQ=255)
                                      aGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTTATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAAc                                                                                                                                                                                             <  1:167653/251‑1 (MQ=255)
                                        gccggCCAAT‑CACCGCCAGCAC‑GCAGCCTGCGGGC‑GATGCC‑CCGGGGAAGGCCAGTGCGCGTGCGAGCACATCTGCCACCAGCAGCGGCTTGGCCCCCAGCAGAATGCCCACTGGCAGTAGTTTGCGCTGTCCGAAGCCTGCCCAGAAGCGCGCCAGAGGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGAACTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGt                                                                                                                                                                                           <  2:305247/247‑1 (MQ=255)
                                        gccggCCAAC‑CACCGCCAGCACTGC‑GCCGGCGGGCAGA‑‑CTGCCGGGGAAGGCCAGTGCGCGGGCGAGTCCAGCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGt                                                                                                                                                                                           <  2:417939/247‑1 (MQ=255)
                                             ccAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGccc                                                                                                                                                                                      >  2:417566/1‑250 (MQ=255)
                                               aaT‑CAGCGCCAGCACTGG‑CCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGGGTACATCTGCCAGCAGCAGCAGCGTGGCCCCCAGCAGCAGGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCTATAAAC‑GCCTCC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCg                                                                                                                                                                                    <  2:480266/251‑1 (MQ=255)
                                                aT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGCCCTGCCACGCTGACGCCCGCCACAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGt                                                                                                                                                                                   >  2:152084/1‑251 (MQ=255)
                                                   cagcgccagcgcTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGACCAGTGCGCGTGCGAGTACATCTGCCAACAGCATCAGATTGGCCCCCAGCAGCATGCTCACTGGCAGTCCGTTGCGCTGATCGAAGCCTGCCAAGAAGCGCGCAAGATGTGGCACCAGCCGACCGATAAAC‑GCCACC‑GGACCTCCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACACATTGCTGACCAGACGCACCCTCgcc                                                                                                                                                                                  >  2:352871/1‑248 (MQ=255)
                                                         catcACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCACAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GTACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCa                                                                                                                                                                           <  1:170926/248‑1 (MQ=255)
                                                            cACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTc                                                                                                                                                                        >  2:91177/1‑251 (MQ=255)
                                                             aCTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACCTCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCCCGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTcc                                                                                                                                                                       >  2:60617/1‑251 (MQ=255)
                                                                gc‑gcCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTccc                                                                                                                                                                      <  1:349777/249‑1 (MQ=255)
                                                                       gCGGGCAGAC‑CT‑CCGGGGAAGGCCAGTGCGCGTGCGAGGACCTCTGCCCGAAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCTGTCCGTTGCGCTGATCGAAGCCTGCCCAGAGGCGCGCCAGATGTGGCTCCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCCCTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTa                                                                                                                                                                <  2:275679/249‑1 (MQ=255)
                                                                        ggtgccagtcct‑cCGGGGCAGGCCCGTGCGCGTGCGAGTACATCTGCCAACAGCCTCAGCGTGCCCCCCAGCAGCTTGCTCACTGGCAGTCCGTTGCGCTGACCGAAGCCTGCCCAGAAGCGCGCCAGAGGTGGCACCAGCAGACCGATAACC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATg                                                                                                                                                              <  2:386581/241‑1 (MQ=255)
                                                                                     ccGGGGAGGGCCCGTGCGCGTGCGAGTACATCCGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCACAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTg                                                                                                                                                    <  2:268333/249‑1 (MQ=255)
                                                                                     ccGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTc                                                                                                                                                  <  1:505107/251‑1 (MQ=255)
                                                                                     ccGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTc                                                                                                                                                  >  1:465644/1‑251 (MQ=255)
                                                                                      cGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCCTGCTCCCTGGCAGTACGTTGCACTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑TCCACC‑GGCCCTGCCCCGCTGACGCCCGCACCACCCAGAAGCAGCACTAACATATTGAGGACACAACGAAACCTCGTCAGGTTCCCTCCCCGCGTATGGGAGGTgaggtcc                                                                                                                                                 >  2:164414/1‑245 (MQ=255)
                                                                                       ggggAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGc                                                                                                                                                >  1:366981/1‑251 (MQ=255)
                                                                                                 agtgcgcgTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTg                                                                                                                                        >  1:319440/4‑249 (MQ=255)
                                                                                                 agtgcgcgTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTg                                                                                                                                        >  1:252803/4‑249 (MQ=255)
                                                                                                   tgcgcgTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGGAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGc                                                                                                                                     >  1:370654/2‑250 (MQ=255)
                                                                                                   tgcgcgTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGAACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGCTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGGGGTTGAGcg                                                                                                                                    >  2:311214/2‑250 (MQ=255)
                                                                                                   ggcgcgTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATTCTCACTGGCAGTACGTTGCGCTGATCGAAGCCGGCCCATAAGCGCGCCAGATGTGGCGCCAGCAGACCGATAAAC‑GCCCCC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCATCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCa                                                                                                                                    <  1:493240/250‑1 (MQ=255)
                                                                                                           cGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGCGGTTGAGCAGGTTCAGt                                                                                                                            >  1:110834/1‑251 (MQ=255)
                                                                                                              gTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGa                                                                                                                         <  1:123359/251‑1 (MQ=255)
                                                                                                                aCATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGAtt                                                                                                                       <  1:60617/251‑1 (MQ=255)
                                                                                                                    ccgacgtccgcagcggcgggccccccTGCAGAATGCTCTCGGGCAGTCCGTGGCTCTGATGAAACCCTCCCCAGAGGCGCGCTATATGTGGCACCAGAAGACAGATAAAC‑GCCACC‑GTACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGcc                                                                                                                   <  2:5984/230‑1 (MQ=255)
                                                                                                                      gCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGcc                                                                                                                   >  1:498380/1‑249 (MQ=255)
                                                                                                                        cagcagCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGCTGACCAAACGTAGCCTCGTCCGGTTCACCCCCAGCGTATGGGCGGTGCTGTCGCTTGGGTTGAGCAGGTTCAGTTGATTCGCcagc                                                                                                                >  2:18250/1‑250 (MQ=255)
                                                                                                                         agcagcATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGATGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGGTCAGTTGAATCGCcagc                                                                                                                >  2:243341/1‑249 (MQ=255)
                                                                                                                              caTCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAac                                                                                                          >  1:300725/1‑250 (MQ=255)
                                                                                                                                tCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAAACAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACg                                                                                                       >  1:65956/1‑251 (MQ=255)
                                                                                                                                 cAGCGTGGCCCCCAGCAGCAGGCGCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGa                                                                                                      <  2:319440/251‑1 (MQ=255)
                                                                                                                                      tGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGCCCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGAGCGGTCGCTGAGGTTGGGCAGGTTCAGATGCTTCGCCAGCAGCAACACGCCAggg                                                                                                 >  2:335122/1‑251 (MQ=255)
                                                                                                                                       ggCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAgggg                                                                                                >  2:235123/1‑250 (MQ=255)
                                                                                                                                        gCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGCTAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCAATCCCAGCGTAAGGGGGGTGCTGTCGTTGCGGTTGAGCAGGTTCAGTTGCTTCGCCAGCAGCAacccccccgggcc                                                                                               >  2:472665/1‑240 (MQ=255)
                                                                                                                                        gCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGAc                                                                                               <  1:113622/251‑1 (MQ=255)
                                                                                                                                         cccccAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAACGTATGGGCGGTGCTGTCGCTGAGGCTGTGCAGGGTCAGTTGATTCGCCAGCAGCAACACGACAggggcg                                                                                              >  1:138489/1‑248 (MQ=255)
                                                                                                                                          ccccAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGCTAAAC‑GCCACC‑GGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGCGGTTGAGCAGGGTCCGTTGATTCGCCAGCCGCAACCCGACAGGGAcgg                                                                                             >  2:186218/1‑249 (MQ=255)
                                                                                                                                              agcagcATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGt                                                                                         <  1:417566/251‑1 (MQ=255)
                                                                                                                                                cagcaTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGGCTGCAGTg                                                                                        >  1:455337/1‑250 (MQ=255)
                                                                                                                                                  gcaTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGacc                                                                                     <  1:152084/251‑1 (MQ=255)
                                                                                                                                                         cacgggcagtactgtgcgcggagagAAGCCTGGCCGGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGCCCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGaccaccac                                                                                <  2:310013/225‑1 (MQ=255)
                                                                                                                                                            tGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAACC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAagca                                                                                                                                                                                                                     >  1:221519/1‑113 (MQ=255)
                                                                                                                                                            tGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCAAGATGTGGCACCAGCAGACCGATAAAC‑GCCCCC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAagca                                                                                                                                                                                                                     <  2:221519/113‑1 (MQ=255)
                                                                                                                                                               cAGTACGTTGCGCTGATAGAAGCCTGCCAAGAAGCGCGCCAGATGTGGCCCCAGCAGACCGCTAACC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACAAGACACAGCACTAACATATTGATGACCAAACGTAGCCTCGTAAGGTTCACTCCCAGCGTATGGGCGGGGCAGTCGCTGAGGATGAGCAGGTTCAGTTGATTCGCACGCAGCAACACGACAGGGACTGCAGttgccacaacaggcaaa                                                                         >  2:244353/1‑234 (MQ=255)
                                                                                                                                                                            tGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATc                                                           >  1:100059/1‑251 (MQ=255)
                                                                                                                                                                             gATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCGCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACAt                                                            <  1:260891/249‑1 (MQ=255)
                                                                                                                                                                               ccGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTg                                                        <  1:189931/250‑1 (MQ=255)
                                                                                                                                                                                   aGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGcc                                                      >  1:141990/1‑249 (MQ=255)
                                                                                                                                                                                    gccgccccagagacgcgccgaaggtgGCAACAGCAGACCTATAACC‑GCCACC‑GCACCAGCCACGCTGACGCACGCACCAACCAGAAGCAGCCCTAACATATTGAGGACCAAACGTAGCCTCGTTAGGTTCACTCCCAGCTTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGc                                                   <  2:344933/228‑1 (MQ=255)
                                                                                                                                                                                       tGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCggg                                                >  1:412878/1‑251 (MQ=255)
                                                                                                                                                                                           cAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTg                                            <  1:235123/251‑1 (MQ=255)
                                                                                                                                                                                                   cgcCAGATGTGGCCCCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTc                                    <  1:360095/251‑1 (MQ=255)
                                                                                                                                                                                                      cAGATGTGGCTCCAGCAGACCGATAAAC‑GCCACC‑GGACCGGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGc                                 <  1:52706/251‑1 (MQ=255)
                                                                                                                                                                                                             ggcacaaccagACCGA‑AAACAGCCCCC‑GGACCTGCCACGCGGACGCCCGCACCAACCAGAAGCAGCATGAACATATTGAGGACCAAACGTAGCCGCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCGGCCAGCGGGCGTGGGACACTCCGCCTgccag                          <  2:377677/243‑1 (MQ=255)
                                                                                                                                                                                                                 cacgcgaaCCGATAACC‑CCCCCC‑GGCCCTCCCACGCCGACGCCCGCCCCAACCAGAAGCAGCACTAACATATGGACGCCCAACCGGAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGccagccag                      <  2:397692/244‑1 (MQ=255)
                                                                                                                                                                                                                     cagACCGATAAAC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGGGGGACACTCCGCCTGCCCGCCAGTaa                   >  2:102921/1‑250 (MQ=255)
                                                                                                                                                                                                                      agACCGATAAAC‑GCCACC‑GGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTaaaa                 <  1:18250/251‑1 (MQ=255)
                                                                                                                                                                                                                              aaaC‑GCCACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGt         >  1:473332/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cACC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCAt   <  1:91177/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aCC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTCGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAgg                                                                                                                                           >  2:454072/1‑114 (MQ=255)
                                                                                                                                                                                                                                      aCC‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAgg                                                                                                                                           <  1:454072/114‑1 (MQ=255)
                                                                                                                                                                                                                                       cccgaaCCTGCCACGCTGACGCACGCCCCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATg  <  1:491137/246‑1 (MQ=255)
                                                                                                                                                                                                                                       cc‑GGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGGTCAGTTGATTCGCCAGCAGCACCACGACCGGGACTGCAGTGACCACCACAGGCAAGAGCTGCCCGACATCCTGCCAGCCGGCGATGGACACGCCGCCTGCCAGACAGCAAAAGAtaacgaacacatg  >  2:104214/1‑238 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
AATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAT‑CAGCGCCAGCACTGC‑GCCTGCGGGCAGAT‑CT‑CCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAAC‑GCCACC‑GGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATG  >  NC_000913/4512406‑4512881

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