Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,221 T→C 5.5% 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.5% 1043.6 / 34.7 331L132P (CTG→CCG) ‡yhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base T (145/158);  new base C (11/7);  total (163/168)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-02

Only 100 of 335 total aligned reads displayed.

TCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTT  >  NC_000913/3368073‑3368370
                                                                                                                                                    |                                                                                                                                                     
tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATTTTGAAATTCACCAGACTATGTCTGGCGAATTAACACCATTGGTAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTg                                                                                                                                                      <  1:5827041/150‑1 (MQ=255)
 ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt                                                                                                                                                     >  2:3592951/1‑150 (MQ=255)
   tcAAAAAACTTACGCCAAATGAGTCTTTTTTACCTCTTTTTTCAACTTGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTACGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                                                   <  2:6039167/149‑1 (MQ=255)
   gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACCTGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                                                   <  2:1780108/150‑1 (MQ=255)
   gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACCTGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                                                   <  2:2021048/150‑1 (MQ=255)
   gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                                                   <  2:1564122/150‑1 (MQ=255)
   gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                                                   <  1:5123558/150‑1 (MQ=255)
     aaacaaCTTCCGCCACATGAGGCTTTTTTACCTCTTTTTTCAACTGGTTGAAATTCACCAGCCACTGTCGGGCGACTGAACCCCATTGGGACTTCCGCAGACACCGTCGGCAGAATTCCCATCGGCAAAGATTTGGCAGCCGCTGACTGc                                                                                                                                                 <  2:2952147/150‑1 (MQ=255)
     aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                                 <  2:6626790/150‑1 (MQ=255)
     aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                                 <  1:2426539/150‑1 (MQ=255)
      aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGcc                                                                                                                                                >  1:6268272/1‑150 (MQ=255)
         aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCaaa                                                                                                                                             <  1:4935070/150‑1 (MQ=255)
            ttACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCt                                                                                                                                          >  2:2231320/1‑150 (MQ=255)
            ttACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCt                                                                                                                                          >  1:273775/1‑150 (MQ=255)
               cGCCACATGAGGCTTTTTTACCTCTTTTTTCAACCTGTTGAAATTCACCAGACACTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTc                                                                                                                                       <  2:1859730/150‑1 (MQ=255)
                          gCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAg                                                                                                                            >  1:4484594/1‑150 (MQ=255)
                          gCTTTTTTACCTCTTTTTTAAACATGTTGAAATACACCAGAAATTGTGTGGCGATTTAACACCTTTGGGAATTCCGCTGCCACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAg                                                                                                                            <  1:2639622/150‑1 (MQ=255)
                             tttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCTAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTc                                                                                                                         <  1:6178902/150‑1 (MQ=255)
                                ttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                      <  2:2095399/150‑1 (MQ=255)
                                ttACCTCTTTACTCAACAAGTTGAAATTCACCAGACAATCTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAAATCCCATCGGCAAAGATTTGGCAGACGCTGGCTGCCAAATCTTTCCCCTAGCCCCGGTcccc                                                                                                                      >  2:3991787/1‑148 (MQ=255)
                                  aCCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCt                                                                                                                    >  2:4036292/1‑150 (MQ=255)
                                    ctctTTTTTCAACATGTTGAATTTCACCAGACATTGTCTGGCGAATTACCACCATTGGGAATTCCGCAGACACCGTCTGCAGAATGCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                  <  1:3340381/150‑1 (MQ=255)
                                    ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCACAGATTTGGCAGACGCTGTCTGCCAAAGCTTTCCCCTTGCCCGGGTCCACCTAc                                                                                                                  >  2:3442283/1‑150 (MQ=255)
                                    ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                  <  2:718590/150‑1 (MQ=255)
                                    ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                  <  2:4616775/150‑1 (MQ=255)
                                      ctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGt                                                                                                                >  1:1021786/1‑150 (MQ=255)
                                       ttttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTc                                                                                                               >  1:1983170/1‑150 (MQ=255)
                                               aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACt                                                                                                       >  2:5019165/1‑150 (MQ=255)
                                                 aTGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTTCCCAGGTCCACCTACGTCCGTTTACttt                                                                                                     >  2:4763461/1‑150 (MQ=255)
                                                  tGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAGTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTc                                                                                                    <  2:5616637/150‑1 (MQ=255)
                                                  ggtgttaaTTCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTACGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGACGCTGTAAGCCAAATCTTTCCCCCTGCCCAGGGCCACCAACGTCCGTTTACTTTc                                                                                                    >  1:1410263/7‑150 (MQ=255)
                                                   gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                                                   <  1:4350991/150‑1 (MQ=255)
                                                    ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                                                  >  2:1487662/1‑150 (MQ=255)
                                                      gAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGtt                                                                                                <  1:838927/150‑1 (MQ=255)
                                                      gAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCGGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTTCTTTCCGtt                                                                                                >  1:1254816/1‑150 (MQ=255)
                                                         aTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaa                                                                                             <  1:6421374/150‑1 (MQ=255)
                                                          ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCGCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTCCGTCCGTTTACTTTCCGTTaaaa                                                                                            >  2:2407531/1‑150 (MQ=255)
                                                          ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACTTTCCGTTaaaa                                                                                            >  1:4708964/1‑150 (MQ=255)
                                                           tCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaaa                                                                                           <  2:4375822/150‑1 (MQ=255)
                                                            cACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAc                                                                                          <  2:5197002/150‑1 (MQ=255)
                                                              ccAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGc                                                                                        <  1:5703949/150‑1 (MQ=255)
                                                               cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCAGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCGTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGcg                                                                                       >  2:3222675/1‑149 (MQ=255)
                                                                aGACTATGTCTGGCGAATTAACTCCTTTGGGAAGTCCGCAGACACCGTCTGGAGAAGTCCGATCGGCAACGATTTGGCAGACGCTGTCTGCCTAATCTTTCCCCTTGCCCAGGTCCACCTTCGTCCGTTTACTTTCCGTTAAAAACGCTg                                                                                      <  2:5344420/150‑1 (MQ=255)
                                                                aGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAATGATTTGGCAGACGCTGTGTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTGCGGTTACTGTCCGGTAAAAACGCTg                                                                                      >  1:4525084/1‑150 (MQ=255)
                                                                aGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCCGGTCCACCTACGTCCGTGTACTTTCCGGTAACAACGCTg                                                                                      >  2:6419185/1‑150 (MQ=255)
                                                                 gACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCATGCCCAGGTCCACCTACGTCCGTTTACTTTCAGTTAAAAACGCTGa                                                                                     >  1:5515209/1‑150 (MQ=255)
                                                                     aTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGcc                                                                                 <  1:3948978/150‑1 (MQ=255)
                                                                          tGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAg                                                                            >  2:1438191/1‑150 (MQ=255)
                                                                              gAATTAACACCATTGGGAATTCCGCAGACACCGTATGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                        <  1:4622241/150‑1 (MQ=255)
                                                                               aaTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCtttt                                                                       <  1:5664502/150‑1 (MQ=255)
                                                                                aTTAACCCCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCAAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                      <  1:2904208/150‑1 (MQ=255)
                                                                                aTTAACACCCTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                      <  1:671667/150‑1 (MQ=255)
                                                                                aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                      <  2:4244861/150‑1 (MQ=255)
                                                                                aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                      <  2:865473/150‑1 (MQ=255)
                                                                                aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCCGACGCCCGCAGCTTTTa                                                                      <  1:3592951/150‑1 (MQ=255)
                                                                                 ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAACGATTTGGCAGACGCTGTCTGCCAACTCTTTCCCCTGGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                     <  1:1117683/150‑1 (MQ=255)
                                                                                 ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                     <  1:3858347/150‑1 (MQ=255)
                                                                                 ttAACACCATGGGGAATTCCGCAGACACCGTCTGCAGAAGTCCCATCGGCAAAGTTTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                     <  1:3447681/150‑1 (MQ=255)
                                                                                      acCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaca                                                                >  2:1002473/1‑148 (MQ=255)
                                                                                      acCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaa                                                                >  2:4933701/1‑150 (MQ=255)
                                                                                       ccATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCGGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                               <  1:1761455/150‑1 (MQ=255)
                                                                                         cttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                             <  2:1137149/149‑1 (MQ=255)
                                                                                         aTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                             >  2:1067419/1‑150 (MQ=255)
                                                                                         aTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCATTTATGAAAAaga                                                             >  1:3642754/1‑150 (MQ=255)
                                                                                          ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAagag                                                            <  1:1482174/150‑1 (MQ=255)
                                                                                            gggAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGAc                                                          <  1:4034974/150‑1 (MQ=255)
                                                                                              gAATTCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCAGACGCGGTCGGCCAAATCTTTCCCCGTGCCCCGGTCCACCAACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAGGAAAAAGAGACGc                                                        >  1:5425991/1‑150 (MQ=255)
                                                                                              gAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGc                                                        >  2:5865947/1‑150 (MQ=255)
                                                                                               aaTTCTGCAGACGGTGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCt                                                       >  2:5089879/1‑150 (MQ=255)
                                                                                               aaTTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCt                                                       <  1:6350904/150‑1 (MQ=255)
                                                                                                aTTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCACACGCGGTCTGCCCAATCTTTACCCAATCCCAGGTCCACCTACCTCCCTTTACATTAAGTTAAAAAAGCTGACCCCCCCAGCTTTTATTAAAAAGAGACGCTc                                                      >  1:5213538/1‑150 (MQ=255)
                                                                                                 ttCCGCAGACACCGTCTGCAGAATTCCCATCGGCAATGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACTTTCCGTTAAAAAAGCTGACGCCCGCAGCTGTTATGAAAAAGAGACGCTcc                                                     >  1:6596431/1‑150 (MQ=255)
                                                                                                      cAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt                                                <  2:53009/150‑1 (MQ=255)
                                                                                                         acacCGTCTCCATAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCAAAAACTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGc                                             <  1:3265053/150‑1 (MQ=255)
                                                                                                           acCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTg                                           <  1:1187064/150‑1 (MQ=255)
                                                                                                               tCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCt                                       <  2:1983170/150‑1 (MQ=255)
                                                                                                               tCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACATTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCt                                       >  2:1639349/1‑150 (MQ=255)
                                                                                                                   cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCTGGTCCACCTACGTCCGTTTACTTTCCGTTAAAATCGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                                   <  1:4244483/150‑1 (MQ=255)
                                                                                                                   cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                                   <  1:1809070/150‑1 (MQ=255)
                                                                                                                      aaTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGcc                                <  1:4190247/150‑1 (MQ=255)
                                                                                                                      aaTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGcc                                <  1:2845773/150‑1 (MQ=255)
                                                                                                                           ccATCGGCAAAGATTTGGCAGACGCTTTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCtt                           >  2:6212122/1‑149 (MQ=255)
                                                                                                                           ccATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTa                           <  2:3305765/150‑1 (MQ=255)
                                                                                                                             aTCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTaga                         >  2:266942/1‑150 (MQ=255)
                                                                                                                                  cAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGc                    <  2:2473312/150‑1 (MQ=255)
                                                                                                                                     aGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCaaa                 >  2:3580845/1‑150 (MQ=255)
                                                                                                                                      gATTTGGCAGACGCTGTCTGCCAAATATTTCCCCTTGCCCAGGTCCACCAACGTCCGTTTAAATTACGGTAAAAAAGCTGACGCCCCCAGCTTTTATGAAAAAAAGACGCTCCGCGGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAt                >  1:2344521/1‑150 (MQ=255)
                                                                                                                                       aTTTGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGGAAAAGAGACGCTCCGCTGGGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAtt               >  2:812379/1‑150 (MQ=255)
                                                                                                                                       aTTTGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAtt               >  1:911732/1‑150 (MQ=255)
                                                                                                                                        tttGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg              >  1:2080686/1‑150 (MQ=255)
                                                                                                                                           ggCGGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGa           <  2:4522134/150‑1 (MQ=255)
                                                                                                                                           ggCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGa           <  2:3161337/150‑1 (MQ=255)
                                                                                                                                           ggCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGa           <  2:5675573/150‑1 (MQ=255)
                                                                                                                                            gCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTg                                      <  2:5549559/122‑1 (MQ=255)
                                                                                                                                            gCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAc          >  2:2431656/1‑150 (MQ=255)
                                                                                                                                             cAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAcc         >  2:2933684/1‑150 (MQ=255)
                                                                                                                                             cAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAcc         >  1:1690864/1‑150 (MQ=255)
                                                                                                                                               gACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCa       >  2:3920575/1‑150 (MQ=255)
                                                                                                                                                aCGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAg      >  1:5903873/1‑150 (MQ=255)
                                                                                                                                                    tGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAAAGCGGCAAATTGCGACCCAGtttt  >  2:2785461/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTT  >  NC_000913/3368073‑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: