Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,387,881 A→C M47L (ATG→CTG)  cynS → cyanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,387,8810AC99.8% 1754.3 / NA 450M47L (ATG→CTG) cynScyanase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (193/256);  total (193/257)
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 450 total aligned reads displayed.

AGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTATCGCCCTGACATTAGTCACCGCCGCCATGTTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCATGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGCGACCCGACCATCTAC  >  NC_008463/3387734‑3388027
                                                                                                                                                   |                                                                                                                                                  
agTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCGTCGTCACCGCCGCCCTg                                                                                                                                                  >  1:5690880/1‑150 (MQ=255)
 gTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTgc                                                                                                                                                 <  2:2291810/150‑2 (MQ=255)
   cATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTc                                                                                                                                               <  2:5459869/150‑1 (MQ=255)
   cATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTc                                                                                                                                               <  1:2036638/150‑1 (MQ=255)
    aTCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCg                                                                                                                                              >  2:10574889/1‑150 (MQ=255)
       ccATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGcc                                                                                                                                           <  1:6798374/150‑1 (MQ=255)
         aTGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAg                                                                                                                                         <  2:9627298/150‑1 (MQ=255)
          tGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGc                                                                                                                                        >  2:280108/1‑150 (MQ=255)
           gCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCa                                                                                                                                       >  2:4556142/1‑150 (MQ=255)
            cAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCAt                                                                                                                                      >  2:2508435/1‑150 (MQ=255)
             aaCACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATg                                                                                                                                     >  2:7084342/1‑150 (MQ=255)
              acacAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGc                                                                                                                                    >  1:7753843/1‑150 (MQ=255)
                 caGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCt                                                                                                                                 >  2:7189906/1‑150 (MQ=255)
                     cAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCAGGCCCTGCCa                                                                                                                             <  1:534032/150‑1 (MQ=255)
                       aGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCaga                                                                                                                           >  1:10341449/1‑150 (MQ=255)
                       aGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCaga                                                                                                                           >  2:5832970/1‑150 (MQ=255)
                        gtgtCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCagag                                                                                                                          >  2:9208147/1‑150 (MQ=255)
                         tgtCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAgg                                                                                                                         >  1:6404221/1‑150 (MQ=255)
                           tCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcg                                                                                                                       <  1:1877290/150‑1 (MQ=255)
                           tCCACCAACGCGCGCCTGCAACTGGCCGGGACCGTGGTGCTGCACATGGCCCGCCTGGGCCTGTCCTGGCCGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTGCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcg                                                                                                                       >  2:1163701/1‑150 (MQ=255)
                               ccAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcgg                                                                                                                   <  2:2613307/150‑1 (MQ=255)
                                cAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcggc                                                                                                                  >  2:6733754/1‑150 (MQ=255)
                                cAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcggc                                                                                                                  >  2:7834907/1‑150 (MQ=255)
                                 aaCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGgcggcggcg                                                                                                                 >  1:2327796/1‑150 (MQ=255)
                                 aaCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGgcggcggcg                                                                                                                 >  2:3960416/1‑150 (MQ=255)
                                  aCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGgcgc                                                                                                                >  1:3629957/1‑150 (MQ=255)
                                     cgcgcCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCa                                                                                                             <  2:4200907/150‑1 (MQ=255)
                                     cgcgcCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCa                                                                                                             >  2:3497619/1‑150 (MQ=255)
                                        gcCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAgg                                                                                                          <  1:3236474/150‑1 (MQ=255)
                                         ccAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGt                                                                                                         >  1:9617913/1‑150 (MQ=255)
                                          cAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTc                                                                                                        >  1:1121240/1‑150 (MQ=255)
                                            gCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgc                                                                                                      <  1:5220753/150‑1 (MQ=255)
                                             cAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgcc                                                                                                     >  2:142064/1‑150 (MQ=255)
                                               aCTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgccgc                                                                                                   >  1:1649910/1‑150 (MQ=255)
                                                 tGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGccc                                                                                                 <  2:7507346/150‑1 (MQ=255)
                                                    ccGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAAc                                                                                              <  2:8790807/150‑1 (MQ=255)
                                                      gagaCCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTc                                                                                            <  1:10035389/150‑1 (MQ=255)
                                                          ccGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGcc                                                                                        <  2:32983/150‑2 (MQ=255)
                                                           cGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGcct                                                                                       <  1:81934/150‑3 (MQ=255)
                                                           cGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGcct                                                                                       >  1:8965088/1‑148 (MQ=255)
                                                            gTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTg                                                                                      <  2:8748061/150‑1 (MQ=255)
                                                              ggTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGga                                                                                    <  1:4057700/150‑1 (MQ=255)
                                                                 ccTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacga                                                                                 <  2:4929631/150‑1 (MQ=255)
                                                                 ccTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacga                                                                                 <  1:3835667/150‑1 (MQ=255)
                                                                 ccTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacga                                                                                 >  1:4181454/1‑150 (MQ=255)
                                                                      aCAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGacgacgacg                                                                            <  1:9103350/150‑1 (MQ=255)
                                                                      aCAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGACCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCTCAAGGTCGCCGCCCAACTCGGACTGGacgacgacg                                                                            >  2:8967111/1‑150 (MQ=255)
                                                                         aGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCgg                                                                         >  2:4723387/1‑150 (MQ=255)
                                                                         aGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCgg                                                                         >  2:7646144/1‑150 (MQ=255)
                                                                          ggCCCGCCTGGGCCTGTCCGGGCAGGTTCGGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACTCGGt                                                                        <  2:9554974/150‑1 (MQ=255)
                                                                                  tGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgc                                                                <  1:10961062/150‑1 (MQ=255)
                                                                                   gggCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgct                                                               <  2:9279527/150‑1 (MQ=255)
                                                                                   gggCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgct                                                               >  2:10471125/1‑150 (MQ=255)
                                                                                   gggCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgct                                                               <  2:1812020/150‑1 (MQ=255)
                                                                                     gCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGTCGACGACGCGGTGCTCCtgctgc                                                             >  2:202099/1‑150 (MQ=255)
                                                                                      ccTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCa                                                            <  2:2448966/150‑1 (MQ=255)
                                                                                       cTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCaa                                                           <  1:7006196/150‑1 (MQ=255)
                                                                                         gTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAg                                                         >  2:8480908/1‑150 (MQ=255)
                                                                                            cTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCAt                                                      <  1:3519158/150‑1 (MQ=255)
                                                                                              ggCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATcc                                                    <  2:9391451/150‑1 (MQ=255)
                                                                                                cAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCg                                                  <  1:8553610/150‑1 (MQ=255)
                                                                                                    aTCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTgc                                              <  2:5537611/150‑1 (MQ=255)
                                                                                                        ggCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCACGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCgg                                          <  1:432137/150‑1 (MQ=255)
                                                                                                         gCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgc                                         <  1:6189504/150‑1 (MQ=255)
                                                                                                         gCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgc                                         <  1:1412168/150‑1 (MQ=255)
                                                                                                         gCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgc                                         <  2:6627874/150‑1 (MQ=255)
                                                                                                          ccGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgca                                        >  1:10667056/1‑150 (MQ=255)
                                                                                                          ccGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgca                                        <  2:1303874/150‑1 (MQ=255)
                                                                                                           cGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcag                                       <  1:4815724/150‑1 (MQ=255)
                                                                                                           cGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGACGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcag                                       <  2:1093256/150‑1 (MQ=255)
                                                                                                             tcGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcagca                                     <  1:3098049/149‑1 (MQ=255)
                                                                                                             aCGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcagca                                     <  1:4586323/150‑1 (MQ=255)
                                                                                                               ggCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCACCTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcagcagt                                   >  1:9370989/1‑148 (MQ=255)
                                                                                                               ggCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCAtt                                   >  1:8833707/1‑150 (MQ=255)
                                                                                                                gCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTc                                  >  1:4405357/1‑150 (MQ=255)
                                                                                                                 cACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTcc                                 <  2:2954929/150‑1 (MQ=255)
                                                                                                                   ccGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCcg                               >  1:9798616/1‑150 (MQ=255)
                                                                                                                    cGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCcgg                              >  1:5736416/1‑150 (MQ=255)
                                                                                                                      gtctCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCcggcg                            >  2:8096815/4‑150 (MQ=255)
                                                                                                                        ctCGCCCAGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCggcggc                          >  1:7658818/2‑150 (MQ=255)
                                                                                                                         tCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCa                         <  1:7991982/150‑1 (MQ=255)
                                                                                                                          cGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCAt                        >  1:2890147/1‑150 (MQ=255)
                                                                                                                          cGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCAt                        <  2:10847994/150‑1 (MQ=255)
                                                                                                                             ccTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTcc                     <  1:8411185/150‑1 (MQ=255)
                                                                                                                             ccTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTcc                     <  1:6060048/150‑1 (MQ=255)
                                                                                                                              cTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCCGCATTCCCGGCGGCAGTccc                    >  2:8338109/1‑150 (MQ=255)
                                                                                                                              cTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTccc                    <  1:9934236/150‑1 (MQ=255)
                                                                                                                              cTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTccc                    >  1:10496020/1‑150 (MQ=255)
                                                                                                                                gaccgtcgTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAg                  <  2:10168018/143‑1 (MQ=255)
                                                                                                                                 accttcgTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGc                 <  1:5915264/144‑1 (MQ=255)
                                                                                                                                   cttcgTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGcga               <  2:3769062/146‑1 (MQ=255)
                                                                                                                                    ttcgTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGcgac              <  1:2597877/147‑1 (MQ=255)
                                                                                                                                      cgTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGcgaccc            <  2:6122193/149‑1 (MQ=255)
                                                                                                                                      cgTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGcgaccc            <  1:7310082/149‑1 (MQ=255)
                                                                                                                                        tCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGcgacccga          <  1:3335456/150‑1 (MQ=255)
                                                                                                                                         cACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGcgacccgac         <  1:2713215/150‑1 (MQ=255)
                                                                                                                                          aCCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGcgacccgacc        <  1:7040871/150‑1 (MQ=255)
                                                                                                                                           ccgccgCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGCGACCCGACCa       >  2:6356414/1‑150 (MQ=255)
                                                                                                                                             gccgccCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGCGACCCGACCATc     <  2:10405275/150‑1 (MQ=255)
                                                                                                                                              ccgccCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGCGACCCGACCATCt    >  2:6855546/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTATCGCCCTGACATTAGTCACCGCCGCCATGTTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCATGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAGCGACCCGACCATCTAC  >  NC_008463/3387734‑3388027

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: