Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,221 T→C 6.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 6.5% 1075.2 / 47.9 346L132P (CTG→CCG) ‡yhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base T (144/167);  new base C (18/4);  total (171/175)
Fisher's exact test for biased strand distribution p-value = 1.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-02

Only 100 of 354 total aligned reads displayed.

CTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTTTT  >  NC_000913/3368072‑3368370
                                                                                                                                                     |                                                                                                                                                     
   tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGTATTAACACCATTTGGAATTCCGCAGACACCGTCTGCAGAATTCACAAAGGCAAAGATTTTGGAGACGCTggc                                                                                                                                                    >  1:3312440/1‑148 (MQ=255)
    gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGGCt                                                                                                                                                   >  2:6432373/1‑150 (MQ=255)
      aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                                 <  2:3626480/150‑1 (MQ=255)
       aaaaaCTTACGCCAAATGAGTCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCTAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGcc                                                                                                                                                <  1:407674/150‑1 (MQ=255)
       aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGcc                                                                                                                                                >  1:3750598/1‑150 (MQ=255)
        aaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCa                                                                                                                                               <  2:2129790/150‑1 (MQ=255)
         aaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGCCAATGTCTGGCGAATTAACACCAATGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCaa                                                                                                                                              <  2:3012167/150‑1 (MQ=255)
           aCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAAt                                                                                                                                            <  1:993925/150‑1 (MQ=255)
                 gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                                                      >  1:2846563/1‑150 (MQ=255)
                      aTGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAAGGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTg                                                                                                                                 <  2:6548204/150‑1 (MQ=255)
                        gAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGcc                                                                                                                               >  2:3035302/1‑150 (MQ=255)
                           gCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAg                                                                                                                            >  1:6320795/1‑150 (MQ=255)
                             ttttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGt                                                                                                                          >  2:560477/1‑150 (MQ=255)
                               ttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTcc                                                                                                                        <  1:4731396/150‑1 (MQ=255)
                                tttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCa                                                                                                                       <  2:6973532/150‑1 (MQ=255)
                                 ttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                      <  1:359103/150‑1 (MQ=255)
                                 ttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                      >  1:3853484/1‑150 (MQ=255)
                                 ttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                      >  1:1136482/1‑150 (MQ=255)
                                     ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCCTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                  <  2:4232577/150‑1 (MQ=255)
                                     ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                  >  1:3472308/1‑150 (MQ=255)
                                     ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                  <  2:3904302/150‑1 (MQ=255)
                                     ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCACAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                  <  1:1447302/150‑1 (MQ=255)
                                      tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACg                                                                                                                 <  2:4109749/150‑1 (MQ=255)
                                       ctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCGGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGt                                                                                                                >  2:5198369/1‑150 (MQ=255)
                                        ttttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTc                                                                                                               >  2:2005853/1‑150 (MQ=255)
                                        ttttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTc                                                                                                               >  2:2364885/1‑150 (MQ=255)
                                         tttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTcc                                                                                                              <  1:6708275/150‑1 (MQ=255)
                                            ttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATGGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAACGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGtt                                                                                                           <  2:2799670/150‑1 (MQ=255)
                                                aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACt                                                                                                       >  2:4050149/1‑150 (MQ=255)
                                                   tGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTc                                                                                                    <  2:6644779/150‑1 (MQ=255)
                                                    gTTGAAATTCACCAGACACTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                                                   <  2:7162902/150‑1 (MQ=255)
                                                     ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                                                  >  2:2016124/1‑150 (MQ=255)
                                                     ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTACCCATGCCCAGGTCCACCTACGTCCGGTTACTTTACg                                                                                                  >  1:2520691/1‑150 (MQ=255)
                                                      tGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGt                                                                                                 <  2:4686018/150‑1 (MQ=255)
                                                      tGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTGTACTTTCCGt                                                                                                 >  1:7149475/1‑150 (MQ=255)
                                                       gataTTCACCAGATAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCTTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGCTCCACCTACGTACGTTTACTTTCCGtt                                                                                                <  2:5472321/147‑1 (MQ=255)
                                                        aaaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTa                                                                                               >  2:4181463/1‑150 (MQ=255)
                                                        aaaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTa                                                                                               <  1:5339005/150‑1 (MQ=255)
                                                        aaaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTAGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTa                                                                                               >  2:2444351/1‑150 (MQ=255)
                                                           ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCGGCCAAATCTGTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACTTTCCGTTcaac                                                                                            >  2:3408866/1‑149 (MQ=255)
                                                             cACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAc                                                                                          <  1:1257987/150‑1 (MQ=255)
                                                             cACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAc                                                                                          <  1:2843049/150‑1 (MQ=255)
                                                                cAGACAAGGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                       <  2:140727/150‑1 (MQ=255)
                                                                       tGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGccc                                                                                <  2:6286570/150‑1 (MQ=255)
                                                                        gTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCg                                                                               <  1:5217799/150‑1 (MQ=255)
                                                                           tGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTTCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCAGACGCCCGccg                                                                            >  1:3388389/1‑148 (MQ=255)
                                                                             gCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCt                                                                          >  2:2330891/1‑150 (MQ=255)
                                                                             gCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCGGCCAAATCTTTCCCCTTGCCCAGGTCCACCAACGTCCGGTTACTGTCCGTTAAAAACGCTGACGCCCGCCGCt                                                                          >  1:4150196/1‑150 (MQ=255)
                                                                               gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                        >  2:3787816/1‑150 (MQ=255)
                                                                               gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                        <  2:6130496/150‑1 (MQ=255)
                                                                               gAATTAACACCATTGGGAATACCGCAGACACCGTCTGCAGAATGCCCATCGGCAACGATTTGGCAGACGATGTCTGCAAAATCATTCCCCTTGCCCTGGTCCACCGACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                        <  1:5794996/150‑1 (MQ=255)
                                                                                 aTTACCACCATTGGGAAGTCCGCAGACACCGTCTGCAGAATTCCCATCGGCACAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                      <  1:1864282/150‑1 (MQ=255)
                                                                                 aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                      <  1:797994/150‑1 (MQ=255)
                                                                                  ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAATATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                     <  2:1592176/150‑1 (MQ=255)
                                                                                  ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                     <  2:950166/150‑1 (MQ=255)
                                                                                     acacCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaa                                                                  <  1:2007016/150‑1 (MQ=255)
                                                                                       acCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGACGCTGTCTGCCAAAAATTTACCCATTCCCAGGTCCACCTACGTCCCGTTAAATTACGTTAAAAAAGCAGACGCCCGCAGCTTTTATGaaaa                                                                >  1:3032336/1‑150 (MQ=255)
                                                                                        ccATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                               <  1:1574700/150‑1 (MQ=255)
                                                                                        ccATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                               <  2:6016475/150‑1 (MQ=255)
                                                                                             gggAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGAc                                                          >  2:174161/1‑150 (MQ=255)
                                                                                              ggAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACg                                                         >  2:1226642/1‑150 (MQ=255)
                                                                                               gAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGc                                                        >  2:3584272/1‑150 (MQ=255)
                                                                                                aaTTCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGTTTTGGCAGACGCTGTCTGCCAAATCTTTACCCATGCCCAGGTCCACCTACGTCCGTTTAATTTCCGGTAAAAACGCTGACGCCCGCCGCTGTTATGAAAAAGAGTCGca                                                       >  1:5502747/1‑149 (MQ=255)
                                                                                                  ttCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                      <  1:1064072/85‑1 (MQ=37)
                                                                                                  ttCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTGGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTcc                                                     <  1:4943133/150‑1 (MQ=255)
                                                                                                    ccGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGAAGCTGTCTGCCAAAACTTTACCCATGCCCAGGTCCACCTACGTCCGGTTACTTTGCGTTAAAAAAGCTGACGCCCCCCGCGTTTATGAAAAAGAGAAGCTCCGc                                                   >  1:2196621/1‑150 (MQ=255)
                                                                                                      gCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtg                                                 >  2:5802164/1‑150 (MQ=255)
                                                                                                      gCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCAGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGcag                                                 >  2:4877554/1‑148 (MQ=255)
                                                                                                        aGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGACGCTGTCGGCCAAAAATTTCCCCTTGCCCCGGGCCACCAACCTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAAAGACGCACCGcgggg                                               >  1:5135698/1‑146 (MQ=255)
                                                                                                           cacCGTCTGCAGAATTCCCATCGGCAAAGCTTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCt                                            <  2:2668329/150‑1 (MQ=255)
                                                                                                            cccGTCTGCCGCATTCCCATCGGCAACGATTGGGCAGCCGCGGTCTGCCAAATCTTTCCCCTCGCCCAGGTCCACCTACGTCCGTTTCCGGTCCGCTAAAAACGCTGACGCCCGCCGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTg                                           <  1:7225997/149‑1 (MQ=255)
                                                                                                             ccGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTgg                                          >  2:884106/1‑150 (MQ=255)
                                                                                                              cGTCTGCAGAATTCCCATCGGCAAAGATGTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt                                         <  1:652867/150‑1 (MQ=255)
                                                                                                              cGTCTGAAGAATTCCCATCTGCAAAGATTTGGCAGACTCTGTCTGACAAATCTTTCCCCTTTCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCTGCTTTTATGAAAAAGAGACGCTCCTCTGTGGCTGGt                                         <  2:3312440/150‑1 (MQ=255)
                                                                                                                tCTGCAGAATTCCCATCGGCAAAGATGTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCt                                       <  1:5862644/150‑1 (MQ=255)
                                                                                                                  tGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt                                     >  2:583073/1‑150 (MQ=255)
                                                                                                                    cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAGGAGACGCTCCGCTGTGGCTGGTCTGTTc                                   <  2:2486599/150‑1 (MQ=255)
                                                                                                                    cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCAAAATCTTTCCCCTTGCCCTGGTCCACCTACTTCCGTTTACTTTCCGTTATAATCGCTGACGCCCGCAGCTTTTATGAAATAGTGTCGCTCAGCTGTGGATGGTCTGTTc                                   <  2:775773/150‑1 (MQ=255)
                                                                                                                             cATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTag                          <  1:4627859/150‑1 (MQ=255)
                                                                                                                             cATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTag                          >  2:5949030/1‑150 (MQ=255)
                                                                                                                             cATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCtcg                          >  2:1103126/1‑148 (MQ=255)
                                                                                                                                ggggaaaGATTTGGCAGACAGCGTCTGCCGAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGc                       >  1:4828065/5‑150 (MQ=255)
                                                                                                                                 ggCAAAGATTTGGCAGACGCTGTCTGCCAAATCGTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCg                      <  1:2765990/150‑1 (MQ=255)
                                                                                                                                  gCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCgg                     <  2:4422190/150‑1 (MQ=255)
                                                                                                                                       gATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAt                <  1:6642237/150‑1 (MQ=255)
                                                                                                                                        aTTTGGCAGACGCTGTCTGCCAAAAATTTACCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCATACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCCCAGTGGCTGGTATGTTAGCCCCCTAGATCGGCaaaat               >  1:2364533/1‑148 (MQ=255)
                                                                                                                                        aTTTGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAtt               >  1:4482622/1‑150 (MQ=255)
                                                                                                                                        aTTTGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAtt               >  1:1638805/1‑150 (MQ=255)
                                                                                                                                         tttGGCAGACAGCGTCTGCCAAATCTTTGCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg              >  1:6964399/1‑150 (MQ=255)
                                                                                                                                            ggCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAAGCTGACGCCCGCAGCTTTTATGAAAAAGATACGCTCCGCTGTGGCTGGTCTGTTAGCCAGCTAGAGCGGCAAATTGCGa           >  1:4297089/1‑150 (MQ=255)
                                                                                                                                             gCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCGGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAGTGCGAc          >  2:6824186/1‑150 (MQ=255)
                                                                                                                                              cAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTTTGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAcc         >  1:2419153/1‑150 (MQ=255)
                                                                                                                                              cAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAcc         <  2:5094921/150‑1 (MQ=255)
                                                                                                                                               aGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAccc        >  1:5640672/1‑150 (MQ=255)
                                                                                                                                                gACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCa       <  1:6142125/150‑1 (MQ=255)
                                                                                                                                                gACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCa       >  2:4771415/1‑150 (MQ=255)
                                                                                                                                                   gCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGtt    >  1:6924563/1‑150 (MQ=255)
                                                                                                                                                   gCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGtt    >  2:5311789/1‑150 (MQ=255)
                                                                                                                                                     tGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGtttt  >  2:4519160/1‑150 (MQ=255)
                                                                                                                                                     tGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCGAGAGCGGCAAATTGCGACCCAGtttt  >  1:3231253/1‑150 (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: