Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,554,503 G→C 6.0% E26Q (GAG→CAG)  PP_3139 → group 1 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,554,5030GC 6.0% 491.9 / 18.3 151E26Q (GAG→CAG) PP_3139group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (78/64);  new base C (5/4);  total (83/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

Only 100 of 159 total aligned reads displayed.

GCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGC  >  NC_002947/3554354‑3554644
                                                                                                                                                     |                                                                                                                                             
gCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTgg                                                                                                                                               <  2:2562754/150‑1 (MQ=255)
 cccGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGa                                                                                                                                              >  2:1830075/1‑150 (MQ=255)
  ccGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAg                                                                                                                                             <  2:2440815/150‑1 (MQ=255)
    gTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAgcg                                                                                                                                           >  1:711299/1‑150 (MQ=255)
    gTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAgcg                                                                                                                                           >  1:1479368/1‑150 (MQ=255)
        tCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCaa                                                                                                                                       <  1:1140644/150‑1 (MQ=255)
          gACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGATTGCGCGCCGAGTCAGTACTCGAAAGTCAGCCACCGTTTCGTACGCCACGCAATTTTGGCGCTGGAGCGCCAAgg                                                                                                                                     >  1:1594848/1‑150 (MQ=255)
           aCAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAggg                                                                                                                                    <  2:1396416/150‑1 (MQ=255)
           aCAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAggg                                                                                                                                    >  1:2099191/1‑150 (MQ=255)
               ccTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTg                                                                                                                                <  1:1488684/150‑1 (MQ=255)
               ccTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTg                                                                                                                                >  1:1879180/1‑150 (MQ=255)
                 tgaagtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAg                                                                                                                              >  2:2841558/1‑150 (MQ=255)
                   aagtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCTCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGt                                                                                                                            <  2:275986/150‑1 (MQ=255)
                    agtgaagGTTGCCAGGGGAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTg                                                                                                                           <  2:711299/150‑1 (MQ=255)
                     gtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGc                                                                                                                          >  2:2486313/1‑150 (MQ=255)
                     gtgaagGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCAATTTTTTATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGc                                                                                                                          <  1:561189/150‑1 (MQ=255)
                           gTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGta                                                                                                                    >  1:457280/1‑150 (MQ=255)
                              gCCATGGAAGGTTTGTCAACAATGGGATTTGGATTTGAATTCGTATGCGCATTGTTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAg                                                                                                                 <  2:603335/150‑1 (MQ=255)
                                cAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGcc                                                                                                               >  2:2084897/1‑150 (MQ=255)
                                cAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGcc                                                                                                               >  2:477011/1‑150 (MQ=255)
                                 aGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCaa                                                                                                                        >  1:461780/1‑140 (MQ=255)
                                  gggAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAg                                                                                                                                      >  1:2549411/1‑125 (MQ=255)
                                   ggAAGGTTGGTCATAAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTg                                                                                                            >  1:2409792/1‑150 (MQ=255)
                                        gTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgg                                                                                                       <  1:2878937/150‑1 (MQ=255)
                                           ggTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtg                                                                                                    <  2:457280/150‑1 (MQ=255)
                                            gTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAATGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtgg                                                                                                   >  2:1660957/1‑150 (MQ=255)
                                            gTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtgg                                                                                                   <  1:477011/150‑1 (MQ=255)
                                                tcaAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATTCTGGCGCAGGGGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACg                                                                                               <  1:2236585/150‑1 (MQ=255)
                                                    gCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAgg                                                                                                                             <  2:1592535/116‑1 (MQ=255)
                                                    gCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGa                                                                                           >  2:1455324/1‑150 (MQ=255)
                                                        gTTTGGATCTGGATTCGGATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTg                                                                                       >  1:2375520/1‑150 (MQ=255)
                                                          ttGGATCTGGATTCGTATGCGCATTTCTTACTTCATGAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGa                                                                                                                               <  2:525310/108‑1 (MQ=255)
                                                              aTCTGGATTCGTATGCGCATTGCTTACTTCGTCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGAt                                                                                 >  2:1487567/1‑150 (MQ=255)
                                                              aTCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGAt                                                                                 >  2:1523849/1‑150 (MQ=255)
                                                               tCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGt                                                                                        <  2:167254/142‑1 (MQ=255)
                                                               tCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGt                                                                                        >  1:167254/1‑142 (MQ=255)
                                                                 tGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGcc                                                                              >  1:1582361/1‑150 (MQ=255)
                                                                  ggATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGtt                                                                                                                                 >  2:1693132/1‑98 (MQ=255)
                                                                     ttCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAgg                                                                          >  2:404563/1‑150 (MQ=255)
                                                                           tGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCaa                                                                    >  1:977434/1‑150 (MQ=255)
                                                                           tGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCaa                                                                    <  2:977434/150‑1 (MQ=255)
                                                                            gcgcATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAg                                                                   <  1:404563/150‑1 (MQ=255)
                                                                                  tGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACg                                                             <  1:1157755/150‑1 (MQ=255)
                                                                                   gCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgg                                                            <  1:2921872/150‑1 (MQ=255)
                                                                                     ttACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgggg                                                          >  1:1297999/1‑150 (MQ=255)
                                                                                     ttACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgggg                                                          >  2:3004613/1‑150 (MQ=255)
                                                                                     ttACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgggg                                                          <  2:2409792/150‑1 (MQ=255)
                                                                                     ttACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACgggg                                                          <  1:1347767/150‑1 (MQ=255)
                                                                                       aCTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGt                                                                                        <  1:305139/118‑1 (MQ=255)
                                                                                       aCTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCa                                                        >  2:1136614/1‑150 (MQ=255)
                                                                                         ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAg                                                      >  2:2829278/1‑150 (MQ=255)
                                                                                         ttCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAg                                                      >  1:2563444/1‑150 (MQ=255)
                                                                                           catcaATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTTGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGAc                                                    <  1:163431/150‑1 (MQ=255)
                                                                                            atcaatcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCaa                                                                    <  2:2762307/133‑1 (MQ=255)
                                                                                            atcaatcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCaa                                                                    >  1:2762307/1‑133 (MQ=255)
                                                                                             tcaatcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGAccc                                                  >  2:361544/1‑150 (MQ=255)
                                                                                                atcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgg                                                                                                       <  2:449198/94‑1 (MQ=255)
                                                                                                atcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGa                                                                                                 >  1:695046/1‑100 (MQ=255)
                                                                                                atcaGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCt                                               >  1:1654483/1‑150 (MQ=255)
                                                                                                   aGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGCCGCAAGCGAACGGGGCAAAAACCGCGACg                                            >  2:17506/1‑150 (MQ=255)
                                                                                                    gTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtg                                                                                                    <  1:207347/93‑1 (MQ=255)
                                                                                                    gTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCgtgggtg                                                                                                    >  2:207347/1‑93 (MQ=255)
                                                                                                        ccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGAt                                                                                 <  2:475203/108‑1 (MQ=255)
                                                                                                        ccGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTg                                       >  2:2965731/1‑150 (MQ=255)
                                                                                                          gAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGAccc                                                  >  1:2797533/1‑137 (MQ=255)
                                                                                                               tCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTg                                >  2:1660691/1‑150 (MQ=255)
                                                                                                               tCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTCCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTg                                >  2:1996791/1‑150 (MQ=255)
                                                                                                                 aGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTggg                              <  2:2563444/150‑1 (MQ=255)
                                                                                                                 aGCCACAGTTTCATCCGTCGCCAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTggg                              <  2:1582361/150‑1 (MQ=255)
                                                                                                                  gCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGc                             >  1:1202738/1‑150 (MQ=255)
                                                                                                                   ccACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCt                            >  1:2138450/1‑150 (MQ=255)
                                                                                                                   ccACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCt                            >  1:143019/1‑150 (MQ=255)
                                                                                                                    cacaGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCt                                                                                                             >  2:2393798/1‑68 (MQ=255)
                                                                                                                    cacaGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCt                                                                                                             <  1:2393798/68‑1 (MQ=255)
                                                                                                                    cacaGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTc                           <  1:2084897/150‑1 (MQ=255)
                                                                                                                     acaGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCaa                                                                    >  1:1562494/1‑108 (MQ=255)
                                                                                                                       aGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAg                        <  2:1878841/150‑1 (MQ=255)
                                                                                                                       aGTTTCATCCGTCGCGAAAGACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAg                        >  2:3053038/1‑150 (MQ=255)
                                                                                                                        gTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAgg                       <  1:1830075/150‑1 (MQ=255)
                                                                                                                              tCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGt                                 <  1:1770437/134‑1 (MQ=255)
                                                                                                                              tCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGgttgt                 >  1:1708984/1‑150 (MQ=255)
                                                                                                                              tCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGgttgt                 <  1:920955/150‑1 (MQ=255)
                                                                                                                               ccGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAGGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGgttgtt                <  1:1660957/150‑1 (MQ=255)
                                                                                                                                  tCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTTCAAGGTGTGCTCAAGGGTGTGTTGGc             >  2:2193275/1‑150 (MQ=255)
                                                                                                                                  tCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGc                             <  2:2825018/134‑1 (MQ=255)
                                                                                                                                  tCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGc             >  1:1087851/1‑150 (MQ=255)
                                                                                                                                  tCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACATATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGgtt                   >  1:487454/1‑144 (MQ=255)
                                                                                                                                       aaaTACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGa        >  2:49030/1‑150 (MQ=255)
                                                                                                                                         cccttGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACg      >  1:2716467/5‑150 (MQ=255)
                                                                                                                                         cccttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCaa                                                                    <  1:431776/84‑1 (MQ=255)
                                                                                                                                          tACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGt     <  2:2138450/150‑1 (MQ=255)
                                                                                                                                          ccttGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGa                                                     >  2:1228779/4‑102 (MQ=255)
                                                                                                                                          ccttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTAGGGCAAGGTggg                              >  1:2462925/4‑125 (MQ=255)
                                                                                                                                           cttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCaa                                                       >  2:2135066/3‑99 (MQ=255)
                                                                                                                                           cttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCaa                                                       <  1:2135066/97‑1 (MQ=255)
                                                                                                                                           cttGGCGCTCCAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTg    >  1:986765/3‑150 (MQ=255)
                                                                                                                                                    ggAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTAc                                             >  2:238573/1‑100 (MQ=255)
                                                                                                                                                    ggAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGc  <  1:915970/143‑1 (MQ=255)
                                                                                                                                                    ggAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGc  >  2:915970/1‑143 (MQ=255)
                                                                                                                                                     gAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGc                             <  1:1998444/115‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GCCCGTAGTCGACAGCCTGAAGTGAAGGTTGCCAGGGAAGGTTGGTCATCAAGCGGGTTTGGATCTGGATTCGTATGCGCATTGCTTACTTCATCAATCAGTACCCGAAAGTCAGCCACAGTTTCATCCGTCGCGAAATACTGGCGCTGGAGCGCCAAGGGGTTGAGGTGCAACGTATAGCCCTGCGTGGGTGGGACGCCGAAGTGCAGGATGCCGAGGACGCAAGCGAACGGGGCAAGACCCGCTACGTATTGCAAGGTGGGCTCAAGGGGTTGTTGGCACCGACGTGGC  >  NC_002947/3554354‑3554644

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