Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,412,796 A→G 5.8% R183R (CGA→CGG) ‡ PP_3020 → serine/threonine protein phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,412,7960AG 5.8% 774.2 / 22.6 243R183R (CGA→CGG) ‡PP_3020serine/threonine protein phosphatase
Reads supporting (aligned to +/- strand):  ref base A (105/123);  new base G (8/6);  total (113/130)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 248 total aligned reads displayed.

GGCCTGCTGCATGCAGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAAC  >  NC_002947/3412647‑3412945
                                                                                                                                                     |                                                                                                                                                     
ggCCTGCTGCATGCAGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGa                                                                                                                                                       <  1:5080605/150‑1 (MQ=255)
ggCCTGCTGCATGCAGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGa                                                                                                                                                       <  1:5832976/150‑1 (MQ=255)
    tgctgcATGCAGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGccc                                                                                                                                                   >  2:2226019/1‑150 (MQ=255)
    tgctgcATGCAGACAGCCCGTGTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGccc                                                                                                                                                   <  2:868851/150‑1 (MQ=255)
      ctgcATGCAGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCTCGCAGGCCctg                                                                                                                                                 <  1:6155391/150‑1 (MQ=255)
        gcatgcaGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCctgct                                                                                                                                               <  2:2105906/150‑1 (MQ=255)
        gcatgcaGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGAGCAGGGCTCGCAGGCCctgca                                                                                                                                               <  1:2626544/150‑2 (MQ=255)
        gcatgcaGACAGCCCGTTTGCCGACTGGACAGACTTACACTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCctgct                                                                                                                                               <  2:6513995/150‑1 (MQ=255)
         catgcaGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCctgctg                                                                                                                                              <  1:2124311/150‑1 (MQ=255)
         catgcaGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCctgctg                                                                                                                                              <  2:3972433/150‑1 (MQ=255)
          atgcaGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCTCGCAGGCCctgcacc                                                                                                                                             <  1:3589739/150‑4 (MQ=255)
          atgcaGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTgg                                                                                                                                             <  2:1201246/150‑1 (MQ=255)
          atgcaGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTgg                                                                                                                                             >  1:5538445/1‑150 (MQ=255)
             caGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCTCGCAGGCCctgcaccggt                                                                                                                                          <  1:4573141/150‑7 (MQ=255)
                aCAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGcc                                                                                                                                       <  2:81112/150‑1 (MQ=255)
                      cGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACAccc                                                                                                                                 >  2:5490805/1‑150 (MQ=255)
                        tttGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCg                                                                                                                               <  2:1699137/150‑1 (MQ=255)
                           gCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTg                                                                                                                            <  1:6118795/150‑1 (MQ=255)
                            ccGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGc                                                                                                                                        >  2:2690510/1‑137 (MQ=255)
                            ccGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGc                                                                                                                                        <  1:2690510/137‑1 (MQ=255)
                              gACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTg                                                                                                                         <  2:3373738/150‑1 (MQ=255)
                              gACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTg                                                                                                                         >  2:1261258/1‑150 (MQ=255)
                                cTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGa                                                                                                                       <  2:2793945/150‑1 (MQ=255)
                                cTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGa                                                                                                                       <  1:3637611/150‑1 (MQ=255)
                                cTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGa                                                                                                                       <  1:3266906/150‑1 (MQ=255)
                                 tGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAg                                                                                                                      <  2:542468/150‑1 (MQ=255)
                                 tGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAg                                                                                                                      <  2:2935222/150‑1 (MQ=255)
                                 tGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAg                                                                                                                      <  2:219839/150‑1 (MQ=255)
                                  ggACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAgg                                                                                                                     <  2:4763685/150‑1 (MQ=255)
                                    aCAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGcc                                                                                                                   >  1:1743083/1‑150 (MQ=255)
                                     cAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCa                                                                                                                  >  1:3840578/1‑150 (MQ=255)
                                            aCGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAgctgct                                                                                                           <  1:4547874/150‑1 (MQ=255)
                                            aCGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAgctgct                                                                                                           >  1:1159464/1‑150 (MQ=255)
                                            aCGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAgctgct                                                                                                           <  2:4363846/150‑1 (MQ=255)
                                             cGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGctgctg                                                                                                          <  1:5121852/150‑1 (MQ=255)
                                                tCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGc                                                                                                       >  2:2267609/1‑150 (MQ=255)
                                                tCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGc                                                                                                       >  1:750221/1‑150 (MQ=255)
                                                      ctggAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAAtgtg                                                                                                 <  1:1032358/150‑1 (MQ=255)
                                                       tggAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAgg                                                                                                                     >  2:1108310/1‑129 (MQ=255)
                                                       tggAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAAtgtgt                                                                                                <  2:3828917/150‑1 (MQ=255)
                                                                       cccTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACggg                                                                                >  1:2043116/1‑150 (MQ=255)
                                                                        ccTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACATCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACgggg                                                                               >  1:6067485/1‑150 (MQ=255)
                                                                        ccTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACgggg                                                                               >  2:1499356/1‑150 (MQ=255)
                                                                              gTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTgg                                                                         <  1:1121432/150‑1 (MQ=255)
                                                                              gTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTgg                                                                         >  1:1374210/1‑150 (MQ=255)
                                                                              gTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTgg                                                                         <  2:5535962/150‑1 (MQ=255)
                                                                              gTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTgg                                                                         <  1:4719487/150‑1 (MQ=255)
                                                                              ctGCGCCACGTGAGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTgg                                                                         <  2:6272854/149‑1 (MQ=255)
                                                                               tGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTggg                                                                        >  1:2295578/1‑150 (MQ=255)
                                                                                  gcCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCa                                                                     >  2:5635894/1‑150 (MQ=255)
                                                                                   ccAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAg                                                                    >  2:1372411/1‑150 (MQ=255)
                                                                                    cAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGc                                                                   >  1:4887985/1‑150 (MQ=255)
                                                                                     aaGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCg                                                                  <  2:2350137/150‑1 (MQ=255)
                                                                                      aGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCgg                                                                 >  2:5139108/1‑150 (MQ=255)
                                                                                        tgtgCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGcc                                                               <  1:3953749/150‑1 (MQ=255)
                                                                                            ccAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCAttt                                                           <  1:3519784/150‑1 (MQ=255)
                                                                                             cAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTc                                                          >  1:2998396/1‑150 (MQ=255)
                                                                                               gTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCTTTTCAc                                                        <  2:5800523/150‑1 (MQ=255)
                                                                                               gTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCAc                                                        >  2:5348300/1‑150 (MQ=255)
                                                                                               gTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCAc                                                        >  1:572375/1‑150 (MQ=255)
                                                                                                 ggTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCt                                                      <  2:1772000/150‑1 (MQ=255)
                                                                                                   tCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTg                                                    <  1:3935272/150‑1 (MQ=255)
                                                                                                   tCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTg                                                    >  2:3681800/1‑150 (MQ=255)
                                                                                                   tCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTg                                                    >  1:6409038/1‑150 (MQ=255)
                                                                                                    cgcgTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTgg                                                   <  2:1330662/150‑1 (MQ=255)
                                                                                                    cgcgTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTgg                                                   <  2:320048/150‑1 (MQ=255)
                                                                                                         cGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGAc                                              <  1:3095760/150‑1 (MQ=255)
                                                                                                           gcgcCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACAt                                            >  1:2893210/1‑150 (MQ=255)
                                                                                                               cTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATgcgc                                        >  2:1914093/1‑150 (MQ=255)
                                                                                                                 gAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCAc                                      >  1:4877093/1‑150 (MQ=255)
                                                                                                                 gAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCAc                                      >  2:4020542/1‑150 (MQ=255)
                                                                                                                  aaaGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCAcc                                     >  2:6711598/1‑150 (MQ=255)
                                                                                                                    agcgggggACACCCACCCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCtt                                   <  2:5555729/147‑1 (MQ=255)
                                                                                                                     gagGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTg                                  <  2:4549809/150‑1 (MQ=255)
                                                                                                                       gggggACACCCAACCGGTGCAGGGCCGGCGAGCCCTGCTGGTCGGCCACACCCCGGGGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCAGTTCACCTTGGTCGACAGGCGCACCTTGca                                >  2:6829076/1‑150 (MQ=255)
                                                                                                                        ggggACACCTCACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGGGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTGCACCGTGGTCGCCATGCGCACCGTGcac                               >  2:4884345/1‑149 (MQ=255)
                                                                                                                        ggggACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGcaa                               <  2:6121133/150‑1 (MQ=255)
                                                                                                                        ggggACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGcaa                               <  2:6615330/150‑1 (MQ=255)
                                                                                                                         gggACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGcaac                              <  2:4871429/150‑1 (MQ=255)
                                                                                                                          ggACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGcaaca                             >  1:1413611/1‑150 (MQ=255)
                                                                                                                           gACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGcaacaa                            >  2:4659062/1‑150 (MQ=255)
                                                                                                                             cacCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGt                          <  1:5877043/150‑1 (MQ=255)
                                                                                                                              acCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTc                         >  1:6835306/1‑150 (MQ=255)
                                                                                                                               cccAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCa                        <  1:6202710/150‑1 (MQ=255)
                                                                                                                                 cAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGc                      >  2:6050045/1‑150 (MQ=255)
                                                                                                                                 cAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGc                      >  1:6539319/1‑150 (MQ=255)
                                                                                                                                  aaCCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGcc                     >  2:2158623/1‑150 (MQ=255)
                                                                                                                                  aaCCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGcc                     >  2:6851626/1‑150 (MQ=255)
                                                                                                                                   aCCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGccc                    >  1:1551685/1‑150 (MQ=255)
                                                                                                                                      ggTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCgg                 <  2:3532750/150‑1 (MQ=255)
                                                                                                                                      ccagCAGGGCTCGCAGGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCgg                 >  1:986202/4‑150 (MQ=255)
                                                                                                                                        tGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGc               >  1:4947485/1‑150 (MQ=255)
                                                                                                                                        tGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGc               <  1:2226019/150‑1 (MQ=255)
                                                                                                                                         gCAGGGCTCGCAGGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGcc              >  2:4420065/1‑150 (MQ=255)
                                                                                                                                         gCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGcc              <  1:5490805/150‑1 (MQ=255)
                                                                                                                                          cAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCa             <  1:824190/150‑1 (MQ=255)
                                                                                                                                           aGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAg            >  2:4580557/1‑150 (MQ=255)
                                                                                                                                            gggCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTgg                                                                         <  1:5526551/88‑1 (MQ=255)
                                                                                                                                             ggCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGAt          >  2:4540880/1‑150 (MQ=255)
                                                                                                                                                cTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATgca       <  1:5718528/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCCTGCTGCATGCAGACAGCCCGTTTGCCGACTGGACAGACTTACGCTCCTGGCTGGAACTGGATGACGACCCTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAAC  >  NC_002947/3412647‑3412945

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: