Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,221 T→C 5.4% L132P (CTG→CCG) ‡ yhcG → DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,2210TC 5.4% 1613.2 / 53.6 529L132P (CTG→CCG) ‡yhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base T (237/242);  new base C (19/9);  total (272/257)
Fisher's exact test for biased strand distribution p-value = 7.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-02

Only 100 of 533 total aligned reads displayed.

CTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTT  >  NC_000913/3368072‑3368370
                                                                                                                                                     |                                                                                                                                                     
 tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTg                                                                                                                                                      >  2:3784652/1‑150 (MQ=255)
    gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                                                   <  1:3317440/150‑1 (MQ=255)
      aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                                 >  2:9132580/1‑150 (MQ=255)
      aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                                 >  1:3779377/1‑150 (MQ=255)
      aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                                 >  2:11223378/1‑150 (MQ=255)
       aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGcc                                                                                                                                                <  2:7744529/150‑1 (MQ=255)
       aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGcc                                                                                                                                                >  2:12221796/1‑150 (MQ=255)
        aaaaCTTACGCCACATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCa                                                                                                                                               <  2:6515882/150‑1 (MQ=255)
         aaaCTTACGCCAAATGAGGATTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCaa                                                                                                                                              <  2:7863016/150‑1 (MQ=255)
         aaaCCTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCaa                                                                                                                                              <  2:5203687/150‑1 (MQ=255)
          aaCTCACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCaaa                                                                                                                                             <  2:10311296/150‑1 (MQ=255)
             ttACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCt                                                                                                                                          <  1:10494401/150‑1 (MQ=255)
             ttACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCt                                                                                                                                          >  2:8441037/1‑150 (MQ=255)
               aCGCCAAATGCGGCATTTGTACCTCTTTTTTCACCATGCTGAAATTCACCAGACAAGGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCTGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCttt                                                                                                                                        <  2:4442282/150‑1 (MQ=255)
                cGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTc                                                                                                                                       <  2:2434648/150‑1 (MQ=255)
                 gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                                                      >  1:4053384/1‑150 (MQ=255)
                 gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                                                      >  2:1771909/1‑150 (MQ=255)
                  ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCGGTCTGCCAAATCTTTccc                                                                                                                                     <  2:10813576/150‑1 (MQ=255)
                   cAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCGGTATGCCAAATCTTTaccc                                                                                                                                    >  1:11392557/1‑150 (MQ=255)
                       tGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAAGTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGc                                                                                                                                <  2:12095923/150‑1 (MQ=255)
                       tGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAAATCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGACGCTGGCTGCCAAATATTTCCCCTTGc                                                                                                                                >  1:5594856/1‑150 (MQ=255)
                       tGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAAATCACCAGACAATGTCTGGCCAATTAACACCATTGGGAATTCCGCCGACACCGTCAGCAGACTTCCCAACGGCAAACATATGGCACACGCTGTCTGCCAAAACGTTCCCCATGc                                                                                                                                >  1:10215793/1‑150 (MQ=255)
                            cTTTTTAACCTCTTTTTTCAACCTGTTGAATTTCACCAGACAAGGTCTGGCTAATTAACACATTTGGTAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAgg                                                                                                                           <  1:5464615/150‑1 (MQ=255)
                                         tttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGTAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTcc                                                                                                              <  1:1505507/150‑1 (MQ=255)
                                            ttCAACATGTTGAAATTCACCAGACAATGTCGGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGtt                                                                                                           <  1:6687659/150‑1 (MQ=255)
                                                 cATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACtt                                                                                                      >  1:6105844/1‑150 (MQ=255)
                                                   tGTTGAAAGTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTc                                                                                                    <  2:8055510/150‑1 (MQ=255)
                                                    gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                                                   <  1:4936373/150‑1 (MQ=255)
                                                     ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                                                  >  1:10645729/1‑150 (MQ=255)
                                                     ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                                                  <  1:10891586/150‑1 (MQ=255)
                                                     ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAAGTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                                                  <  1:8752958/150‑1 (MQ=255)
                                                     ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCAGTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                                                  <  1:7777476/150‑1 (MQ=255)
                                                      tGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTAATTTCCGt                                                                                                 >  1:11858376/1‑150 (MQ=255)
                                                        aaaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTCCGTCCGTTTACTTTCCGtaa                                                                                               >  2:10502170/1‑148 (MQ=255)
                                                          agtCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaa                                                                                             <  1:2998836/148‑1 (MQ=255)
                                                          aTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaa                                                                                             >  1:521933/1‑150 (MQ=255)
                                                          aTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaa                                                                                             >  2:10976963/1‑150 (MQ=255)
                                                           ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaa                                                                                            >  1:11014390/1‑150 (MQ=255)
                                                           ttCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCCTCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCCGGTCCCCCTACGTCCGTTTACTGTCCGTTaaaa                                                                                            >  2:3424144/1‑150 (MQ=255)
                                                              aCCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACg                                                                                         <  2:4704663/150‑1 (MQ=255)
                                                               ccAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGc                                                                                        <  2:566245/150‑1 (MQ=255)
                                                                cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCAGCAGACACCGTCTGAAGAATTACCATCGGCAAAGACTAGGCAGACGGTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCAACCTCCGTTTACTTTCCGTTAAAAAAGCt                                                                                       >  1:7435565/1‑150 (MQ=255)
                                                                   aCAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGAc                                                                                    <  2:1433953/150‑1 (MQ=255)
                                                                        gTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCg                                                                               >  2:3127543/1‑150 (MQ=255)
                                                                           tGGCGAATTAACCCCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAg                                                                            <  1:8219388/150‑1 (MQ=255)
                                                                            ggCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGGTAAAAACGCTGACGCCcgccgc                                                                           >  2:8199540/1‑150 (MQ=255)
                                                                            ggCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCAGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGGCCGGTTAATTTCCGTTAAAAACGCTGACGCCCGCAGc                                                                           >  1:8928110/1‑150 (MQ=255)
                                                                              cGAATTAACACCATTGGGAATTCCGCAGACACCGTCGGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTGGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCtt                                                                         <  1:10053101/150‑1 (MQ=255)
                                                                               gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                        <  2:7103214/150‑1 (MQ=255)
                                                                               gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                        <  2:8792509/150‑1 (MQ=255)
                                                                               gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                        >  2:6219466/1‑150 (MQ=255)
                                                                                 aTTAACACCATTGGGATGTCCGCAGACACCGTCTGCAGAATGCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTTCGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                      <  1:9334407/150‑1 (MQ=255)
                                                                                      cacCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaa                                                                 <  2:6668261/150‑1 (MQ=255)
                                                                                      cacCAGTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaa                                                                 <  1:7991759/150‑1 (MQ=255)
                                                                                        ccATTGGTAATTCCGCACACACCGTCGGACTAATTCCCAGCGGCATATATTTGGCAGACGCTGTCTGCCAAATATTTCCCCTTGCCCTGGTCCACCTACTTCCGTTTACTTTCCTTTAAAAACGCTGACGCCCGCAGCTTTTATTaaaaa                                                               <  1:4253906/150‑1 (MQ=255)
                                                                                        ccATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                               <  1:1867450/150‑1 (MQ=255)
                                                                                        ccATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                               >  2:6888708/1‑150 (MQ=255)
                                                                                         cATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCCAAAATTTCCCCATGCCCCGGTCCACCTACGTCCGGTTACTTTACGGTAAAAAAGCTGACGCCCGCAGCTTTTATGtaaaat                                                              >  1:1169467/1‑149 (MQ=255)
                                                                                         cATTGGGAATGCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAag                                                              <  1:7898034/150‑1 (MQ=255)
                                                                                          aTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAATATTTTGCAGAAGCTGTCTGCCAAATATTTACCCTTGCCCAGGTCCACCAACGTCCGTTTAATTTACCTTAAAAAACCTGACGCCCGCAGCGTTTATGaaaaata                                                             >  1:4732861/1‑148 (MQ=255)
                                                                                          aTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGTCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                             <  1:6204424/150‑1 (MQ=255)
                                                                                           ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAagag                                                            <  1:9105120/150‑1 (MQ=255)
                                                                                           ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTAATTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGCAAAagag                                                            >  2:4241120/1‑150 (MQ=255)
                                                                                              ggAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACg                                                         <  1:6556277/150‑1 (MQ=255)
                                                                                                  ttCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTTGCAGACGCTGTCTGCCAAAACATTACCCATGCCCACGTCCACCAACCTCCGTTTACATTCCGTTAAAAAAGCATACGCCCCCCGCATTTATGAAAAATAGACGCTcc                                                     >  1:5965772/1‑150 (MQ=255)
                                                                                                     cGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCt                                                  >  2:4982163/1‑150 (MQ=255)
                                                                                                       cAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCAACGTCCGTTTACTTTCCGGTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt                                                >  1:6944805/1‑150 (MQ=255)
                                                                                                            acCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTg                                           <  1:1663686/150‑1 (MQ=255)
                                                                                                              cGTCTGCAGACTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt                                         <  2:10634432/150‑1 (MQ=255)
                                                                                                              cGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCGGTCTGCCAAATCTTTCCCCTCGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt                                         <  1:9963650/150‑1 (MQ=255)
                                                                                                               gTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTc                                        <  1:7610468/150‑1 (MQ=255)
                                                                                                                  tGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt                                     >  2:11430789/1‑150 (MQ=255)
                                                                                                                   gCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGtt                                    >  2:7624223/1‑150 (MQ=255)
                                                                                                                    cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                                   >  2:4005606/1‑150 (MQ=255)
                                                                                                                     aGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAGGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCg                                  >  2:9813403/1‑150 (MQ=255)
                                                                                                                       aaTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGcc                                <  1:11129942/150‑1 (MQ=255)
                                                                                                                       aaTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGcc                                <  1:5035384/150‑1 (MQ=255)
                                                                                                                            ccATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTa                           >  2:2690635/1‑150 (MQ=255)
                                                                                                                              aTCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                       >  2:1181071/1‑88 (MQ=255)
                                                                                                                              aTCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTaga                         >  2:643613/1‑150 (MQ=255)
                                                                                                                               tCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTagag                        >  2:8568583/1‑150 (MQ=255)
                                                                                                                                cGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGc                       <  1:445912/150‑1 (MQ=255)
                                                                                                                                  gCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCgg                     <  1:720195/150‑1 (MQ=255)
                                                                                                                                   cAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGc                    <  2:11659787/150‑1 (MQ=255)
                                                                                                                                    aaaGATTTGGCAGACGCTGTCTGCCAAAAAATTACCCTTGCCCAGGTCCACCTACGTCCCTTTAAATTACCTTAAAAAAGCAGAACCCCCCAGCATTTATGAAAAAGAGAAGCTCCGCAGTGGCTGGTCTGTTCGCCAGCAAGGGCGGCa                   >  1:3668806/1‑150 (MQ=255)
                                                                                                                                      cgatccggCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGGCCACCTACGTCCGTTTACTTTCCGTTAACCACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCGCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCaaa                 <  2:937984/144‑1 (MQ=255)
                                                                                                                                        aTTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAtt               <  1:5778404/150‑1 (MQ=255)
                                                                                                                                        aTTTGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAAGCTGACGCCCGCAGCTTTTATGTAAAAGTGACGCTCCGCTGTGGCTGGTCTGTTAGCCACCGAGAGCGGCAAAtt               >  1:2142942/1‑150 (MQ=255)
                                                                                                                                         tttGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCATGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg              <  2:5573713/150‑1 (MQ=255)
                                                                                                                                         tttGGAAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTTCTTTACGTTAAAATCGCTGACGCCGGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCATATTg              <  2:6008/150‑1 (MQ=255)
                                                                                                                                         gttGGCAGACGCTGTCTGCCAAATCTTTCCCCTCGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg              <  2:2879346/149‑1 (MQ=255)
                                                                                                                                          ttGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc             <  1:2770089/150‑1 (MQ=255)
                                                                                                                                          ttGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc             <  2:2138905/150‑1 (MQ=255)
                                                                                                                                          ttGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc             >  2:3626571/1‑150 (MQ=255)
                                                                                                                                           tGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCg            <  1:8516107/150‑1 (MQ=255)
                                                                                                                                             gCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAc          <  1:2377210/150‑1 (MQ=255)
                                                                                                                                             gCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAc          >  1:9637773/1‑150 (MQ=255)
                                                                                                                                               aGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAccc        >  2:9520384/1‑150 (MQ=255)
                                                                                                                                                  cGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGt     <  1:12221796/150‑1 (MQ=255)
                                                                                                                                                     tGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGtttt  <  2:2330754/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTT  >  NC_000913/3368072‑3368370

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