Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,011,944 A→T 5.8% D290V (GAT→GTT) 
D291E (GAT→GAA
yigM →
metR ←
putative inner membrane EamA‑like transporter
methionine biosynthesis regulon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,011,9440AT 5.8% 235.9 / 3.4 69D290V (GAT→GTT) 
D291E (GAT→GAA
yigM
metR
putative inner membrane EamA‑like transporter
methionine biosynthesis regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (32/33);  new base T (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGA  >  NC_000913/4011681‑4012226
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                           
gCAATGGGGCGTGCTGGTGTTGCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAg                                                                                                                                                                                                                                                                     <  2:189332/288‑1 (MQ=255)
         cATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACGTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAAc                                                                                                                                                                                                                                                            <  1:262084/288‑1 (MQ=255)
           ttCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCCCGCCGGTGCTCCTGGCCTCACCGTGGGTGCCTCGTAAGTGGGTTCCCACCGGCCCTTCACAACCGGCACATGATCTCAGGCTTGATTCGGCGCTGTGCGAATACACGc                                                                                                                                                                                                                                                          >  2:262084/1‑288 (MQ=255)
            gcTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCtt                                                                                                                                                                                                                                                        <  2:475681/288‑1 (MQ=255)
                  gtTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTaa                                                                                                                                                                                                                                                  <  2:268541/289‑1 (MQ=255)
                      cTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTCCAGGCGCGCTGGTGATCCTGGCCTCACTGTTGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGGGATTGCGCGCTGCGCGAATAAACGCTTCAGTACCTg                                                                                                                                                                                                                                               >  2:422159/1‑288 (MQ=255)
                      cTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATaacc                                                                                                                                                                                                                                                            >  2:105394/1‑273 (MQ=255)
                      cTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAAc                                                                                                                                                                                                                                                            <  1:105394/275‑1 (MQ=255)
                          gCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTg                                                                                                                                                                                                                                           >  1:139454/1‑288 (MQ=255)
                               gtggCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACgct                                                                                                                                                                                                                                      >  1:271692/1‑288 (MQ=255)
                                                 tatttgttggggAACTACGGCGCGACGCAGGTGGACGCGGGAACGCGGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGTCCAACGTTTATTACAGGCGCGCTGGTGATCCGGGCCTCACTGTGGGGGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCAcgg                                                                                                                                                                                                                   <  2:484799/278‑1 (MQ=255)
                                                    tttATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACTCTGCCCGCCCTCGCGCACggcggc                                                                                                                                                                                                               >  2:432872/1‑289 (MQ=255)
                                                     ttATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGTTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGgcggcg                                                                                                                                                                                                              <  2:114683/289‑1 (MQ=255)
                                                      tATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGt                                                                                                                                                                                                             >  1:533069/1‑289 (MQ=255)
                                                              aCTCCGGCGGGACGCAGGTGGAGGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCgg                                                                                                                                                                                                     <  2:525184/289‑1 (MQ=255)
                                                                  cGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGAGGAGCGCAGGCCGGATTGCGCGCTGTGCGAATCAACGCTTCCGTAACTGGCTGACGCTTCTACCCTGCGCGCCCGGCGGCGTaaacgcggccc                                                                                                                                                                                                  >  2:231053/1‑278 (MQ=255)
                                                                  cGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCAACAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGTTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGCTGATCGCAGGCGTGATTGCGCGCTGATGGAACAAACGCATCCGTAAACGGCTGACGCTGCTTGCCATGCACACCGGCGGCTGACAGTGGGccc                                                                                                                                                                                                  >  2:88745/1‑286 (MQ=255)
                                                                     cgcgACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCcac                                                                                                                                                                                               <  2:65260/288‑1 (MQ=255)
                                                                       cgACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTTGCTGCCGCTGCTCGCCATCGCGCACGGCGGCGTACCGTCGGCCCCACaa                                                                                                                                                                                             >  2:234176/1‑288 (MQ=255)
                                                                          cGCAGGTGGGAGCCGGAAGGCTGGGCGTTATGAATAATAGGAACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTAGTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGcc                                                                                                                                                                                          <  2:271692/288‑1 (MQ=255)
                                                                                tGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGcc                                                                                                                                                                                    <  1:215049/288‑1 (MQ=255)
                                                                                    cGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAggg                                                                                                                                                                               >  1:86672/1‑289 (MQ=255)
                                                                                       cGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCtt                                                                                                                                                                                        >  2:547298/1‑277 (MQ=255)
                                                                                       cGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCtt                                                                                                                                                                                        <  1:547298/277‑1 (MQ=255)
                                                                                       cGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGttt                                                                                                                                                                            <  2:326980/289‑1 (MQ=255)
                                                                                       cGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGtt                                                                                                                                                                             <  2:49012/288‑1 (MQ=255)
                                                                                           aCGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGt                                                                                                                                                                        <  1:422159/289‑1 (MQ=255)
                                                                                           aCGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGt                                                                                                                                                                        <  1:432629/289‑1 (MQ=255)
                                                                                             gCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTca                                                                                                                                                                      <  2:86672/289‑1 (MQ=255)
                                                                                               tGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTcac                                                                                                                                                                     <  1:231053/288‑1 (MQ=255)
                                                                                                         gAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGAc                                                                                                                                                             >  1:528032/1‑286 (MQ=255)
                                                                                                           ataataTGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTg                                                                                                                                                         >  1:95654/1‑288 (MQ=255)
                                                                                                            taagaTGCGCGTTCCGGCAGGGCTGCTGGTAAACCTGGCTAGCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGgc                                                                                                                                                       <  2:426142/289‑1 (MQ=255)
                                                                                                                 tGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCa                                                                                                                                                  <  1:88745/289‑1 (MQ=255)
                                                                                                                     cGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCAGGGTTTTGGTCACCACCAGACCCTGGCGCCCAAAAc                                                                                                                                              >  1:291225/1‑288 (MQ=255)
                                                                                                                      gTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAAct                                                                                                                                             <  1:328647/289‑1 (MQ=255)
                                                                                                                       ttCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAAct                                                                                                                                             <  2:533069/288‑1 (MQ=255)
                                                                                                                        tCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTc                                                                                                                                                                       >  2:304967/1‑261 (MQ=255)
                                                                                                                        tCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTc                                                                                                                                                                       <  1:304967/261‑1 (MQ=255)
                                                                                                                                ggCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTTACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCtcctac                                                                                                                                    >  1:406122/1‑288 (MQ=255)
                                                                                                                                ggCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCcaca                                                                                                                                                                                              <  1:318548/230‑1 (MQ=255)
                                                                                                                                ggCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCcaca                                                                                                                                                                                              >  2:318548/1‑230 (MQ=255)
                                                                                                                                     ctgGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGTTCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAAt                                                                                                                               >  2:190212/1‑288 (MQ=255)
                                                                                                                                          aaaCCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTa                                                                                                                         >  1:85379/1‑289 (MQ=255)
                                                                                                                                             ccTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAgcg                                                                                                                      <  1:537260/289‑1 (MQ=255)
                                                                                                                                                   tATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCgg                                                                                                                                                                                                     <  2:258420/204‑1 (MQ=255)
                                                                                                                                                   tATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCgg                                                                                                                                                                                                     >  1:258420/1‑204 (MQ=255)
                                                                                                                                                                  cgcacacgagCCAAGGTTTATTCCCGGCGCGTGGGTGTTGCTGCACTCACTTTGGGTGCAGCTTAGGTGGTTCACTCCGCGCTCTTCACAAGCGGCAAATGATCGCAGGTGTGTTTGCGCGCTGAGCGAATAAACGCTTCCGGAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCcgcgcg                                                                                                  <  2:231294/279‑1 (MQ=255)
                                                                                                                                                                                   ttattaCAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGTTCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATcaa                                                                                 <  1:190212/288‑1 (MQ=255)
                                                                                                                                                                                     attaCAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCa                                                                                                                                                                                                                       <  2:178947/153‑1 (MQ=255)
                                                                                                                                                                                     attaCAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCa                                                                                                                                                                                                                       >  1:178947/1‑153 (MQ=255)
                                                                                                                                                                                     attaCAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATcaacaa                                                                              >  1:179000/1‑289 (MQ=255)
                                                                                                                                                                                      ttaCAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAAt                                                                             <  1:145113/289‑1 (MQ=255)
                                                                                                                                                                                            gcgcgcTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCCTCTGAATCAACAATAAGgtt                                                                       >  1:521271/1‑288 (MQ=255)
                                                                                                                                                                                                              ccTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCcgcg                                                                                                    >  1:63145/1‑242 (MQ=255)
                                                                                                                                                                                                              ccTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCcgcg                                                                                                    <  2:63145/242‑1 (MQ=255)
                                                                                                                                                                                                              ccTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAAt                                                                             >  2:512873/1‑265 (MQ=255)
                                                                                                                                                                                                              ccTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAAt                                                                             <  1:512873/265‑1 (MQ=255)
                                                                                                                                                                                                              ccTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGt                                                     >  1:74011/1‑289 (MQ=255)
                                                                                                                                                                                                                tCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGa                                                   <  1:432872/289‑1 (MQ=255)
                                                                                                                                                                                                                  aCTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCCCAAAACTCTCCACTACCCAATGCGGTCGCACGGCAATACCCATCCGCCCGGCAACCATCTGAATCCACAATAAGGTGTTTACCACCCTTTTCAGTGACg                                                 >  2:172253/1‑289 (MQ=255)
                                                                                                                                                                                                                            cATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGc                                        >  1:175795/1‑288 (MQ=255)
                                                                                                                                                                                                                                       ggTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTgaa                            >  2:73706/1‑289 (MQ=255)
                                                                                                                                                                                                                                          tCGCTCCGCGCTCTTCACAAACGGCAGATGTTCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGAcc                                                                                                                                                            >  2:41435/1‑159 (MQ=255)
                                                                                                                                                                                                                                          tCGCTCCGCGCTCTTCACAAACGGCAGATGTTCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGAcc                                                                                                                                                            <  1:41435/159‑1 (MQ=255)
                                                                                                                                                                                                                                                             acgtgCAGATGATCGCGGCCGCGGTTGTGCGCTGAGCGAATAAACGCTTCCGCATATCGCTGACGCTGTTCGCCACAGCGCACGGGGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCa       <  2:528032/285‑1 (MQ=255)
                                                                                                                                                                                                                                                              aCGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAg      >  1:445429/1‑289 (MQ=255)
                                                                                                                                                                                                                                                               cGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAg      >  2:365586/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                 gCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATCCCCATCCGCGCGGCAACCATCTGAATCAACAATACGGTGTTATCGACGCTTTTACGTGACGGGCCGACGCCTGCCGGCCGAAGCAAATGACGCCCGCCATCCacttc   >  2:246558/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                   ggatgagCGCAGGCGGGATTGCGCGCTGGGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGa  <  2:139454/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                           
GCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGCCCAGGGTTTTGGTCACCACCAGACCCTGGCGCTCAAAACTCTCTACTACCCAATGCGGTAGCGCGGCAATACCCATCCGCGCGGCAACCATCTGAATCAACAATAAGGTGTTATCGACGCTTTTCAGTGACGGGCTGACGCCTGCCGGCTGAAGAAAATGCCGCCAGACATCCAGTCGA  >  NC_000913/4011681‑4012226

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: