Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,739,601 T→C 13.8% E156G (GAG→GGG)  tyrA ← fused chorismate mutase T/prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,739,6010TC13.8% 143.0 / 2.3 60E156G (GAG→GGG) tyrAfused chorismate mutase T/prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (30/21);  new base C (5/4);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCAGAAACCATTGGTATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATT  >  NC_000913/2739383‑2739826
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   
gCTCCAGAAACCATTGGTATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCCGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGAAATTTGCCAAAAAATTACTccaaaaagtg                                                                                                                                                                                                                           >  2:92810/1‑220 (MQ=255)
gCTCCAGAAACCATTGGTATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACACGCATGGCCTGTAATGGCCCATTTTTCCCTCATGCCAGATCGACCAGAATAC‑AATCTTTCGTTAAAGTCGGATATTTGCCACAAAACTGCTCaagacctgg                                                                                                                                                                                                                           >  2:464603/1‑221 (MQ=255)
 gcccggAAACACTTGGTTTTCTTCCGTATTACGGCCAGCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACg                                                                                                                                                                                                                             <  2:79272/221‑1 (MQ=255)
    cAGAAACCTTTGGTATGTTTCCGGTTTTCGTCCATCCCCCCAGCCCGCAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGAtt                                                                                                                                                                                                                       <  2:162152/229‑1 (MQ=255)
     aGAAACCTTTGGTATCCTTCCGGTTTACGTCCATCACCCCAGACCACAACTTGCTTTGCCAGGCTTCCGCTTTCCGGACCGAACATCGGGTGTAGCCCCAGCCCCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTg                                                                                                                                                                                                                      <  2:383929/229‑1 (MQ=255)
      gAAACCATTGGTATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAAATTGCTCAGTAACGTGGATTgg                                                                                                                                                                                                                     >  1:219386/1‑229 (MQ=255)
        acccATTGTTTAGCTTCCGGTTCACCTCCATCCCACCAAGCCTCTACTTGCGTTTCCAGGCTACCGCTGGCCGGACCGAGCAGCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCCCCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGca                                                                                                                                                                                                                   <  2:361068/227‑1 (MQ=255)
        aaccagggtggtggttCCGGGTTACGGCCATCCCACCAGAACGAAACTTGCTTGGCCAGGCGTCCGCGGTCCGGAACGAACATCGGGTGTAACCCCAGCACCGGACCAGCAAGCGCCCCCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGTATTGGca                                                                                                                                                                                                                   <  2:404495/215‑1 (MQ=255)
                     cTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATTCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATCCCAGATCAAACCGAATAC‑AAGCTTTCGCGAAAGGCGGTAATTTTGCAATAACTTGCTCAATAACCGGGGTGGGGACACCaacaaaca                                                                                                                                                                                                       >  2:641947/1‑228 (MQ=255)
                     cTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCCTTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGCCGGTAATTTGCCAATAACTTGCCAATGAACGTGGCATGGCACAATaacaaaca                                                                                                                                                                                                       >  2:17986/1‑228 (MQ=255)
                     cTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGTTAAATGCGGTAATTTTCCAATAACTGGCTCAGTAACGGGGATAGGCTCACTAACAATCAc                                                                                                                                                                                                      >  1:385888/1‑229 (MQ=255)
                     cTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAAAGTGGATTGGCACACAAAAAATCAc                                                                                                                                                                                                      >  1:271982/1‑229 (MQ=255)
                               aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCa                                                                                                                                                                                            <  1:17986/229‑1 (MQ=255)
                                     ctccccccaaccctaaaCTTGCTTTTCCAGGCTACCGCCTTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGca                                                                                                                                                                                      <  2:322331/214‑1 (MQ=255)
                                      tCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCCGTAAAGGCGGTACATTGCCAATAACTTTCCCAGCAACGTGGATTGGCACACTACAAATCACCCTTCCCGCATCGGcaa                                                                                                                                                                                     >  2:613984/1‑229 (MQ=255)
                                       cacaCCAGCCCGCAACTTGTTTTGCCAGGGTACAGCTTCCCGGACCGAACCTCGGGTGTAGCCCCGGCCCCTGACCAGCATGCGCCACCAGCATGGCGTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTTGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGcaa                                                                                                                                                                                     <  2:690118/228‑1 (MQ=255)
                                              gACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATACCTTGCCCAGTAACATGGATTGGCACCCTCACAATCACCATTCCGGCATCGGCAACA‑ATATc                                                                                                                                                                             >  1:585371/1‑228 (MQ=255)
                                                 cacaACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAAACTTCTCAGTAACGTGGGTTGGCACACTAACAATCAACATTCCGGCATCGGCAAAA‑atataacac                                                                                                                                                                         >  2:132132/1‑224 (MQ=255)
                                                   caACTTGCTTTGCCAGGCTACCGGTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGc                                                                                                                                                                       <  1:191945/229‑1 (MQ=255)
                                                    accgtGCCTTGCCAGGCTACCGCTGTCCGGACCGAACATGGGGTGTAGCCCCAGCAACGGAACATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGACAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCt                                                                                                                                                                      <  1:689473/225‑1 (MQ=255)
                                                    aaCTTGGCTTGCCATGCTACCGCTGGCCGGACCGGACACAGGTGGTAGCCCCAGGACCGGGCCATCATGCGCCTCCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAATTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCt                                                                                                                                                                      <  1:393801/229‑1 (MQ=255)
                                                      cTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGTCCAGTAACATGGATTGGCCCACTACAAATCACCATTCCGGCATCGGCCAAC‑ATATCAGCCGCCCg                                                                                                                                                                    >  1:733141/1‑229 (MQ=255)
                                                              gCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCCGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCTACA‑AAATCAGCCGCTCGacaccat                                                                                                                                                             >  1:690630/1‑222 (MQ=255)
                                                              gCCAGGCTACCGCTGTCCGGACCGAACATCGAGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATCGCCTGTAATGGCCCATTTTTTACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGGTCAATAACATGGATTGGCACACTAACAAACACCATTCCCGGATCGCCAACA‑AAACAAACCGCTCGCTCCCAg                                                                                                                                                             >  1:656483/1‑228 (MQ=255)
                                                                  ggCTACCGGTTTCCCGACCGAACATCGGGTTTACCCCAACGACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATg                                                                                                                                                        <  1:259649/229‑1 (MQ=255)
                                                                  ggCTACCGCTTTCCGGACCGAACAGCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATg                                                                                                                                                        <  1:480929/229‑1 (MQ=255)
                                                                             tCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGTGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGACCCCTAACAATCACCATTCCGTCCTCGGCAACA‑ATATCAGCCCCTCGCTccacgccaagtgactccaaaa                                                                                                                                             >  2:126489/1‑210 (MQ=255)
                                                                             tCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGGAACGGTGATTGGGACACTAACAATAACCATTCCCGGATC‑GCAACACAAATAAGCCGCTCTCACCCAGCCATATTTCGCCaaaa                                                                                                                                             >  2:397951/1‑229 (MQ=255)
                                                                                  aCCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCACCA‑ATATCAGCCGCCCGATCCCAGTCATGTTGCTCCAGAATCCg                                                                                                                                         >  1:43919/1‑228 (MQ=255)
                                                                                       aCATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCCGTCATGTTGCTCCAGAATCCCCACCTg                                                                                                                                   >  1:411437/1‑229 (MQ=255)
                                                                                       aCATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCAAGTTGCTCCCGAATCCGCACCTg                                                                                                                                   >  1:512100/1‑229 (MQ=255)
                                                                                          tCGGGTGGACCCCCAGCCCCGGCCCATCCTGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATa                                                                                                                                <  2:231212/229‑1 (MQ=255)
                                                                                                   gCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTTACAATAACTTGCTCAGTAACGTGGATTGGCACCATACCAATCACCATTCCGCGCTCAGCAACA‑ATATCCACCGCCCGGACCCCGTCACGTTCCTCCACAAACCGCAcccggaaaacccaga                                                                                                                        >  2:674161/1‑215 (MQ=255)
                                                                                                      cgcgCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGt                                                                                                                    <  1:690142/226‑1 (MQ=255)
                                                                                                      ccggctccGGACCATCCTGCGCACCCAGCATGGCCTGTAATGGCCCTTTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGt                                                                                                                    <  2:264581/223‑1 (MQ=255)
                                                                                                      cccgccccGGACCTTCATGCGCCACCAGCATGGGCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCCTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGt                                                                                                                    <  2:96086/223‑1 (MQ=255)
                                                                                                       cAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTTCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCCGCCCCCCGGTCCCAGTCCTGTTGCTCCAGAAACCGCCCCTGAAAAACCCAAAggag                                                                                                                    >  2:101006/1‑226 (MQ=255)
                                                                                                        aGCACCGGACCATCATCCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTCTCCAATAAATTGACCAGTAACAGTGATTGGCACACCAACAATCACCCTACCGGCCTCAGCAACA‑AAAATAGCCGACCCATCCCAGTCAAGTTATTCCAGAATCCGAAcaccatccacccaaagaggac                                                                                                                  >  2:507622/1‑209 (MQ=255)
                                                                                                         gCACCGGACCATCATGTGCCACCAGAATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAA‑AATCTTTCGGTAAAGGCCGTAATTTGCCAATAACTTGTTCAGAAACTTCAATTTGCCCCCTCAAAATAACCATTTCGCGATCGGCAACA‑AAATAAACCGCCCACACCCAGTCAAGTTGCCCCAGAAACCCCACCCAAAAACCCCagagagcaag                                                                                                                 >  2:39035/1‑223 (MQ=255)
                                                                                                               gACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCACCA‑ATATAAACCGCCCCCTCCCCCCCATGTTGCCCCACAAACCGCACCCCATAAACCGAAAGGGCCAGCAtct                                                                                                            >  2:271037/1‑228 (MQ=255)
                                                                                                               gACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAAAGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATAAGCCGCCCGATCCCCATCATGTTGCCCCAGAATCCACACCTGAAAACCCGGGAGGGCCAGCAtcc                                                                                                            >  2:54623/1‑227 (MQ=255)
                                                                                                                       tggcccaccgataTGGCCTGTTATGGCCCATTTTGCACTGATGCCAAATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCCCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACa                                                                                                  <  2:271982/217‑1 (MQ=255)
                                                                                                                       tGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCCa‑aa                                                                                                   >  1:642232/1‑227 (MQ=255)
                                                                                                                            ccccaGCATGGCCTGTAATGGTCCATTTTTCACTGAAGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCg                                                                                              <  1:132132/226‑1 (MQ=255)
                                                                                                                            caccaGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAAGCg                                                                                              >  1:647547/1‑228 (MQ=255)
                                                                                                                                         tGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGGAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGaccgc                                                                                >  1:311225/1‑229 (MQ=255)
                                                                                                                                                ggccATTTTTTCTTGATGCCAGA‑CTACCCAAAAACAAATCGTTCGGTAAGGGTGGTAATTTGCCAATAACGTGCCCATAAACGGGGATGGGCACCCAAACAAGCCCCATTCCGGCAGCGGAAACA‑ATACCAGCCCCTCGATCCCAGTCATGTTGCCCCAGAATCCGCACCGGATACCCCGAGAGGGTCAGCATCTTCCCGA‑ACAGGCGTCCCATCTGaccgccaccgc                                                                          <  2:43919/226‑1 (MQ=255)
                                                                                                                                                 cccATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAAAAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGGTAACCCGAGAGGGTCAGCCTCTTCTCGA‑ACAGGCGTCCCATCTGacccccaccccc                                                                         >  1:711349/1‑228 (MQ=255)
                                                                                                                                                     tttttCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCCCGATCCCAGTCATGTTGCTCCAGAATCCGAACCTGATAACCCCAGAAGGTCAGCATCTTCTCCA‑AAAGGCGTCCCATCTAcccgcccccgcccccga                                                                    >  1:736454/1‑224 (MQ=255)
                                                                                                                                                                tGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAAAAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGACCCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGGGTCCCATCTGACCCCAACCGCCGACGATAAACACCgac                                                         >  1:736554/1‑227 (MQ=255)
                                                                                                                                                                         gACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAAAATCCCCACCTGCAAAACCGCAAGGGGCCGCATCTTCTCGACACA‑GCGTCCCACCTGAACGCCACCGCCGCCGGTAACCCACGGACACAgaga                                                 >  2:326181/1‑228 (MQ=255)
                                                                                                                                                                                           tCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGCCAAACGGTTTTaaa                               >  1:225850/1‑229 (MQ=255)
                                                                                                                                                                                               tAAAGCGGGTAATTTGCAAATAACTTGCTCAGATACGTGGATTGGCACACTAACAACCCCCATTCCGGCATCGGAAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGTCCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCt                           <  1:641948/229‑1 (MQ=255)
                                                                                                                                                                                               aaaaGGCGGTAATTTGCAAATAACTTGCCCAGTAACGTGGATTGGCCCACTAACAATCCCCATTCCGGCATGGGCAACA‑ATACCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCCCGA‑ACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCt                           <  1:464603/228‑1 (MQ=255)
                                                                                                                                                                                                  aGGCGGTAATTTGCCAATACCTTGCTAAGTAACGGGGATTGGCCAACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCCCCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTg                        <  1:92810/229‑1 (MQ=255)
                                                                                                                                                                                                            ttGCCAATAACTTGCTCAGTACCGGGGATTGGCACCCTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGACCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCAc              <  1:101006/229‑1 (MQ=255)
                                                                                                                                                                                                            ttGCCAATAACTTGCCTAGTAACGGGGATGGGCACACTAACAATCACCATTCCGGCACCGGCAACA‑ATACCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCAc              <  1:5150/229‑1 (MQ=255)
                                                                                                                                                                                                              gCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACAGGACGCAGTGCCGGACAAAGGGTTTTAAATCCTTTGTCGttttacctg            >  2:730351/1‑224 (MQ=255)
                                                                                                                                                                                                                  aTAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGGTTTAAATCCTTTGCCGTTTTCACTggcg         >  1:296778/1‑226 (MQ=255)
                                                                                                                                                                                                                         gCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGAtt  >  1:276710/1‑228 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   
GCTCCAGAAACCATTGGTATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAC‑AATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACA‑ATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGA‑ACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATT  >  NC_000913/2739383‑2739826

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