Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,387,854 A→C I38L (ATC→CTC)  cynS → cyanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,387,8540AC99.8% 1855.8 / NA 478I38L (ATC→CTC) cynScyanase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (207/270);  minor base G (1/0);  total (208/270)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 485 total aligned reads displayed.

CCCTTCCCTCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTATCGCCCTGACATTAGTCACCGCCGCCATGTTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCATGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCG  >  NC_008463/3387706‑3388001
                                                                                                                                                    |                                                                                                                                                   
cccttccctCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGtctc                                                                                                                                                   <  1:4878531/151‑1 (MQ=255)
cccttccctCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGtctc                                                                                                                                                   <  2:2334602/151‑1 (MQ=255)
   gtccctCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACCGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGcc                                                                                                                                                <  1:3800310/150‑1 (MQ=255)
   gtccctCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCACCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGcc                                                                                                                                                <  1:8196474/150‑1 (MQ=255)
    tccctCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGccc                                                                                                                                               >  2:10100180/1‑151 (MQ=255)
     ccctCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCt                                                                                                                                              >  1:6048656/1‑151 (MQ=255)
       ctCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGa                                                                                                                                            >  1:7349401/1‑151 (MQ=255)
       ctCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGa                                                                                                                                            <  2:1881956/151‑1 (MQ=255)
        tCACAACCACCCGCAAGGAGAGTCATCCCTTGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGAc                                                                                                                                           >  1:11131866/1‑151 (MQ=255)
        tCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGAc                                                                                                                                           >  1:11221826/1‑151 (MQ=255)
         cacaACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGAcc                                                                                                                                          >  1:6734764/1‑150 (MQ=255)
            aaCCCCCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCttc                                                                                                                                       <  2:207253/151‑2 (MQ=255)
            aaCCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCttc                                                                                                                                       <  2:3168240/151‑2 (MQ=255)
               caccCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTc                                                                                                                                    <  2:8670895/151‑1 (MQ=255)
                  ccGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAcc                                                                                                                                 >  2:6544271/1‑151 (MQ=255)
                  ccGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAcc                                                                                                                                 >  2:7228417/1‑151 (MQ=255)
                   cGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccg                                                                                                                                <  2:5982194/151‑1 (MQ=255)
                    gCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccgc                                                                                                                               >  1:7606503/1‑151 (MQ=255)
                     cAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccgcc                                                                                                                              >  1:6888620/1‑151 (MQ=255)
                     cAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccgcc                                                                                                                              >  2:10605334/1‑151 (MQ=255)
                      aaGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccgccg                                                                                                                             <  2:3769486/151‑1 (MQ=255)
                      aaGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccgccg                                                                                                                             <  1:3769334/151‑1 (MQ=255)
                       aGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCGGTCCTGGCCGGATCTGGCCGACGGCACCGGGCTCGCCCTGACCTTCGTCACcgccgc                                                                                                                            >  2:6668917/1‑151 (MQ=255)
                         gagaGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCgccgccc                                                                                                                          <  1:7914886/151‑2 (MQ=255)
                           gagTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTg                                                                                                                        >  1:3548861/1‑151 (MQ=255)
                            agTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTgc                                                                                                                       <  1:1083402/151‑2 (MQ=255)
                             gTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTgct                                                                                                                      <  1:1913662/151‑3 (MQ=255)
                              tCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTc                                                                                                                     <  2:1633290/151‑1 (MQ=255)
                                   ccATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCa                                                                                                                <  2:488382/151‑1 (MQ=255)
                                   ccATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCa                                                                                                                <  2:3377602/151‑1 (MQ=255)
                                   ccATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCa                                                                                                                <  1:9207568/151‑1 (MQ=255)
                                    cATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAg                                                                                                               <  2:274396/151‑1 (MQ=255)
                                     aTGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGc                                                                                                              <  1:1148244/151‑1 (MQ=255)
                                        cAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATg                                                                                                           <  2:11187830/151‑1 (MQ=255)
                                         aaCACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACAGCCGCCCTGCTCGGCCAGCATGc                                                                                                          >  2:10997503/1‑151 (MQ=255)
                                           cacaGCCAAGTGTCCCCCAACGCGCGCCAGCAACAGGCCGAGACCGAGGTCATGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGccc                                                                                                        >  1:8534068/1‑151 (MQ=255)
                                             ccgCCAAGTGTCCCCCACCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTg                                                                                                      <  2:10314196/149‑1 (MQ=255)
                                             caGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTg                                                                                                      <  2:3351777/151‑1 (MQ=255)
                                              aGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGc                                                                                                     >  2:7352123/1‑151 (MQ=255)
                                              aGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGc                                                                                                     >  1:3954371/1‑151 (MQ=255)
                                                 cAAGTGTCCCCCATCGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCag                                                                                                  <  1:5672465/151‑1 (MQ=255)
                                                 cAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGGCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCag                                                                                                  >  2:2282429/1‑151 (MQ=255)
                                                   aGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCCGCATGCCCTGCCagag                                                                                                >  2:2556597/1‑151 (MQ=255)
                                                    gtgtCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAgg                                                                                               <  2:119075/151‑1 (MQ=255)
                                                      gtCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcg                                                                                             >  2:7146876/1‑151 (MQ=255)
                                                       tCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcgg                                                                                            <  1:1391250/151‑1 (MQ=255)
                                                           ccAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCCGCAGGCCCTGCCAGAggcggcggc                                                                                        >  2:7041532/1‑151 (MQ=255)
                                                            cAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGgcggcggcg                                                                                       >  2:5413900/1‑151 (MQ=255)
                                                            cAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGgcggcggcg                                                                                       >  1:2745994/1‑151 (MQ=255)
                                                                 cgcgcCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCaa                                                                                  <  2:295722/151‑1 (MQ=255)
                                                                   cgcCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAgg                                                                                >  2:10352732/1‑151 (MQ=255)
                                                                    gcCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCGTGCTCGGCCAGCTTGCCCTGCCAGTGGCGGCGGCGCGCTAGGt                                                                               >  2:9482994/1‑151 (MQ=255)
                                                                     ccAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTc                                                                              >  1:6545053/1‑151 (MQ=255)
                                                                      cAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcg                                                                             <  2:10547984/151‑1 (MQ=255)
                                                                         cAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgccg                                                                          >  1:7147528/1‑151 (MQ=255)
                                                                         cAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgccg                                                                          <  2:8581351/151‑1 (MQ=255)
                                                                          aaCTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCAGGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgccgc                                                                         <  2:4296582/151‑1 (MQ=255)
                                                                           aCTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCgccgcc                                                                        >  1:8793329/1‑151 (MQ=255)
                                                                            cTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGCCGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGccc                                                                       <  1:5265218/151‑1 (MQ=255)
                                                                            cTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGccc                                                                       <  2:7811785/151‑1 (MQ=255)
                                                                               gCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAAc                                                                    >  2:3522004/1‑151 (MQ=255)
                                                                                   agaCCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCgg                                                                >  2:8397759/1‑151 (MQ=255)
                                                                                    gaCCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGc                                                               >  2:1932361/1‑151 (MQ=255)
                                                                                     aCCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGcc                                                              <  2:10806759/151‑2 (MQ=255)
                                                                                          ggTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgac                                                         >  2:7733004/1‑151 (MQ=255)
                                                                                           gTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacg                                                        <  2:8056370/151‑1 (MQ=255)
                                                                                               tGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGACGCCGCCCAACACGGCCTGgacgacga                                                    <  1:9908983/151‑1 (MQ=255)
                                                                                                gaccAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacgacgac                                                   <  1:2375852/148‑1 (MQ=255)
                                                                                                gAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacgacgac                                                   <  2:1489542/151‑1 (MQ=255)
                                                                                                gAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacgacgac                                                   <  2:6241381/151‑1 (MQ=255)
                                                                                                  aCAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCCTGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGAcgc                                                 <  2:5674132/151‑1 (MQ=255)
                                                                                                  aCAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGAcgc                                                 >  1:10024784/1‑151 (MQ=255)
                                                                                                  aCAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGAcgc                                                 <  1:9067789/151‑1 (MQ=255)
                                                                                                   cAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGAcgcg                                                >  2:11394662/1‑151 (MQ=255)
                                                                                                    aaGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCgg                                               <  1:3350486/151‑1 (MQ=255)
                                                                                                     aGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGt                                              >  1:9795513/1‑151 (MQ=255)
                                                                                                      ggCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTg                                             >  2:6390991/1‑151 (MQ=255)
                                                                                                       gCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGc                                            <  2:5226512/151‑1 (MQ=255)
                                                                                                          cGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTcc                                         >  2:7268777/1‑151 (MQ=255)
                                                                                                              tGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgct                                     <  2:4586243/151‑1 (MQ=255)
                                                                                                               gggCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgctg                                    >  1:4605638/1‑151 (MQ=255)
                                                                                                                     gTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGc                              <  2:9705329/151‑1 (MQ=255)
                                                                                                                     gTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGc                              <  2:9266821/151‑1 (MQ=255)
                                                                                                                     gTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGc                              <  1:9502367/151‑1 (MQ=255)
                                                                                                                          ggCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATccc                         <  1:10511999/151‑1 (MQ=255)
                                                                                                                             aGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCt                      <  1:10220551/151‑1 (MQ=255)
                                                                                                                             aGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCt                      <  1:4964873/151‑1 (MQ=255)
                                                                                                                                  cTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGcgcg                 <  2:8919141/151‑1 (MQ=255)
                                                                                                                                   tGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCgg                <  2:9022091/151‑1 (MQ=255)
                                                                                                                                    ggCCGCCGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCCGAGGCGGCGGCGCGCAAGGTCTCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgc               <  2:1377332/151‑1 (MQ=255)
                                                                                                                                     gCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgca              <  1:6155254/151‑1 (MQ=255)
                                                                                                                                       cGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcagc            >  1:10712062/1‑151 (MQ=255)
                                                                                                                                       cGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcagc            >  1:10350964/1‑151 (MQ=255)
                                                                                                                                       cGACGGCACCGGTCTCGCCCTGACCGTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcagc            >  2:8170379/1‑151 (MQ=255)
                                                                                                                                         aCGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCAt          <  1:3008374/151‑1 (MQ=255)
                                                                                                                                             cACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTccc      >  1:8884521/1‑151 (MQ=255)
                                                                                                                                             cACCGGACTCGCCCAGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCAAGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAACGCATCCCGCTGCGCGGCAGCATTccc      <  1:4204433/151‑1 (MQ=255)
                                                                                                                                                cGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCcggc   <  1:1710990/151‑1 (MQ=255)
                                                                                                                                                cGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCcggc   >  2:4792383/1‑151 (MQ=255)
                                                                                                                                                 ggTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCcggcg  >  1:6378520/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCTTCCCTCACAACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTATCGCCCTGACATTAGTCACCGCCGCCATGTTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCATGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCG  >  NC_008463/3387706‑3388001

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: