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→TTT) pppAbifunctional 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

CGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTG  >  NC_000913/3113731‑3114193
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
cGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGGTATCGAGGTCGATGACACTCGCGG‑CACTCACCCAGCCGGGTACTATCATCAc                                                                                                                                                                                                                                   >  1:784635/1‑238 (MQ=255)
  cccATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGCTCGAGGTCGGTGACACTCCCGG‑CAAACAGACACACGGGAAAAatcataaacga                                                                                                                                                                                                                                >  2:210113/1‑233 (MQ=255)
  cccATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCCCGG‑CAATCAGCAAGGGGGGAAAAatcaaacccc                                                                                                                                                                                                                                 >  2:99630/1‑232 (MQ=255)
  cccATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGGCGATGACACCCGCGG‑CAAACAGCCAGGCGGATAAAAAAATCAcccc                                                                                                                                                                                                                                >  2:368596/1‑237 (MQ=255)
  cccATGCCGAATGCTTCTTTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACACCGGGCAGCCATTGGTGATCGAGGTCGATGACACACGCGG‑CAATCAGCCAGGCGGATAATATCATCACCg                                                                                                                                                                                                                                 >  1:55589/1‑238 (MQ=255)
     aTGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGAAACTCGCGG‑CAATAAGCCAGGGGGTGAATATCCTCACCGccc                                                                                                                                                                                                                              >  1:357387/1‑237 (MQ=255)
       gCCTAATGCTTCTTTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGc                                                                                                                                                                                                                            >  1:154419/1‑238 (MQ=255)
              ttttttttACGCAG‑ACTTTTCCGCGTATTCCAACGCAGAGAGAAAAAGGTGACAAACCCCCCCAGGCCGCCGGTGACGGCATCTTGTAGCGTGGGGGAACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGCCCACAATACGCCCTGAGTAAAACCATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATcc                                                                                                                                                                                                                    <  2:310881/234‑1 (MQ=255)
                  cttTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGCCGGGAAATATCATCACCGCCAGCCCCCATCCACtt                                                                                                                                                                                                                >  1:158801/1‑239 (MQ=255)
                  cttTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGCCGATGACAATCGCGG‑CAATCAGCCAGGCGGTAAATATCATCACCGCCAGCCCCCATCCACtt                                                                                                                                                                                                                >  1:31542/1‑239 (MQ=255)
                  cttTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGGGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAAGCGGTAAATATCATCACCCCCAGCCCCCATCCACtt                                                                                                                                                                                                                >  1:731934/1‑239 (MQ=255)
                        gCAGAATTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGGGATAAACCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGc                                                                                                                                                                                                          <  1:672072/239‑1 (MQ=255)
                            aaCTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCCCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCCCTTTCTGGccag                                                                                                                                                                                                       >  1:20598/1‑238 (MQ=255)
                                      gCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCg                                                                                                                                                                                             >  1:131660/1‑238 (MQ=255)
                                      gCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGAATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGGAAAAAAATCGTGCAGCCAGTGGGGACCAGGGGCGCGGCAATCCGCGG‑CATACAGCCAGGCGGTTAATATCATAACCGCCCCCCCCCCACCCATTACTGGgcaaaccaaaacca                                                                                                                                                                                             >  2:686632/1‑232 (MQ=255)
                                        cttcccaagcaagagagTAAAAGGTGAAAAACCCCACCAGGCCGCCGGTGACTGCATCTTGTAGCGTGAGGGGCCCCGGCGGCGCCCATGCCGCAATCAGTCC‑TGTCCAAAATCCGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGc                                                                                                                                                                                            <  1:587484/227‑1 (MQ=255)
                                           tCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGTCGGTGAATATAATCACCGCCAGCCCCCATCCACTTTCTGGCCCGCCCAGACTCGCCAGc                                                                                                                                                                                        >  1:753157/1‑239 (MQ=255)
                                           tCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGc                                                                                                                                                                                        >  1:245551/1‑239 (MQ=255)
                                             aacccaagagatAAAAAAGGTAAAACTCCCCCCAGGCCGCCGGTTACTGCATTGTGTCGCGGGGGGGCACTCGCCCGCCCCCATCCCGCAATC‑GTCCTTGCCCAAAACACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGCCGATGACCCTCGCGG‑CAATCAGCCAGGCGGAAAATACCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaa                                                                                                                                                                                      <  2:154419/233‑1 (MQ=255)
                                                       agTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGt                                                                                                                                                                            >  1:103153/1‑239 (MQ=255)
                                                         tAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGAATCGCCAGCAAAAACG‑AAGTGCt                                                                                                                                                                         >  1:304622/1‑239 (MQ=255)
                                                         aaaaaaGTGATAAATCCCCCCGGGCACCGGGTAACTGCATTTTGTAGCGGGGGCGGCCTCTGCGGCGCCCAGCCCCCAATCATTCC‑TGTCCCAAATCCCCCCTGAGTAAAAACATGGGGCACCCATGGGTGATCGAGGGCGATGACACCCGCGG‑CACTCACCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTg                                                                                                                                                                           <  2:460983/237‑1 (MQ=255)
                                                          aaaaaGTGATAAACCCCACCGGGACGCCGGTGACCGCATCTTGTAGGGTGAGCGGACTTTGCTGCGCCCATGCCGCAATCATTCC‑TGTCCACAATCCGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGACCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCt                                                                                                                                                                         <  1:368596/239‑1 (MQ=255)
                                                          aaaaaGTGATAAACCCCACCAGGACCCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCGGCTGCGCCAATGCCGCAAGCAGTCC‑GGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCACCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTg                                                                                                                                                                           <  1:99630/237‑1 (MQ=255)
                                                                 gATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCCCCGCCAGCCCCCATCCACTTTCTGGCCAGCCCAGACTCGCCAGCAAAAAAGCAAGTGCGGCCaaaa                                                                                                                                                                  >  1:600686/1‑237 (MQ=255)
                                                                 gATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATa                                                                                                                                                                  >  1:507378/1‑239 (MQ=255)
                                                                 gATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGGGCTGTCAATa                                                                                                                                                                  >  1:467535/1‑239 (MQ=255)
                                                                 gATAAATCCCACCAGGACGCCGGTGACCGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCCCTTTCTGGCCAGACCAGACTCGCCCGCCAAAAAGCAAGGGCTGTCAATa                                                                                                                                                                  >  1:290102/1‑239 (MQ=255)
                                                                       tCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCATCAAAAAAGCAAGCGCTGTCAATAACTcc                                                                                                                                                             >  1:496394/1‑238 (MQ=255)
                                                                                gcccccGGTGAGTTCATTTTTGAGCGGGAGCGGCCTCTGCTGCCCCCAGCCCGAAACGGGCCC‑TGTCCACAATACGCCTGGAGTAAAAACACCGGGCAGCCATTGGAGATCGAGGTCGATGA‑ACCCGGGGACAATCAGCAAGGCGGAAAATACACTCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGAt                                                                                                                                                   <  2:73179/235‑1 (MQ=255)
                                                                                  cGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCCGCAAAAAAGCAAGTGCTTTCAATAACTCCACCACCGGCTa                                                                                                                                                  >  1:638010/1‑238 (MQ=255)
                                                                                     cGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCCCTTTCTGGCCAGACCAGACCCCCCAGCAAAAAACAAAAGGCTGGCATAAACTCCACCAGCggtaaca                                                                                                                                                >  2:804493/1‑232 (MQ=255)
                                                                                     cGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCAAGGCGGATAATATCATCACCGCCAGCCCCCATCCACGTTCTGGCCAGACCAGACTCGCCCGCAAAAAAGCAAATGCTGTCAATAACCACACCAGCGGGTAACgg                                                                                                                                              >  2:380637/1‑238 (MQ=255)
                                                                                                        cGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGCCACTCGCGG‑CCATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCACAAAAGCAAGTGCTGTCAATAACTCCACCCGCGGATAACGCTTGCTGATTTTCGCCTGAc                                                                                                                           >  1:33932/1‑239 (MQ=255)
                                                                                                          tGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGAGCTGTCAATaa                                                                                                                                                                 >  1:652318/1‑199 (MQ=255)
                                                                                                                            cccATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCCCTTTCTGGCCCGACCAGACTCCCCACCAAAAAACAAAGAGCTGTCAATAACTCCACCAGCGGAAAACGCTGGCGGATTTTCGCCTGACAGTCCAGGCacaccccattc                                                                                                       >  2:81317/1‑229 (MQ=255)
                                                                                                                               aTGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCACAAAAGCCAGTGCTGTCAATAACTCCACCAGCGGAAAAAGCTTGCTGATTTTCGCCTGCCCATCGGGGCGGCCCCCTTTGacc                                                                                                    >  2:739675/1‑237 (MQ=255)
                                                                                                                                 gCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAg                                                                                                                                                                                 >  1:307353/1‑160 (MQ=255)
                                                                                                                                 gCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTGGCTGATTTTCGCCTGACAGTCCCGGCAGCGCCCTTTGAGCa                                                                                                   >  1:528980/1‑238 (MQ=255)
                                                                                                                                            tCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCCCCAGCGGATAACGCTTGCTGATTTTCGCCTGCCATCCACGGCAGCGCCCTTTGAGGATCAGCCAGgaga                                                                                       >  1:715692/1‑239 (MQ=255)
                                                                                                                                                 gtcccaaaaaCGCCATAAGTAAAACCGTCGGGCCGCCTTTGGTGACCGAGGTCGATGACCCTCGCGG‑CAATCACCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAg                                                                                   <  1:229914/230‑1 (MQ=255)
                                                                                                                                                 gcccAAAATACCCCCTGAGTAAAAACAGCGGGCACCCATTGGTGACCGAGGTCGATGACACCGGCGG‑CAAGCGACCAGGCGGAAAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACACCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAg                                                                                   <  1:802082/237‑1 (MQ=255)
                                                                                                                                                            gCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGTGCACCAGCCAGGCGACAAGCGGAATATTg                                                                         >  1:457969/1‑238 (MQ=255)
                                                                                                                                                                 gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAGGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGCCAGTCGCGGCAGCGCCCTTTGAGCATAAGCCAGGAGAAAAGCGGAAAATTTTCAcct                                                                   >  2:787783/1‑237 (MQ=255)
                                                                                                                                                                 gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTTACga                                                                   >  1:558029/1‑238 (MQ=255)
                                                                                                                                                                 gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGt                                                                   >  1:727898/1‑239 (MQ=255)
                                                                                                                                                                 gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCAACAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGAAGCGCCCTTTGAGCATCAGCCAGGAGAACACCGGAAAATGGTCACga                                                                   >  1:235747/1‑238 (MQ=255)
                                                                                                                                                                 gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGGCGATGACACTTGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCACCGGATAACGCTTGCTGATTTTCGCCTGACAGTAGCGCCAACGCCCGTTTAGCTTAAGCCAGGAAGACAGAGGCATATTTACACga                                                                   >  1:573661/1‑238 (MQ=255)
                                                                                                                                                                         aCAGCGGGCCGCCATTGGTGAGCGAGGTCGAGGCCACTCGCGG‑ACATCAGCAAGGCGGATAATATCATCCCCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGAt                                                           <  1:695992/239‑1 (MQ=255)
                                                                                                                                                                            tCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGCCAGTCCCGGCAGCGCCCTTTGAGCCTCAGCCCGGGGAACAACGGAATATTGTAACGTATCCGGATgg                                                         >  1:407318/1‑238 (MQ=255)
                                                                                                                                                                            tCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCAc                                                                     >  1:778744/1‑226 (MQ=255)
                                                                                                                                                                                              tCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGCGAACAGCCGACTATTTTCACGAATCCGGATGGTCTTGTGACAATGTGGGCa                                      >  1:382748/1‑239 (MQ=255)
                                                                                                                                                                                              tCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGaa                                                                              >  1:714914/1‑199 (MQ=255)
                                                                                                                                                                                               cGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGGGAACAGCGGAAAATTGTCACGTATCCGGATGGTCGGCTGACCATGAGGGCAg                                     >  1:430405/1‑239 (MQ=255)
                                                                                                                                                                                                    tCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGAGGGTCTGCTGCAAACTTGGGCAGAGCGa                                >  1:367087/1‑239 (MQ=255)
                                                                                                                                                                                                     cGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGCTAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGGGAATCAGCCAGGAGAACAGCGGCAGATTGTAACGTATCCGGATGGCCTGCTGAAAATGAGGGCAGGGGGaa                               >  2:786772/1‑239 (MQ=255)
                                                                                                                                                                                                      gATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGAATCCGGATGGTCTGCTGACAATTTGGGCAGTGCGaag                              >  1:460883/1‑238 (MQ=255)
                                                                                                                                                                                                               cgcgG‑ACATCAACCAGGCGGTAAATTTCACCACCGCCAGCCCCCATCCACTTTCGGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAgcg                     <  1:466773/239‑1 (MQ=255)
                                                                                                                                                                                                                gcgG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGCGGGCAGTGCGAACCCGGTAgcgc                    >  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)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
CGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTA‑TCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCC‑TGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGG‑CAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGACAGTCGCGGCAGCGCCCTTTGAGCATCAGCCAGGAGAACAGCGGAATATTGTCACGTATCCGGATGGTCTGCTGACAATGTGGGCAGTGCGAACGCGGTAGCGCAAGGCTTATTTTTGACTG  >  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: ‑