Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,261,933 T→A 6.3% Q289L (CAG→CTG)  PP_4634 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,261,9330TA 6.3% 1028.8 / 5.6 328Q289L (CAG→CTG) PP_4634alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (144/148);  new base A (24/3);  total (177/151)
Fisher's exact test for biased strand distribution p-value = 6.59e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

ATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACG  >  NC_002947/5261786‑5262081
                                                                                                                                                   |                                                                                                                                                    
aTCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGc                                                                                                                                                    >  2:6694550/1‑150 (MQ=255)
 tCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGcc                                                                                                                                                   <  2:594633/150‑1 (MQ=255)
 tCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGcc                                                                                                                                                   <  2:2794070/150‑1 (MQ=255)
    aaCACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCa                                                                                                                                                <  2:5804873/150‑1 (MQ=255)
    aaCACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCa                                                                                                                                                >  1:2450676/1‑150 (MQ=255)
    aaCACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCGTTGCGGCCCCCCCCTCCAGCGCCCACCCGCAATCCATGATGGGGCCCCGCCCCCCCCCGTCCACCCGCCCCAGCCGcc                                                                                                                                                >  1:308151/1‑149 (MQ=255)
     acacGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCGCCATCCATGATGGTGCCAGGCCCCCCGCCGTCCAGCCGCCCCGGCCGCaa                                                                                                                                               >  2:851378/1‑150 (MQ=255)
     acacGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCAGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCGTCCACCCACCCCAGCCGCaa                                                                                                                                               >  2:7132/1‑150 (MQ=255)
      cacGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAAc                                                                                                                                              <  2:3704527/150‑1 (MQ=255)
      cacGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAAc                                                                                                                                              <  2:5473633/150‑1 (MQ=255)
          gTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACt                                                                                                                                          <  1:3376764/150‑1 (MQ=255)
              ggcgggcgTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCaa                                                                                                                                      <  2:2359184/150‑1 (MQ=255)
               gcgggcgTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCCCGTTTCGGCCCCCCCAACCAGCGCCCACCCACAAAACAAGATGGTTCCAGGCCCCCCCCCATCCACCCACCCCAGCCGCCACCCCAGCaac                                                                                                                                     >  1:4057959/1‑149 (MQ=255)
                    cgTGAAGTTTGCTGCGCTGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCAGCCGCAACCACTGCAAGGCCAc                                                                                                                                <  1:5284200/150‑1 (MQ=255)
                     gTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCAcc                                                                                                                               <  2:5721337/150‑1 (MQ=255)
                      tGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCg                                                                                                                              >  2:1137952/1‑150 (MQ=255)
                      tGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCg                                                                                                                              <  1:3254038/150‑1 (MQ=255)
                      tGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCGTTTCGGCCCCCCCATACAGCGCCCACCCACAAAACAAGAAGGTGGCCCGCCCCCCCCCCACCAGCCACCCCAGCCGCCCCCACAGCAACGCCAccc                                                                                                                              >  1:82443/1‑149 (MQ=255)
                           tttGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCCACCACTGCAAGGCCACCGGCcac                                                                                                                         >  2:6161621/1‑150 (MQ=255)
                             tgctgcGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAg                                                                                                                       <  2:1840075/150‑1 (MQ=255)
                              gctgcGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCCCCGGCCCCAgg                                                                                                                      >  2:6460020/1‑150 (MQ=255)
                                tgcGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCt                                                                                                                    >  2:5142051/1‑150 (MQ=255)
                                tgcGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCGTTGCGGCCCCACGGTCCGGCGCCCACCCACAAACCCGGATGGTGCCAGGCCCCCCGCCGTCCAGCCGCCCCCGCCGCAACCACGGCCAGGCCCCCCGCCCCACGCt                                                                                                                    >  2:2536776/1‑150 (MQ=255)
                                 gcgcAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCGTTGCGGCCCCCCGATCCTGCGCCCCCCCCCAAACCATGATGGTGCCAGGCCCCCCCCCGTCCACCCGCCCCCGCCGCCACCCCCGCAAGGCCCCCCGCCCCAGGcgt                                                                                                                   >  1:5735965/1‑148 (MQ=255)
                                 gcgcAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCGTTGCGGCCCCACGCTCCTGCGCCCCCCCGCAATCCATGATGGTGCCAGGCCCCCCGCCTTACAGCCTCCCCAGCCGCAACCACGGCAAGGCCACCCGCCACAGGCtt                                                                                                                   >  2:1925514/1‑150 (MQ=255)
                                             gcgcgcAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATACTGCGCCCACCCGCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCGCCCCTGCCGCAACCCCTGCAAGGCCACCCGCCACAGGCTTTCGCGCAACCgg                                                                                                       >  2:5629351/1‑150 (MQ=255)
                                                 gcAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCAt                                                                                                   <  2:4330214/150‑1 (MQ=255)
                                                 gcAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCAt                                                                                                   <  2:3411147/150‑1 (MQ=255)
                                                 gcAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCAt                                                                                                   <  2:5007167/150‑1 (MQ=255)
                                                  cAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATg                                                                                                  <  1:5340852/150‑1 (MQ=255)
                                                      gCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGaaa                                                                                              <  1:5624951/150‑1 (MQ=255)
                                                       cTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGAACCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGaaaa                                                                                             >  2:4932197/1‑150 (MQ=255)
                                                         tGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAAc                                                                                           <  1:1014097/150‑1 (MQ=255)
                                                         tGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAAc                                                                                           <  1:3714692/150‑1 (MQ=255)
                                                          gagaTTGCCGGGGCCGCGTTGCGGCCCCACGATCCAGCGCCCACCCACAATCCATGATGGTGCCAGGCCCCCCCCCGTCCCGCCTCCCCGGCCGCAACCACAGCAACGCCCCCCGCCCCAGGCTTTCGCGAAACCGGGCAGGGAAAAAcg                                                                                          >  2:3516618/1‑150 (MQ=255)
                                                                 ccGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCGGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAACGCGAAATgg                                                                                  >  2:1503810/1‑150 (MQ=255)
                                                                  cGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATgg                                                                                  >  2:6097164/1‑150 (MQ=255)
                                                                    gggCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAACTGGcc                                                                                >  2:6346782/1‑150 (MQ=255)
                                                                        cGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAAtt                                                                            >  1:2837336/1‑150 (MQ=255)
                                                                        cGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCACGAAACCGGTCATGGAAAAACGCGAAATGGCCAAtt                                                                            >  1:2938329/1‑150 (MQ=255)
                                                                         gCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAAttt                                                                           >  2:3110413/1‑150 (MQ=255)
                                                                         gCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAAttt                                                                           >  1:7086770/1‑150 (MQ=255)
                                                                          cTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTc                                                                          <  2:5269652/150‑1 (MQ=255)
                                                                           tttGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTcc                                                                         <  1:2097882/150‑1 (MQ=255)
                                                                           tttGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTcc                                                                         >  1:4346401/1‑150 (MQ=255)
                                                                             tGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTc                                                                       <  2:3007190/150‑1 (MQ=255)
                                                                               cGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCgg                                                                     <  2:2216236/150‑1 (MQ=255)
                                                                                ggCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGc                                                                    >  1:4827553/1‑150 (MQ=255)
                                                                                ggCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGACCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGc                                                                    <  2:2101934/150‑1 (MQ=255)
                                                                                 gCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGcc                                                                   <  2:6872826/150‑1 (MQ=255)
                                                                                 gCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGcc                                                                   >  1:4436604/1‑150 (MQ=255)
                                                                                 gCCCCACGATCCTGCGCCCACCCTCAATCCAGGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGcc                                                                   <  1:3760347/150‑1 (MQ=255)
                                                                                   cccACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCga                                                                 <  2:4019064/150‑1 (MQ=255)
                                                                                    ccACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCgag                                                                >  1:221184/1‑150 (MQ=255)
                                                                                    ccACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCgag                                                                >  1:5490269/1‑150 (MQ=255)
                                                                                    ccACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCAGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCgag                                                                >  1:1606371/1‑150 (MQ=255)
                                                                                     cACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCgaga                                                               >  1:7286670/1‑150 (MQ=255)
                                                                                      aCGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGat                                                              <  1:7447384/150‑1 (MQ=255)
                                                                                             tGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGACATGGCCCATTTCCTCGGCCGAGATATCcaca                                                       >  2:4705743/1‑150 (MQ=255)
                                                                                              gcgcCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAg                                                      >  1:4596831/1‑150 (MQ=255)
                                                                                                  ccACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCCATTTCCTCGGCCGAGATATCCACAGGTGc                                                  >  1:5828107/1‑150 (MQ=255)
                                                                                                  ccACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGc                                                  >  1:6646851/1‑150 (MQ=255)
                                                                                                     cccTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAAt                                               <  1:6092437/150‑1 (MQ=255)
                                                                                                      ccTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATg                                              >  1:5000046/1‑150 (MQ=255)
                                                                                                       cTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATgc                                             >  1:7156492/1‑150 (MQ=255)
                                                                                                       cTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATgc                                             >  2:3222020/1‑150 (MQ=255)
                                                                                                        tCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAAttt                                                                           >  1:3873549/1‑119 (MQ=255)
                                                                                                        tCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAAttt                                                                           <  2:3873549/119‑1 (MQ=255)
                                                                                                           aTCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAg                                         <  2:1901341/150‑1 (MQ=255)
                                                                                                           aTCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAg                                         <  2:7013220/150‑1 (MQ=255)
                                                                                                           aTCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAg                                         <  1:6694550/150‑1 (MQ=255)
                                                                                                             ccATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGt                                       <  1:544891/150‑1 (MQ=255)
                                                                                                              cATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTg                                      >  1:6114317/1‑150 (MQ=255)
                                                                                                               atgatgGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGc                                     >  1:2404098/1‑150 (MQ=255)
                                                                                                                   tgGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTc                                 >  1:1784468/1‑150 (MQ=255)
                                                                                                                       gCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACCGGTGCAATGCGCAGGTGCTTTCGCTc                             <  2:2619531/150‑1 (MQ=255)
                                                                                                                       gCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTc                             <  1:4954139/150‑1 (MQ=255)
                                                                                                                        ccAGGCCCCCCTCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCTGCCACAGGCTTTCGCGAAACCGTTCATGGAAAAAAGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCg                            <  2:1293727/150‑1 (MQ=255)
                                                                                                                        ccAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCg                            <  2:910664/150‑1 (MQ=255)
                                                                                                                        ccAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCg                            >  2:243179/1‑150 (MQ=255)
                                                                                                                        ccAGGCCCCCCGCCTTCCAGCCTCCCCTGCCACAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTcgcgcg                            >  1:6557315/1‑150 (MQ=255)
                                                                                                                              cccccGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCAccc                      >  1:1189393/1‑150 (MQ=255)
                                                                                                                              cccccGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCAccc                      >  2:1729728/1‑150 (MQ=255)
                                                                                                                               ccccGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCt                     <  2:4870536/150‑1 (MQ=255)
                                                                                                                                cccGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTg                    >  2:1011797/1‑150 (MQ=255)
                                                                                                                                 ccGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGt                   >  2:4507231/1‑150 (MQ=255)
                                                                                                                                 ccGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGt                   >  1:2628996/1‑150 (MQ=255)
                                                                                                                                 ccGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGt                   >  1:909222/1‑150 (MQ=255)
                                                                                                                                 ccGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGt                   >  2:6083981/1‑150 (MQ=255)
                                                                                                                                  cGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTa                  <  2:94833/150‑1 (MQ=255)
                                                                                                                                   gCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAg                 >  1:2327608/1‑150 (MQ=255)
                                                                                                                                     cTTCCACCCTCCCCACCCGCAACCACCCCAAGCCCACCGGCCACACGCTTTCGCCAAACCCGTCATGCAAAAACCCGAAATGGCCAATTTCCACGGCCGAGATATCCACACGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGAt               <  2:7321272/150‑1 (MQ=255)
                                                                                                                                      ttCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATa              >  1:6427062/1‑150 (MQ=255)
                                                                                                                                        ccAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAAc            >  1:6238742/1‑150 (MQ=255)
                                                                                                                                          aGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAAccc          >  2:3910563/1‑150 (MQ=255)
                                                                                                                                            ccTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCcaa        >  2:1629160/1‑150 (MQ=255)
                                                                                                                                            ccTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCcaa        <  1:3524964/150‑1 (MQ=255)
                                                                                                                                              tCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCca                                                         <  2:2767568/99‑1 (MQ=255)
                                                                                                                                              tCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCca                                                         >  1:2767568/1‑99 (MQ=255)
                                                                                                                                              tCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTGCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCGGTCGATCACCcaaca      >  2:3780372/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCAAACACGGTCTGGCGGGCGTGAAGTTTGCTGCGCAGCCGCTTGCGCGCAAGGCTTGAGATTGCCGGGGCCGCTTTGCGGCCCCACGATCCTGCGCCCACCCTCAATCCATGATGGTGCCAGGCCCCCCGCCTTCCAGCCTCCCCTGCCGCAACCACTGCAAGGCCACCGGCCACAGGCTTTCGCGAAACCGGTCATGGAAAAACGCGAAATGGCCAATTTCCTCGGCCGAGATATCCACAGGTGCAATGCGCAGGTGCTTTCGCTCGGCACCCTGTAGATAACCCAACAAACG  >  NC_002947/5261786‑5262081

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: