Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,387,865 A→C T41T (ACA→ACC cynS → cyanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,387,8650AC99.8% 1797.7 / NA 462T41T (ACA→ACCcynScyanase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (204/257);  total (204/258)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 462 total aligned reads displayed.

AACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTATCGCCCTGACATTAGTCACCGCCGCCATGTTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCATGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAG  >  NC_008463/3387718‑3388011
                                                                                                                                                   |                                                                                                                                                  
aaCCACCCGCAAGGAGAGTCATCCCAGGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCtt                                                                                                                                                  <  2:10642883/150‑1 (MQ=255)
 accaccCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCACCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCttc                                                                                                                                                 <  2:3321733/150‑2 (MQ=255)
   caccCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGt                                                                                                                                               <  2:6806521/150‑1 (MQ=255)
   caccCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGt                                                                                                                                               >  2:4484217/1‑150 (MQ=255)
      ccGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAc                                                                                                                                            <  1:6027199/150‑1 (MQ=255)
        gCACGGAGAGTCATCCCATGCAACCCAGCCACGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccg                                                                                                                                          <  2:2496373/150‑1 (MQ=255)
        gCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccg                                                                                                                                          <  1:1000842/150‑1 (MQ=255)
         cAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccgc                                                                                                                                         <  1:9023691/150‑1 (MQ=255)
          aaGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccgcc                                                                                                                                        >  2:8987191/1‑150 (MQ=255)
            ggAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACcgccgc                                                                                                                                      <  1:207733/150‑1 (MQ=255)
               gagTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCgccgccct                                                                                                                                   >  1:9205658/1‑148 (MQ=255)
                agTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTg                                                                                                                                  >  2:4574175/1‑150 (MQ=255)
                agTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTg                                                                                                                                  <  1:3254248/150‑1 (MQ=255)
                   cATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTc                                                                                                                               >  2:5047388/1‑150 (MQ=255)
                    aTCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCg                                                                                                                              >  1:5733304/1‑150 (MQ=255)
                    aTCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCg                                                                                                                              >  2:2555422/1‑150 (MQ=255)
                    aTCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCg                                                                                                                              >  2:7274820/1‑150 (MQ=255)
                     tCCCATGCAACACAGCCAAGTGTCCCCCACCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCgg                                                                                                                             <  2:2312324/150‑1 (MQ=255)
                     tCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCgg                                                                                                                             <  2:10529910/150‑1 (MQ=255)
                      cccATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGc                                                                                                                            >  2:552950/1‑150 (MQ=255)
                       ccATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGcc                                                                                                                           >  2:1788565/1‑150 (MQ=255)
                        cATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCa                                                                                                                          >  1:1028712/1‑150 (MQ=255)
                         aTGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAg                                                                                                                         <  2:9627298/150‑1 (MQ=255)
                           gCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCa                                                                                                                       <  1:4297419/150‑1 (MQ=255)
                           gCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCa                                                                                                                       <  2:1984107/150‑1 (MQ=255)
                                 caGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCt                                                                                                                 <  1:2547201/150‑1 (MQ=255)
                                   gCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGc                                                                                                               <  1:433960/150‑1 (MQ=255)
                                     cAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCa                                                                                                             <  1:6106593/150‑1 (MQ=255)
                                     cAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCAGGCCCTGCCa                                                                                                             <  1:534032/150‑1 (MQ=255)
                                        gtgtCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCagag                                                                                                          >  2:9208147/1‑150 (MQ=255)
                                         tgtCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAgg                                                                                                         >  1:5654114/1‑150 (MQ=255)
                                         tgtCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAgg                                                                                                         <  1:6271192/150‑1 (MQ=255)
                                         tgtCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAgg                                                                                                         <  2:8411952/150‑1 (MQ=255)
                                          gtCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggc                                                                                                        <  1:4464091/150‑1 (MQ=255)
                                           tCCACCAACGCGCGCCTGCAACTGGCCGGGACCGTGGTGCTGCACATGGCCCGCCTGGGCCTGTCCTGGCCGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTGCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcg                                                                                                       >  2:1163701/1‑150 (MQ=255)
                                             ccccAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggc                                                                                                     <  2:7119270/150‑1 (MQ=255)
                                              cccACCGCGCGCCAGCACCTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcg                                                                                                    <  2:10162322/150‑1 (MQ=255)
                                              cccAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcg                                                                                                    >  1:4587570/1‑150 (MQ=255)
                                                cAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcggc                                                                                                  >  2:6733754/1‑150 (MQ=255)
                                                cAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcggc                                                                                                  >  2:7834907/1‑150 (MQ=255)
                                                cAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcggc                                                                                                  >  2:1425657/1‑150 (MQ=255)
                                                 aaCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACTGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGgcggcggcg                                                                                                 >  1:2325991/1‑150 (MQ=255)
                                                 aaCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGgcggcggcg                                                                                                 >  2:3960416/1‑150 (MQ=255)
                                                 aaCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGgcggcggcg                                                                                                 >  2:8818497/1‑150 (MQ=255)
                                                 aaCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGgcggcggcg                                                                                                 <  1:3899674/150‑1 (MQ=255)
                                                  aCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGgcgc                                                                                                <  1:8313436/150‑1 (MQ=255)
                                                  aCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGgcgc                                                                                                <  2:2816586/150‑1 (MQ=255)
                                                   cgcgcgCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGgcgcg                                                                                               >  1:10464574/1‑150 (MQ=255)
                                                       cgcCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAg                                                                                           <  2:58014/150‑1 (MQ=255)
                                                          cAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTc                                                                                        <  2:8270537/150‑1 (MQ=255)
                                                           aGCACCTGGCCGAGACCTGGGTCCAGAACAAGGCCCGCCTGGTCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcg                                                                                       <  2:10375690/150‑1 (MQ=255)
                                                           aGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcg                                                                                       <  2:3694136/150‑1 (MQ=255)
                                                             cAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgcc                                                                                     <  1:3581831/150‑1 (MQ=255)
                                                              aaCTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgccg                                                                                    >  1:5909131/1‑150 (MQ=255)
                                                               aCTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgccgc                                                                                   >  1:1649910/1‑150 (MQ=255)
                                                                 tGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGccc                                                                                 <  2:7507346/150‑1 (MQ=255)
                                                                  ggCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCa                                                                                <  2:2253118/150‑1 (MQ=255)
                                                                     cGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACt                                                                             >  2:770807/1‑150 (MQ=255)
                                                                           cGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGcct                                                                       >  1:9141804/1‑148 (MQ=255)
                                                                           cGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGcct                                                                       <  1:81934/150‑3 (MQ=255)
                                                                            ggggTCCTGAACAAGGCCCGCCGGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTg                                                                      <  2:8763846/148‑1 (MQ=255)
                                                                             tGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTgg                                                                     >  1:705484/1‑150 (MQ=255)
                                                                              ggTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGga                                                                    <  1:6616740/150‑1 (MQ=255)
                                                                              ggTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGga                                                                    <  1:2818894/150‑1 (MQ=255)
                                                                                tCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacg                                                                  >  1:7724286/1‑150 (MQ=255)
                                                                                 ccTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacga                                                                 >  2:263112/1‑150 (MQ=255)
                                                                                  cTGCACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacgac                                                                <  2:6257155/150‑1 (MQ=255)
                                                                                    gAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacgacga                                                              <  1:3046792/150‑1 (MQ=255)
                                                                                    gAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacgacga                                                              <  1:7375940/150‑1 (MQ=255)
                                                                                       cAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGAcgc                                                           <  2:3194165/150‑1 (MQ=255)
                                                                                         aGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCgg                                                         <  1:516163/150‑1 (MQ=255)
                                                                                         aGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCgg                                                         >  2:4723387/1‑150 (MQ=255)
                                                                                          ggCCCGCCTGGGCCTGTCCGGGCAGGTTCGGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACTCGGt                                                        <  2:9554974/150‑1 (MQ=255)
                                                                                                   gggCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgct                                               >  1:8882272/1‑150 (MQ=255)
                                                                                                   gggCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgct                                               <  1:1253744/150‑1 (MQ=255)
                                                                                                   gggCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgct                                               <  2:9279527/150‑1 (MQ=255)
                                                                                                      ccTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCa                                            <  2:2448966/150‑1 (MQ=255)
                                                                                                       cTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCaa                                           <  1:3790256/150‑1 (MQ=255)
                                                                                                       cTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCaa                                           <  1:9624288/150‑1 (MQ=255)
                                                                                                       cTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCaa                                           <  1:7006196/150‑1 (MQ=255)
                                                                                                         gTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAg                                         <  2:4199231/150‑1 (MQ=255)
                                                                                                                cAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCg                                  <  1:8553610/150‑1 (MQ=255)
                                                                                                                     tCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTgcg                             >  1:7782548/1‑150 (MQ=255)
                                                                                                                         gCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgc                         <  1:1412168/150‑1 (MQ=255)
                                                                                                                         gCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgc                         <  2:6627874/150‑1 (MQ=255)
                                                                                                                         gCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgc                         <  1:8832857/150‑1 (MQ=255)
                                                                                                                         gCCGACGGCACCGGTCTCGCCCTGACCTGCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgc                         <  2:8105353/150‑1 (MQ=255)
                                                                                                                          ccGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgca                        <  2:1303874/150‑1 (MQ=255)
                                                                                                                           cGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcag                       <  1:10544837/150‑1 (MQ=255)
                                                                                                                           cGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcag                       >  1:9593376/1‑150 (MQ=255)
                                                                                                                           cGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGACGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcag                       <  2:1093256/150‑1 (MQ=255)
                                                                                                                             aCGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcagca                     >  1:322154/1‑150 (MQ=255)
                                                                                                                                  aCCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTccc                >  1:9850449/1‑150 (MQ=255)
                                                                                                                                  aCCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTccc                <  1:4179634/150‑1 (MQ=255)
                                                                                                                                       tctCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCcggcgg           >  2:8555579/3‑150 (MQ=255)
                                                                                                                                          cGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCAt        <  2:10847994/150‑1 (MQ=255)
                                                                                                                                            cccTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTc      <  1:1114747/150‑1 (MQ=255)
                                                                                                                                              cTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTccc    >  1:1807097/1‑150 (MQ=255)
                                                                                                                                              cTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTccc    <  2:9436207/150‑1 (MQ=255)
                                                                                                                                               tGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCa   <  1:3396041/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTATCGCCCTGACATTAGTCACCGCCGCCATGTTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCATGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAG  >  NC_008463/3387718‑3388011

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: