Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,113,957 T→A 16.4% L132F (TTA→TTT pppA ← bifunctional prepilin leader peptidase/ methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,113,9570TA16.4% 130.0 / 2.3 62L132F (TTA→TTTpppAbifunctional prepilin leader peptidase/ methylase
Reads supporting (aligned to +/- strand):  ref base T (40/9);  new base A (5/6);  total (47/15)
Fisher's exact test for biased strand distribution p-value = 2.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTG  >  NC_000913/3113731‑3114193
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
cGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGGTATCGAGGTCGATGACACTCGCGGCACTCACCCAGCCGGGTACTATCATCAc                                                                                                                                                                                                                                   >  1:784635/1‑238 (MQ=255)
  cccATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGCTCGAGGTCGGTGACACTCCCGGCAAACAGACACACGGGAAAAatcataaacga                                                                                                                                                                                                                                >  2:210113/1‑233 (MQ=255)
  cccATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCCCGGCAATCAGCAAGGGGGGAAAAatcaaacccc                                                                                                                                                                                                                                 >  2:99630/1‑232 (MQ=255)
  cccATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGGCGATGACACCCGCGGCAAACAGCCAGGCGGATAAAAAAATCAcccc                                                                                                                                                                                                                                >  2:368596/1‑237 (MQ=255)
  cccATGCCGAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACACCGGGCAGCCATTGGTGATCGAGGTCGATGACACACGCGGCAATCAGCCAGGCGGATAATATCATCACCg                                                                                                                                                                                                                                 >  1:55589/1‑238 (MQ=255)
     aTGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGAAACTCGCGGCAATAAGCCAGGGGGTGAATATCCTCACCGccc                                                                                                                                                                                                                              >  1:357387/1‑237 (MQ=255)
       gCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGc                                                                                                                                                                                                                            >  1:154419/1‑238 (MQ=255)
              ttttttttACGCAGACTTTTCCGCGTATTCCAACGCAGAGAGAAAAAGGTGACAAACCCCCCCAGGCCGCCGGTGACGGCATCTTGTAGCGTGGGGGAACTCTGCTGCGCCCATGCCGCAATCAGTCCTGCCCACAATACGCCCTGAGTAAAACCATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATcc                                                                                                                                                                                                                    <  2:310881/234‑1 (MQ=255)
                  cttTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGCCGGGAAATATCATCACCGCCAGCCCCCATCCACtt                                                                                                                                                                                                                >  1:158801/1‑239 (MQ=255)
                  cttTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGCCGATGACAATCGCGGCAATCAGCCAGGCGGTAAATATCATCACCGCCAGCCCCCATCCACtt                                                                                                                                                                                                                >  1:31542/1‑239 (MQ=255)
                  cttTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGGGATCGAGGTCGATGACACTCGCGGCAATCAGCCAAGCGGTAAATATCATCACCCCCAGCCCCCATCCACtt                                                                                                                                                                                                                >  1:731934/1‑239 (MQ=255)
                        gCAGAATTATTCCGGCTATCCAACGCAGAGAGTAAAAAGGGATAAACCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGc                                                                                                                                                                                                          <  1:672072/239‑1 (MQ=255)
                            aaCTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCCCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCCCTTTCTGGccag                                                                                                                                                                                                       >  1:20598/1‑238 (MQ=255)
                                      gCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCg                                                                                                                                                                                             >  1:131660/1‑238 (MQ=255)
                                      gCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGAATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGGAAAAAAATCGTGCAGCCAGTGGGGACCAGGGGCGCGGCAATCCGCGGCATACAGCCAGGCGGTTAATATCATAACCGCCCCCCCCCCACCCATTACTGGgcaaaccaaaacca                                                                                                                                                                                             >  2:686632/1‑232 (MQ=255)
                                        cttcccaagcaagagagTAAAAGGTGAAAAACCCCACCAGGCCGCCGGTGACTGCATCTTGTAGCGTGAGGGGCCCCGGCGGCGCCCATGCCGCAATCAGTCCTGTCCAAAATCCGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGc                                                                                                                                                                                            <  1:587484/227‑1 (MQ=255)
                                           tCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGTCGGTGAATATAATCACCGCCAGCCCCCATCCACTTTCTGGCCCGCCCAGACTCGCCAGc                                                                                                                                                                                        >  1:753157/1‑239 (MQ=255)
                                           tCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGc                                                                                                                                                                                        >  1:245551/1‑239 (MQ=255)
                                             aacccaagagatAAAAAAGGTAAAACTCCCCCCAGGCCGCCGGTTACTGCATTGTGTCGCGGGGGGGCACTCGCCCGCCCCCATCCCGCAATCGTCCTTGCCCAAAACACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGCCGATGACCCTCGCGGCAATCAGCCAGGCGGAAAATACCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaa                                                                                                                                                                                      <  2:154419/233‑1 (MQ=255)
                                                       agTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGt                                                                                                                                                                            >  1:103153/1‑239 (MQ=255)
                                                         tAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGAATCGCCAGCAAAAACGAAGTGCt                                                                                                                                                                         >  1:304622/1‑239 (MQ=255)
                                                         aaaaaaGTGATAAATCCCCCCGGGCACCGGGTAACTGCATTTTGTAGCGGGGGCGGCCTCTGCGGCGCCCAGCCCCCAATCATTCCTGTCCCAAATCCCCCCTGAGTAAAAACATGGGGCACCCATGGGTGATCGAGGGCGATGACACCCGCGGCACTCACCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTg                                                                                                                                                                           <  2:460983/237‑1 (MQ=255)
                                                          aaaaaGTGATAAACCCCACCGGGACGCCGGTGACCGCATCTTGTAGGGTGAGCGGACTTTGCTGCGCCCATGCCGCAATCATTCCTGTCCACAATCCGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGACCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCt                                                                                                                                                                         <  1:368596/239‑1 (MQ=255)
                                                          aaaaaGTGATAAACCCCACCAGGACCCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCGGCTGCGCCAATGCCGCAAGCAGTCCGGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCACCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTg                                                                                                                                                                           <  1:99630/237‑1 (MQ=255)
                                                                 gATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCCCCGCCAGCCCCCATCCACTTTCTGGCCAGCCCAGACTCGCCAGCAAAAAAGCAAGTGCGGCCaaaa                                                                                                                                                                  >  1:600686/1‑237 (MQ=255)
                                                                 gATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATa                                                                                                                                                                  >  1:507378/1‑239 (MQ=255)
                                                                 gATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGGGCTGTCAATa                                                                                                                                                                  >  1:467535/1‑239 (MQ=255)
                                                                 gATAAATCCCACCAGGACGCCGGTGACCGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCCCTTTCTGGCCAGACCAGACTCGCCCGCCAAAAAGCAAGGGCTGTCAATa                                                                                                                                                                  >  1:290102/1‑239 (MQ=255)
                                                                       tCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCATCAAAAAAGCAAGCGCTGTCAATAACTcc                                                                                                                                                             >  1:496394/1‑238 (MQ=255)
                                                                                gcccccGGTGAGTTCATTTTTGAGCGGGAGCGGCCTCTGCTGCCCCCAGCCCGAAACGGGCCCTGTCCACAATACGCCTGGAGTAAAAACACCGGGCAGCCATTGGAGATCGAGGTCGATGAACCCGGGGACAATCAGCAAGGCGGAAAATACACTCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGAt                                                                                                                                                   <  2:73179/235‑1 (MQ=255)
                                                                                  cGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCCGCAAAAAAGCAAGTGCTTTCAATAACTCCACCACCGGCTa                                                                                                                                                  >  1:638010/1‑238 (MQ=255)
                                                                                     cGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCCCTTTCTGGCCAGACCAGACCCCCCAGCAAAAAACAAAAGGCTGGCATAAACTCCACCAGCggtaaca                                                                                                                                                >  2:804493/1‑232 (MQ=255)
                                                                                     cGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCAAGGCGGATAATATCATCACCGCCAGCCCCCATCCACGTTCTGGCCAGACCAGACTCGCCCGCAAAAAAGCAAATGCTGTCAATAACCACACCAGCGGGTAACgg                                                                                                                                              >  2:380637/1‑238 (MQ=255)
                                                                                                        cGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGCCACTCGCGGCCATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCACAAAAGCAAGTGCTGTCAATAACTCCACCCGCGGATAACGCTTGCTGATTTTCGCCTGAc                                                                                                                           >  1:33932/1‑239 (MQ=255)
                                                                                                          tGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGAGCTGTCAATaa                                                                                                                                                                 >  1:652318/1‑199 (MQ=255)
                                                                                                                            cccATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCCCTTTCTGGCCCGACCAGACTCCCCACCAAAAAACAAAGAGCTGTCAATAACTCCACCAGCGGAAAACGCTGGCGGATTTTCGCCTGACAGTCCAGGCacaccccattc                                                                                                       >  2:81317/1‑229 (MQ=255)
                                                                                                                               aTGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCACAAAAGCCAGTGCTGTCAATAACTCCACCAGCGGAAAAAGCTTGCTGATTTTCGCCTGCCCATCGGGGCGGCCCCCTTTGacc                                                                                                    >  2:739675/1‑237 (MQ=255)
                                                                                                                                 gCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAg                                                                                                                                                                                 >  1:307353/1‑160 (MQ=255)
                                                                                                                                 gCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTGGCTGATTTTCGCCTGACAGTCCCGGCAGCGCCCTTTGAGCa                                                                                                   >  1:528980/1‑238 (MQ=255)
                                                                                                                                            tCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCCCCAGCGGATAACGCTTGCTGATTTTCGCCTGCCATCCACGGCAGCGCCCTTTGAGGATCAGCCAGgaga                                                                                       >  1:715692/1‑239 (MQ=255)
                                                                                                                                                 gtcccaaaaaCGCCATAAGTAAAACCGTCGGGCCGCCTTTGGTGACCGAGGTCGATGACCCTCGCGGCAATCACCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAg                                                                                   <  1:229914/230‑1 (MQ=255)
                                                                                                                                                 gcccAAAATACCCCCTGAGTAAAAACAGCGGGCACCCATTGGTGACCGAGGTCGATGACACCGGCGGCAAGCGACCAGGCGGAAAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACACCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAg                                                                                   <  1:802082/237‑1 (MQ=255)
                                                                                                                                                            gCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGTGCACCAGCCAGGCGACAAGCGGAATATTg                                                                         >  1:457969/1‑238 (MQ=255)
                                                                                                                                                                 gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAGGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGCCAGTCGCGGCAGCGCCCTTTGAGCATAAGCCAGGAGAAAAGCGGAAAATTTTCAcct                                                                   >  2:787783/1‑237 (MQ=255)
                                                                                                                                                                 gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTTACga                                                                   >  1:558029/1‑238 (MQ=255)
                                                                                                                                                                 gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGt                                                                   >  1:727898/1‑239 (MQ=255)
                                                                                                                                                                 gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCAACAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGAAGCGCCCTTTGAGCATCAGCCAGGAGAACACCGGAAAATGGTCACga                                                                   >  1:235747/1‑238 (MQ=255)
                                                                                                                                                                 gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGGCGATGACACTTGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCACCGGATAACGCTTGCTGATTTTCGCCTGACAGTAGCGCCAACGCCCGTTTAGCTTAAGCCAGGAAGACAGAGGCATATTTACACga                                                                   >  1:573661/1‑238 (MQ=255)
                                                                                                                                                                         aCAGCGGGCCGCCATTGGTGAGCGAGGTCGAGGCCACTCGCGGACATCAGCAAGGCGGATAATATCATCCCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGAt                                                           <  1:695992/239‑1 (MQ=255)
                                                                                                                                                                            tCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGCCAGTCCCGGCAGCGCCCTTTGAGCCTCAGCCCGGGGAACAACGGAATATTGTAACGTATCCGGATgg                                                         >  1:407318/1‑238 (MQ=255)
                                                                                                                                                                            tCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCAc                                                                     >  1:778744/1‑226 (MQ=255)
                                                                                                                                                                                              tCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGCGAACAGCCGACTATTTTCACGAATCCGGATGGTCTTGTGACAATGTGGGCa                                      >  1:382748/1‑239 (MQ=255)
                                                                                                                                                                                              tCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGaa                                                                              >  1:714914/1‑199 (MQ=255)
                                                                                                                                                                                               cGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGGGAACAGCGGAAAATTGTCACGTATCCGGATGGTCGGCTGACCATGAGGGCAg                                     >  1:430405/1‑239 (MQ=255)
                                                                                                                                                                                                    tCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGAGGGTCTGCTGCAAACTTGGGCAGAGCGa                                >  1:367087/1‑239 (MQ=255)
                                                                                                                                                                                                     cGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGCTAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGGGAATCAGCCAGGAGAACAGCGGCAGATTGTAACGTATCCGGATGGCCTGCTGAAAATGAGGGCAGGGGGaa                               >  2:786772/1‑239 (MQ=255)
                                                                                                                                                                                                      gATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGAATCCGGATGGTCTGCTGACAATTTGGGCAGTGCGaag                              >  1:460883/1‑238 (MQ=255)
                                                                                                                                                                                                               cgcgGACATCAACCAGGCGGTAAATTTCACCACCGCCAGCCCCCATCCACTTTCGGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAgcg                     <  1:466773/239‑1 (MQ=255)
                                                                                                                                                                                                                gcgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACCCGGTAgcgc                    >  1:850892/1‑239 (MQ=255)
                                                                                                                                                                                                                      actcggCCGGGAGGAAAATATCACAACCGCCAGTCCCTATCCACTTTCTGGCCAGAACAGACTCGCCAGCAAAAAAGCAAGTTCTGGCAATAACTCCACCAGCCGATAACCCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGc               <  1:501282/234‑1 (MQ=255)
                                                                                                                                                                                                                      aaTCCGCCAGGGGGAAAATACCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGc               <  1:739675/239‑1 (MQ=255)
                                                                                                                                                                                                                      aaTCAGCCAGGCGAATAAAATCATCACCGCCAGCCCCCAGCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGc               <  1:898074/239‑1 (MQ=255)
                                                                                                                                                                                                                                    ataataTCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGCCAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTg  >  1:591009/1‑238 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
CGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTG  >  NC_000913/3113731‑3114193

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