Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,219 G→A 6.5% T131T (ACG→ACA 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,2190GA 6.5% 2339.3 / 94.3 727T131T (ACG→ACAyhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base G (379/299);  new base A (33/14);  total (414/313)
Fisher's exact test for biased strand distribution p-value = 6.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 737 total aligned reads displayed.

TTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTT  >  NC_000913/3368070‑3368368
                                                                                                                                                     |                                                                                                                                                     
  ctctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCt                                                                                                                                                     >  2:4975641/1‑150 (MQ=255)
    ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt                                                                                                                                                   >  1:7988374/1‑150 (MQ=255)
     tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGCATTCCCATCGGCAAAGATTTGGCAGACGCTGTc                                                                                                                                                  >  2:748688/1‑150 (MQ=255)
      gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                                                 <  1:4442387/150‑1 (MQ=255)
       cAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTg                                                                                                                                                <  1:3183562/150‑1 (MQ=255)
        aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                               >  2:12089302/1‑150 (MQ=255)
        aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                               >  1:9710731/1‑150 (MQ=255)
        aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                               >  2:8066809/1‑150 (MQ=255)
          aaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCa                                                                                                                                             >  2:15115311/1‑150 (MQ=255)
          aaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCa                                                                                                                                             >  1:13676544/1‑150 (MQ=255)
           aaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGccca                                                                                                                                            >  2:2031379/1‑148 (MQ=255)
            aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCaaa                                                                                                                                           <  2:8343784/150‑1 (MQ=255)
            aaCTCACGCCACATGAGGCTTTTTTACCGCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCaaa                                                                                                                                           <  2:9742355/150‑1 (MQ=255)
             aCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAAt                                                                                                                                          <  2:376280/150‑1 (MQ=255)
              cTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATc                                                                                                                                         <  1:4926431/150‑1 (MQ=255)
              cTTACGCCAAAGGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATc                                                                                                                                         <  1:10381475/150‑1 (MQ=255)
                   gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                                                    >  1:15453326/1‑150 (MQ=255)
                     cAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcccc                                                                                                                                  >  2:1566015/1‑150 (MQ=255)
                        aTGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTg                                                                                                                               >  1:1120743/1‑150 (MQ=255)
                                 ttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTcc                                                                                                                      <  2:13864007/150‑1 (MQ=255)
                                 ttttACCTCTTTTTTCAAAATGTTTAAATTCACCAGACATTGTCTGGCGAATTAACACCTTTGGTAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTcc                                                                                                                      <  1:5836066/150‑1 (MQ=255)
                                    tACCTCTTTTTTCAACATGTTGAAATTCACCAGACATTGTCTGGCGAATTAACACCTTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTGGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAcc                                                                                                                   <  1:7106534/150‑1 (MQ=255)
                                      ccTCTTTTTTCAACATGTTGAAATACACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCGGCAGTATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTa                                                                                                                 <  1:6185588/150‑1 (MQ=255)
                                       ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCGTGCCCAGGGCCACCTAc                                                                                                                >  1:6901286/1‑150 (MQ=255)
                                          ttttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTc                                                                                                             >  2:6280029/1‑150 (MQ=255)
                                            ttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCg                                                                                                           <  2:1441505/150‑1 (MQ=255)
                                                  aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACt                                                                                                     >  2:2841926/1‑150 (MQ=255)
                                                    aagTTGAAATTCACCAGAAAATGTATGGCGAATTAACAACATAGGGAATTCCGCAGACACCGTATGCAGAATTACCAACGGCAAAGAATTGGCAGAAGCTGTCTGACAAATCTTTCCCCTTGCCAAGGTCCACCTACGTCCGTTTAAttt                                                                                                   >  2:856472/3‑150 (MQ=255)
                                                      gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                                                 <  1:8990041/150‑1 (MQ=255)
                                                      gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATATTTCCCCTTGCCCAGGTCCACCTACGGCCGGTTACTTTcc                                                                                                 >  2:4682129/1‑150 (MQ=255)
                                                       ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACTTTCCg                                                                                                >  2:10557934/1‑150 (MQ=255)
                                                        tGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGt                                                                                               <  1:10610230/150‑1 (MQ=255)
                                                             ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCCAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTAATTTCCGTTaaaa                                                                                          >  1:8587759/1‑150 (MQ=255)
                                                             ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaa                                                                                          >  1:3369976/1‑150 (MQ=255)
                                                               cACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAc                                                                                        <  1:5177493/150‑1 (MQ=255)
                                                                  cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCCGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                     >  2:14720026/1‑150 (MQ=255)
                                                                  cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                     <  1:15041967/150‑1 (MQ=255)
                                                                  cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAAGTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                     <  1:940366/150‑1 (MQ=255)
                                                                   aGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCATGCCCAGGTCCCCCAACGTCCGTTTACATTCCGTTAAAAAAGCTg                                                                                    >  1:13307429/1‑150 (MQ=255)
                                                                        aTGTCTGGCGAATTAACACCATTGGCAATTCCGCACACACCCACTCCACAATTCCCATCGGCAAACATTTGCCAGACGCTGTCTGCAAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCCTTAAAAACGCTGACGcc                                                                               <  1:13559862/150‑1 (MQ=255)
                                                                          gTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGGCCACCTACGTCCGGTTACTTTCCGTTAAAAACGCTGACGCCCg                                                                             >  2:6486256/1‑150 (MQ=255)
                                                                             tGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCAACGTCCGGTTACATTCCGTTAAAAAACCAGACGCCCCCAg                                                                          >  1:14630249/1‑150 (MQ=255)
                                                                             tGGCGAATTAACACCATTGGGAATTCCGCAGACACAGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAg                                                                          <  2:8845971/150‑1 (MQ=255)
                                                                               gCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCt                                                                        >  2:12924358/1‑150 (MQ=255)
                                                                                 gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                      >  2:310507/1‑150 (MQ=255)
                                                                                 gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                      <  1:9597499/150‑1 (MQ=255)
                                                                                 gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttg                                                                      >  2:7215850/1‑149 (MQ=255)
                                                                                   aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                    <  2:13478006/150‑1 (MQ=255)
                                                                                       acacCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaa                                                                >  2:2261852/1‑150 (MQ=255)
                                                                                          cccttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                             <  1:12170150/147‑1 (MQ=255)
                                                                                          ccATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                             <  2:3989823/150‑1 (MQ=255)
                                                                                            gatGGGAATTCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCCGACGCTGTCTTCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACTTTCCGTTAAAAAAGCTGACGCCCGCAGCTTTTATGAAAAagc                                                           >  1:8709866/3‑149 (MQ=38)
                                                                                            gatGGGAATTCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAAGCTGACGCCCGCAGCTTTTATGAAAAagt                                                           >  1:12579251/3‑149 (MQ=39)
                                                                                             ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAagag                                                          >  2:1010074/1‑150 (MQ=255)
                                                                                              tGGGAATTCCGCAGACACCGTCTGCATAATTCCCATCGGCCAAGCTTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTAAGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                         <  1:10188210/150‑1 (MQ=255)
                                                                                              tGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                         <  1:1422674/150‑1 (MQ=255)
                                                                                              tGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAGGAAACAgaga                                                         >  2:14187458/1‑150 (MQ=255)
                                                                                               gggAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGGCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTGTTATGAAAAAGAGAc                                                        >  2:11486463/1‑150 (MQ=255)
                                                                                               gggAATTCCGCAGACACCGTCAGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACGTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAACAAGAGAc                                                        >  2:345228/1‑150 (MQ=255)
                                                                                                ggAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                     >  2:5204188/1‑120 (MQ=255)
                                                                                                 gAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGACGCGGGCTGCCAAAAATTTAACCATTCCCAGGTCCACCAACGTCCGGTTAAATTACGTTAAAAAAGCTGACGCCCCCAGCGTTTATGTAAAagagtagc                                                      >  1:4673708/1‑146 (MQ=255)
                                                                                                 gAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGc                                                      <  2:449983/150‑1 (MQ=255)
                                                                                                  aaTTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAATGATTTGGCAGACGCGGTCTGCCCAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCGGACGCCCGCAGCTTTTGTGGAAAAGAGACGcg                                                     >  1:8359927/1‑149 (MQ=255)
                                                                                                       cGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCAAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCt                                                <  1:7736343/150‑1 (MQ=255)
                                                                                                         cAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt                                              <  2:12604731/150‑1 (MQ=255)
                                                                                                         cAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt                                              <  1:4109965/150‑1 (MQ=255)
                                                                                                          aGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTATCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgtg                                             <  1:2031379/150‑1 (MQ=255)
                                                                                                              cccGTCGGCAGAAGGCCCAGCGGCACCGATTGGGCAGCCGCTGTCTGCCACATCTTTCCCCTTGCCCAGGTCCCCCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTg                                         <  2:15008819/149‑1 (MQ=255)
                                                                                                              acCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTg                                         >  2:10270484/1‑150 (MQ=255)
                                                                                                                cGTCTGCAGAATTCCCATCGGCAAAGATTGGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt                                       <  1:14058532/150‑1 (MQ=255)
                                                                                                                 gTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTc                                      <  1:3569959/150‑1 (MQ=255)
                                                                                                                  tCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCt                                     <  1:2598939/150‑1 (MQ=255)
                                                                                                                      cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                                 <  1:2890134/150‑1 (MQ=255)
                                                                                                                         aaTTCCCATCGGCCAAGCTTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGcc                              <  1:11811076/150‑1 (MQ=255)
                                                                                                                         aaTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGtt                                  >  2:1044114/1‑146 (MQ=255)
                                                                                                                             cccATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCt                          <  2:11718951/150‑1 (MQ=255)
                                                                                                                               cATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCg                                <  1:1691236/142‑1 (MQ=255)
                                                                                                                               cATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTag                        >  2:9458313/1‑150 (MQ=255)
                                                                                                                                aTCGGCAAAGATTTGGCAGCCGCTGTCTGCCCAATCTTTCCCCGTGCCCAGGGCCCCCTACGGCCGTTTACTTTCCGTTAAAAACGCAGACGCCCGCAGCTTTTCTGAAAAAGCGACGCTCCGCGGGGGCTGGTCTGGTCGCCAGCTaga                       >  2:13269376/1‑150 (MQ=255)
                                                                                                                                 tCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTagag                      <  1:7609687/150‑1 (MQ=255)
                                                                                                                                    gCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCgg                   >  2:9089379/1‑150 (MQ=255)
                                                                                                                                     cAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                             >  2:6917350/1‑107 (MQ=255)
                                                                                                                                      aaaGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGTAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCa                 >  1:10946402/1‑150 (MQ=255)
                                                                                                                                          aTTTGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAtt             >  1:11621280/1‑150 (MQ=255)
                                                                                                                                           tttGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg            <  1:9826784/150‑1 (MQ=255)
                                                                                                                                           tttGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg            >  2:5652298/1‑150 (MQ=255)
                                                                                                                                           tttGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg            >  2:586020/1‑150 (MQ=255)
                                                                                                                                           tttGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCa                             >  1:10038693/1‑133 (MQ=255)
                                                                                                                                           tttGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg            >  1:3089461/1‑150 (MQ=255)
                                                                                                                                            ttGGCAGACAGCGTCTGCCAAATCTTTGCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc           >  2:5382516/1‑150 (MQ=255)
                                                                                                                                             tGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCg          <  2:7681912/150‑1 (MQ=255)
                                                                                                                                             tGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCg          >  2:14856563/1‑150 (MQ=255)
                                                                                                                                              ggCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCACATTGCGa         >  2:10624176/1‑150 (MQ=255)
                                                                                                                                              ggCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGa         <  2:4293156/150‑1 (MQ=255)
                                                                                                                                               gCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAc        >  1:5581992/1‑150 (MQ=255)
                                                                                                                                                 aGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                                                 >  1:12712867/1‑59 (MQ=255)
                                                                                                                                                 aGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAccc      >  1:13017014/1‑150 (MQ=255)
                                                                                                                                                  gACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCa     <  1:9973231/150‑1 (MQ=255)
                                                                                                                                                   aCGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGATTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAg    >  2:9248296/1‑150 (MQ=255)
                                                                                                                                                     gCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGtt  >  2:11218239/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTT  >  NC_000913/3368070‑3368368

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: