Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,219 G→A 6.4% 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.4% 1135.8 / 44.9 343T131T (ACG→ACAyhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base G (150/170);  new base A (18/4);  total (168/175)
Fisher's exact test for biased strand distribution p-value = 1.63e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

Only 100 of 348 total aligned reads displayed.

TTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGTT  >  NC_000913/3368070‑3368368
                                                                                                                                                     |                                                                                                                                                     
ttCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCCGACg                                                                                                                                                       >  2:3128078/1‑150 (MQ=255)
    ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAAATCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt                                                                                                                                                   >  2:2009342/1‑150 (MQ=255)
    ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTAGAAATTCGCCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCCGACACCGTCTGCAGAAGTCCCCTCGGCAACGATTTGGCAGACGCAGt                                                                                                                                                   >  2:5622583/1‑150 (MQ=255)
     tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTc                                                                                                                                                  >  2:7198545/1‑150 (MQ=255)
      gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                                                 <  1:4139575/150‑1 (MQ=255)
      gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGGCt                                                                                                                                                 >  2:6432373/1‑150 (MQ=255)
        aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                               <  2:3626480/150‑1 (MQ=255)
         aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGcc                                                                                                                                              >  1:3246234/1‑150 (MQ=255)
          aaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCa                                                                                                                                             <  2:2129848/150‑1 (MQ=255)
           acaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACACTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCaa                                                                                                                                            <  2:5139456/148‑1 (MQ=255)
           aaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGAAAAGATTTGGCAGACAGCGTCTGCCaa                                                                                                                                            <  2:3014076/150‑1 (MQ=255)
                   gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                                                    >  1:2846563/1‑150 (MQ=255)
                    ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGTCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTccc                                                                                                                                   >  2:2335156/1‑150 (MQ=255)
                    ccAAATGAGGCTTTTTTACCTCTTTGTTCAACATGTTGAAATTAACCAGACAAGGTCTGGCGAAGTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCTGACGCTGTCTGCCAAATCTTTccc                                                                                                                                   >  2:5212708/1‑150 (MQ=255)
                    ccAAATGAGGCTTGTTTACCTCTTGTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTGTccc                                                                                                                                   >  2:6303101/1‑150 (MQ=255)
                               ttttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGt                                                                                                                        >  1:5388054/1‑150 (MQ=255)
                                 ttttACCACTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGCTTTGGCAGACGCGGTCTGCCAAATCTTTCCCCTTGCCCAGGTcc                                                                                                                      <  1:4731397/150‑1 (MQ=255)
                                  tttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCa                                                                                                                     <  2:6973532/150‑1 (MQ=255)
                                  tttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCa                                                                                                                     >  1:3391945/1‑150 (MQ=255)
                                   ttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCTAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                    <  1:4651829/150‑1 (MQ=255)
                                   ttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                    >  1:1136482/1‑150 (MQ=255)
                                    tACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGTCACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAca                                                                                                                   >  2:2044714/1‑149 (MQ=255)
                                     aCCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCt                                                                                                                  >  2:4819429/1‑150 (MQ=255)
                                       ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCCTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                <  2:4232577/150‑1 (MQ=255)
                                       ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                <  2:3904302/150‑1 (MQ=255)
                                       ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCACAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                                <  1:1447302/150‑1 (MQ=255)
                                        tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACg                                                                                                               <  2:4109749/150‑1 (MQ=255)
                                           tttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTcc                                                                                                            <  1:6708275/150‑1 (MQ=255)
                                           tttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTcc                                                                                                            <  1:3732211/150‑1 (MQ=255)
                                              ttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGtt                                                                                                         <  2:2799771/150‑1 (MQ=255)
                                                  aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACt                                                                                                     >  2:4050149/1‑150 (MQ=255)
                                                     tGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTc                                                                                                  <  2:6644779/150‑1 (MQ=255)
                                                      gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                                                 <  1:6641888/150‑1 (MQ=255)
                                                      gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                                                 <  2:5814566/150‑1 (MQ=255)
                                                        tGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGTCGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGt                                                                                               >  1:4217714/1‑150 (MQ=255)
                                                        tGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGt                                                                                               <  1:2578140/150‑1 (MQ=255)
                                                         gAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGCCCACCTACGTCCGTTTACTTTCCGtt                                                                                              >  2:6548654/1‑150 (MQ=255)
                                                          aaaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTa                                                                                             >  2:4181463/1‑150 (MQ=255)
                                                          aaaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCGTGCCCAGGGCCACCTACGTCCGTTTACTTTCCGTTa                                                                                             >  1:4426371/1‑150 (MQ=255)
                                                             ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGGATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaat                                                                                          >  2:538182/1‑149 (MQ=255)
                                                             ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCGGCCAAATCTGTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACTTTCCGTTcaac                                                                                          >  2:3408866/1‑149 (MQ=255)
                                                               cACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAc                                                                                        >  2:4922516/1‑150 (MQ=255)
                                                               cACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAc                                                                                        <  1:2843049/150‑1 (MQ=255)
                                                               cACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTGGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAc                                                                                        <  1:3798041/150‑1 (MQ=255)
                                                                aCCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTACCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaaaag                                                                                       >  1:3943789/1‑148 (MQ=255)
                                                                 ccAGACAATGTCTGGCGAATTAACACCCTGGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGc                                                                                      <  2:180604/150‑1 (MQ=255)
                                                                 ccAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGc                                                                                      >  1:452018/1‑150 (MQ=255)
                                                                  cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                     >  1:5259598/1‑150 (MQ=255)
                                                                  cAGACAAGGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                     <  2:140727/150‑1 (MQ=255)
                                                                     aCAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGAc                                                                                  <  1:3329237/150‑1 (MQ=255)
                                                                              ggCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGGTAAAAACGCTGACGCCCGCAGc                                                                         >  1:5150085/1‑150 (MQ=255)
                                                                                 gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                      <  2:6130496/150‑1 (MQ=255)
                                                                                 gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                      >  2:3787816/1‑150 (MQ=255)
                                                                                 gAATTAACACCATTGGGAATACCGCAGACACCGTCTGCAGAATGCCCATCGGCAACGATTTGGCAGACGATGTCTGCAAAATCATTCCCCTTGCCCTGGTCCACCGACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                      <  1:5794996/150‑1 (MQ=255)
                                                                                    ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                   <  2:950166/150‑1 (MQ=255)
                                                                                    ttAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                   <  1:7184652/150‑1 (MQ=255)
                                                                                     tAACACCATTGGGAATTCCGCAGACACCGTCTGCACAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATg                                                                  <  2:6030835/150‑1 (MQ=255)
                                                                                       acacCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaa                                                                <  1:2007016/150‑1 (MQ=255)
                                                                                       acacCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCGGTCTGCCAAATCTTTACCCATGCCCCGGTCCACCAACGGCCGTTTAATTTACGTTAAAAACGCTGACGCCCCCAGCTTTTATGaa                                                                >  1:853326/1‑150 (MQ=255)
                                                                                         acCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGACGCTGTCTGCCAAAAATTTACCCATTCCCAGGTCCACCTACGTCCCGTTAAATTACGTTAAAAAAGCAGACGCCCGCAGCTTTTATGaaaa                                                              >  1:3032336/1‑150 (MQ=255)
                                                                                          ccATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCGTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                             <  1:1574643/150‑1 (MQ=255)
                                                                                            aTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                           <  1:793313/150‑1 (MQ=255)
                                                                                            aTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTATGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                           >  1:2385563/1‑150 (MQ=255)
                                                                                             ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATATTTCCCCATGCCCAGGTCCACCTACGTCCGTTTAAATTGCGTTTAAAAAGCGGGCGCCCGCAGCTTTTATTAAAAagtg                                                          >  1:4950008/1‑148 (MQ=255)
                                                                                                ggAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACg                                                       >  2:1226424/1‑150 (MQ=255)
                                                                                                 gAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGc                                                      >  2:3584272/1‑150 (MQ=255)
                                                                                                  aaTTCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGTTTTGGCAGACGCTGTCTGCCAAATCTTTACCCATGCCCAGGTCCACCTACGTCCGTTTAATTTCCGGTAAAAACGCTGACGCCCGCCGCTGTTATGAAAAAGAGTCGca                                                     >  1:5502747/1‑149 (MQ=255)
                                                                                                  aaTTCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCt                                                     >  2:2571943/1‑150 (MQ=255)
                                                                                                    ttCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                    <  1:1064072/85‑1 (MQ=37)
                                                                                                      ccGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGAAGCTGTCTGCCAAAACTTTACCCATGCCCAGGTCCACCTACGTCCGGTTACTTTGCGTTAAAAAAGCTGACGCCCCCCGCGTTTATGAAAAAGAGAAGCTCCGc                                                 >  1:2196621/1‑150 (MQ=255)
                                                                                                       cGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCt                                                <  2:2473346/150‑1 (MQ=255)
                                                                                                        gCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtg                                               >  2:5858051/1‑150 (MQ=255)
                                                                                                        gCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAAGCTGTCTGCCAAATATTTCCCCTTGCCCAGGTACACCTAAGTCCGTTTACTTTACGTTAAAAACGCAGACGCCCGCAGCTTTTATGTAAAAGAGACGCTCCGCtg                                               >  1:6237473/1‑150 (MQ=255)
                                                                                                              acCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTg                                         >  2:2683372/1‑150 (MQ=255)
                                                                                                                cGTCTGCAGAATTCCCATCGGCAAATATTTGGCAGACGCTGTCTGCCAAAACTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt                                       <  1:1229841/150‑1 (MQ=255)
                                                                                                                cGTCTGCAGAATTCCCATCGGCAAAGATGTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt                                       <  1:652867/150‑1 (MQ=255)
                                                                                                                cGTCTGAAGAATTCCCATCTGCAAAGATTTGGCAGACTCTGTCTGACAAATCTTTCCCCTTTCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCTGCTTTTATGAAAAAGAGACGCTCCTCTGTGGCTGGt                                       <  2:3312440/150‑1 (MQ=255)
                                                                                                                  tCTGCAGAATTCCCATCGGCAAAGATGTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCt                                     <  1:5862644/150‑1 (MQ=255)
                                                                                                                      cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAGCGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                                 <  2:5973134/150‑1 (MQ=255)
                                                                                                                      cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                                 <  2:2368112/150‑1 (MQ=255)
                                                                                                                      cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCAAAATCTTTCCCCTTGCCCTGGTCCACCTACTTCCGTTTACTTTCCGTTATAATCGCTGACGCCCGCAGCTTTTATGAAATAGTGTCGCTCAGCTGTGGATGGTCTGTTc                                 <  2:775773/150‑1 (MQ=255)
                                                                                                                               cATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCtcg                        >  2:1103126/1‑148 (MQ=255)
                                                                                                                                 tCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTagag                      <  1:696962/150‑1 (MQ=255)
                                                                                                                                  ggggaaaGATTTGGCAGACAGCGTCTGCCGAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGc                     >  1:4828136/5‑150 (MQ=255)
                                                                                                                                  cGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGc                     <  1:3584272/150‑1 (MQ=255)
                                                                                                                                   ggCAAAGATTTGGCAGACGCTGTCTGCCAAATCGTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCg                    <  1:2765990/150‑1 (MQ=255)
                                                                                                                                    gCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCgg                   <  2:4422190/150‑1 (MQ=255)
                                                                                                                                    gCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCgg                   >  2:1580173/1‑150 (MQ=255)
                                                                                                                                          aTTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAtt             >  2:1719676/1‑150 (MQ=255)
                                                                                                                                          aTTTGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAtt             >  2:5483680/1‑150 (MQ=255)
                                                                                                                                          aTTTGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAtt             >  1:1638805/1‑150 (MQ=255)
                                                                                                                                            ttGGCAGACAGCGTCTGCCAAATCTTTGCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc           >  2:4564420/1‑150 (MQ=255)
                                                                                                                                               gCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAc        >  2:4724484/1‑150 (MQ=255)
                                                                                                                                                cAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTTTGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAcc       >  1:2419153/1‑150 (MQ=255)
                                                                                                                                                 aGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGAccc      >  1:5640672/1‑150 (MQ=255)
                                                                                                                                                  gACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAACAAGAGACGCTCCGCGGGGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCa     >  1:1828953/1‑150 (MQ=255)
                                                                                                                                                  gACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCa     <  2:6834248/150‑1 (MQ=255)
                                                                                                                                                  gACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCa     <  1:6142125/150‑1 (MQ=255)
                                                                                                                                                   aCGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAg    >  2:899609/1‑150 (MQ=255)
                                                                                                                                                     gCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAGtt  >  1:6924563/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: