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

GCTCCAGAAACCATTGGTATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATT  >  NC_000913/2739383‑2739826
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   
gCTCCAGAAACCATTGGTATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCCGATCGACCAGAATACAATCTTTCGGTAAAGGCGGAAATTTGCCAAAAAATTACTccaaaaagtg                                                                                                                                                                                                                           >  2:92810/1‑220 (MQ=255)
gCTCCAGAAACCATTGGTATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACACGCATGGCCTGTAATGGCCCATTTTTCCCTCATGCCAGATCGACCAGAATACAATCTTTCGTTAAAGTCGGATATTTGCCACAAAACTGCTCaagacctgg                                                                                                                                                                                                                           >  2:464603/1‑221 (MQ=255)
 gcccggAAACACTTGGTTTTCTTCCGTATTACGGCCAGCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACg                                                                                                                                                                                                                             <  2:79272/221‑1 (MQ=255)
    cAGAAACCTTTGGTATGTTTCCGGTTTTCGTCCATCCCCCCAGCCCGCAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGAtt                                                                                                                                                                                                                       <  2:162152/229‑1 (MQ=255)
     aGAAACCTTTGGTATCCTTCCGGTTTACGTCCATCACCCCAGACCACAACTTGCTTTGCCAGGCTTCCGCTTTCCGGACCGAACATCGGGTGTAGCCCCAGCCCCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTg                                                                                                                                                                                                                      <  2:383929/229‑1 (MQ=255)
      gAAACCATTGGTATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAAATTGCTCAGTAACGTGGATTgg                                                                                                                                                                                                                     >  1:219386/1‑229 (MQ=255)
        acccATTGTTTAGCTTCCGGTTCACCTCCATCCCACCAAGCCTCTACTTGCGTTTCCAGGCTACCGCTGGCCGGACCGAGCAGCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCCCCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGca                                                                                                                                                                                                                   <  2:361068/227‑1 (MQ=255)
        aaccagggtggtggttCCGGGTTACGGCCATCCCACCAGAACGAAACTTGCTTGGCCAGGCGTCCGCGGTCCGGAACGAACATCGGGTGTAACCCCAGCACCGGACCAGCAAGCGCCCCCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGTATTGGca                                                                                                                                                                                                                   <  2:404495/215‑1 (MQ=255)
                     cTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATTCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATCCCAGATCAAACCGAATACAAGCTTTCGCGAAAGGCGGTAATTTTGCAATAACTTGCTCAATAACCGGGGTGGGGACACCaacaaaca                                                                                                                                                                                                       >  2:641947/1‑228 (MQ=255)
                     cTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCCTTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGCCGGTAATTTGCCAATAACTTGCCAATGAACGTGGCATGGCACAATaacaaaca                                                                                                                                                                                                       >  2:17986/1‑228 (MQ=255)
                     cTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGTTAAATGCGGTAATTTTCCAATAACTGGCTCAGTAACGGGGATAGGCTCACTAACAATCAc                                                                                                                                                                                                      >  1:385888/1‑229 (MQ=255)
                     cTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAAAGTGGATTGGCACACAAAAAATCAc                                                                                                                                                                                                      >  1:271982/1‑229 (MQ=255)
                               aCGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCa                                                                                                                                                                                            <  1:17986/229‑1 (MQ=255)
                                     ctccccccaaccctaaaCTTGCTTTTCCAGGCTACCGCCTTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGca                                                                                                                                                                                      <  2:322331/214‑1 (MQ=255)
                                      tCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCCGTAAAGGCGGTACATTGCCAATAACTTTCCCAGCAACGTGGATTGGCACACTACAAATCACCCTTCCCGCATCGGcaa                                                                                                                                                                                     >  2:613984/1‑229 (MQ=255)
                                       cacaCCAGCCCGCAACTTGTTTTGCCAGGGTACAGCTTCCCGGACCGAACCTCGGGTGTAGCCCCGGCCCCTGACCAGCATGCGCCACCAGCATGGCGTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTTGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGcaa                                                                                                                                                                                     <  2:690118/228‑1 (MQ=255)
                                              gACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATACCTTGCCCAGTAACATGGATTGGCACCCTCACAATCACCATTCCGGCATCGGCAACAATATc                                                                                                                                                                             >  1:585371/1‑228 (MQ=255)
                                                 cacaACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAAACTTCTCAGTAACGTGGGTTGGCACACTAACAATCAACATTCCGGCATCGGCAAAAatataacac                                                                                                                                                                         >  2:132132/1‑224 (MQ=255)
                                                   caACTTGCTTTGCCAGGCTACCGGTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGc                                                                                                                                                                       <  1:191945/229‑1 (MQ=255)
                                                    accgtGCCTTGCCAGGCTACCGCTGTCCGGACCGAACATGGGGTGTAGCCCCAGCAACGGAACATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGACAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCt                                                                                                                                                                      <  1:689473/225‑1 (MQ=255)
                                                    aaCTTGGCTTGCCATGCTACCGCTGGCCGGACCGGACACAGGTGGTAGCCCCAGGACCGGGCCATCATGCGCCTCCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAATTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCt                                                                                                                                                                      <  1:393801/229‑1 (MQ=255)
                                                      cTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGTCCAGTAACATGGATTGGCCCACTACAAATCACCATTCCGGCATCGGCCAACATATCAGCCGCCCg                                                                                                                                                                    >  1:733141/1‑229 (MQ=255)
                                                              gCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCCGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCTACAAAATCAGCCGCTCGacaccat                                                                                                                                                             >  1:690630/1‑222 (MQ=255)
                                                              gCCAGGCTACCGCTGTCCGGACCGAACATCGAGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATCGCCTGTAATGGCCCATTTTTTACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGGTCAATAACATGGATTGGCACACTAACAAACACCATTCCCGGATCGCCAACAAAACAAACCGCTCGCTCCCAg                                                                                                                                                             >  1:656483/1‑228 (MQ=255)
                                                                  ggCTACCGGTTTCCCGACCGAACATCGGGTTTACCCCAACGACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATg                                                                                                                                                        <  1:259649/229‑1 (MQ=255)
                                                                  ggCTACCGCTTTCCGGACCGAACAGCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATg                                                                                                                                                        <  1:480929/229‑1 (MQ=255)
                                                                             tCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGTGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGACCCCTAACAATCACCATTCCGTCCTCGGCAACAATATCAGCCCCTCGCTccacgccaagtgactccaaaa                                                                                                                                             >  2:126489/1‑210 (MQ=255)
                                                                             tCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGGAACGGTGATTGGGACACTAACAATAACCATTCCCGGATCGCAACACAAATAAGCCGCTCTCACCCAGCCATATTTCGCCaaaa                                                                                                                                             >  2:397951/1‑229 (MQ=255)
                                                                                  aCCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCACCAATATCAGCCGCCCGATCCCAGTCATGTTGCTCCAGAATCCg                                                                                                                                         >  1:43919/1‑228 (MQ=255)
                                                                                       aCATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCCGTCATGTTGCTCCAGAATCCCCACCTg                                                                                                                                   >  1:411437/1‑229 (MQ=255)
                                                                                       aCATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCAAGTTGCTCCCGAATCCGCACCTg                                                                                                                                   >  1:512100/1‑229 (MQ=255)
                                                                                          tCGGGTGGACCCCCAGCCCCGGCCCATCCTGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATa                                                                                                                                <  2:231212/229‑1 (MQ=255)
                                                                                                   gCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTTACAATAACTTGCTCAGTAACGTGGATTGGCACCATACCAATCACCATTCCGCGCTCAGCAACAATATCCACCGCCCGGACCCCGTCACGTTCCTCCACAAACCGCAcccggaaaacccaga                                                                                                                        >  2:674161/1‑215 (MQ=255)
                                                                                                      cgcgCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGt                                                                                                                    <  1:690142/226‑1 (MQ=255)
                                                                                                      ccggctccGGACCATCCTGCGCACCCAGCATGGCCTGTAATGGCCCTTTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGt                                                                                                                    <  2:264581/223‑1 (MQ=255)
                                                                                                      cccgccccGGACCTTCATGCGCCACCAGCATGGGCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCCTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGt                                                                                                                    <  2:96086/223‑1 (MQ=255)
                                                                                                       cAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTTCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCCGCCCCCCGGTCCCAGTCCTGTTGCTCCAGAAACCGCCCCTGAAAAACCCAAAggag                                                                                                                    >  2:101006/1‑226 (MQ=255)
                                                                                                        aGCACCGGACCATCATCCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTCTCCAATAAATTGACCAGTAACAGTGATTGGCACACCAACAATCACCCTACCGGCCTCAGCAACAAAAATAGCCGACCCATCCCAGTCAAGTTATTCCAGAATCCGAAcaccatccacccaaagaggac                                                                                                                  >  2:507622/1‑209 (MQ=255)
                                                                                                         gCACCGGACCATCATGTGCCACCAGAATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATAAAATCTTTCGGTAAAGGCCGTAATTTGCCAATAACTTGTTCAGAAACTTCAATTTGCCCCCTCAAAATAACCATTTCGCGATCGGCAACAAAATAAACCGCCCACACCCAGTCAAGTTGCCCCAGAAACCCCACCCAAAAACCCCagagagcaag                                                                                                                 >  2:39035/1‑223 (MQ=255)
                                                                                                               gACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCACCAATATAAACCGCCCCCTCCCCCCCATGTTGCCCCACAAACCGCACCCCATAAACCGAAAGGGCCAGCAtct                                                                                                            >  2:271037/1‑228 (MQ=255)
                                                                                                               gACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAAAGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATAAGCCGCCCGATCCCCATCATGTTGCCCCAGAATCCACACCTGAAAACCCGGGAGGGCCAGCAtcc                                                                                                            >  2:54623/1‑227 (MQ=255)
                                                                                                                       tggcccaccgataTGGCCTGTTATGGCCCATTTTGCACTGATGCCAAATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCCCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACa                                                                                                  <  2:271982/217‑1 (MQ=255)
                                                                                                                       tGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCCaaa                                                                                                   >  1:642232/1‑227 (MQ=255)
                                                                                                                            ccccaGCATGGCCTGTAATGGTCCATTTTTCACTGAAGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCg                                                                                              <  1:132132/226‑1 (MQ=255)
                                                                                                                            caccaGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAAGCg                                                                                              >  1:647547/1‑228 (MQ=255)
                                                                                                                                         tGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGGAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGaccgc                                                                                >  1:311225/1‑229 (MQ=255)
                                                                                                                                                ggccATTTTTTCTTGATGCCAGACTACCCAAAAACAAATCGTTCGGTAAGGGTGGTAATTTGCCAATAACGTGCCCATAAACGGGGATGGGCACCCAAACAAGCCCCATTCCGGCAGCGGAAACAATACCAGCCCCTCGATCCCAGTCATGTTGCCCCAGAATCCGCACCGGATACCCCGAGAGGGTCAGCATCTTCCCGAACAGGCGTCCCATCTGaccgccaccgc                                                                          <  2:43919/226‑1 (MQ=255)
                                                                                                                                                 cccATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAAAAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGGTAACCCGAGAGGGTCAGCCTCTTCTCGAACAGGCGTCCCATCTGacccccaccccc                                                                         >  1:711349/1‑228 (MQ=255)
                                                                                                                                                     tttttCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCCCGATCCCAGTCATGTTGCTCCAGAATCCGAACCTGATAACCCCAGAAGGTCAGCATCTTCTCCAAAAGGCGTCCCATCTAcccgcccccgcccccga                                                                    >  1:736454/1‑224 (MQ=255)
                                                                                                                                                                tGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAAAAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGACCCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGGGTCCCATCTGACCCCAACCGCCGACGATAAACACCgac                                                         >  1:736554/1‑227 (MQ=255)
                                                                                                                                                                         gACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAAAATCCCCACCTGCAAAACCGCAAGGGGCCGCATCTTCTCGACACAGCGTCCCACCTGAACGCCACCGCCGCCGGTAACCCACGGACACAgaga                                                 >  2:326181/1‑228 (MQ=255)
                                                                                                                                                                                           tCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGCCAAACGGTTTTaaa                               >  1:225850/1‑229 (MQ=255)
                                                                                                                                                                                               tAAAGCGGGTAATTTGCAAATAACTTGCTCAGATACGTGGATTGGCACACTAACAACCCCCATTCCGGCATCGGAAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGTCCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCt                           <  1:641948/229‑1 (MQ=255)
                                                                                                                                                                                               aaaaGGCGGTAATTTGCAAATAACTTGCCCAGTAACGTGGATTGGCCCACTAACAATCCCCATTCCGGCATGGGCAACAATACCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCCCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCt                           <  1:464603/228‑1 (MQ=255)
                                                                                                                                                                                                  aGGCGGTAATTTGCCAATACCTTGCTAAGTAACGGGGATTGGCCAACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCCCCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTg                        <  1:92810/229‑1 (MQ=255)
                                                                                                                                                                                                            ttGCCAATAACTTGCTCAGTACCGGGGATTGGCACCCTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGACCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCAc              <  1:101006/229‑1 (MQ=255)
                                                                                                                                                                                                            ttGCCAATAACTTGCCTAGTAACGGGGATGGGCACACTAACAATCACCATTCCGGCACCGGCAACAATACCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCAc              <  1:5150/229‑1 (MQ=255)
                                                                                                                                                                                                              gCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACAGGACGCAGTGCCGGACAAAGGGTTTTAAATCCTTTGTCGttttacctg            >  2:730351/1‑224 (MQ=255)
                                                                                                                                                                                                                  aTAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGGTTTAAATCCTTTGCCGTTTTCACTggcg         >  1:296778/1‑226 (MQ=255)
                                                                                                                                                                                                                         gCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGAtt  >  1:276710/1‑228 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   
GCTCCAGAAACCATTGGTATGCTTCCGGTTTACGTCCATCACACCAGACCACAACTTGCTTTGCCAGGCTACCGCTGTCCGGACCGAACATCGGGTGTAGCCCCAGCACCGGACCATCATGCGCCACCAGCATGGCCTGTAATGGCCCATTTTTCACTGATGCCAGATCGACCAGAATACAATCTTTCGGTAAAGGCGGTAATTTGCCAATAACTTGCTCAGTAACGTGGATTGGCACACTAACAATCACCATTCCGGCATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATT  >  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: