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% 1642.6 / NA 424T41T (ACA→ACCcynScyanase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (197/226);  total (198/226)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 424 total aligned reads displayed.

AACCACCCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTATCGCCCTGACATTAGTCACCGCCGCCATGTTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCATGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAG  >  NC_008463/3387718‑3388011
                                                                                                                                                   |                                                                                                                                                  
 accaccCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCttc                                                                                                                                                 >  2:7975585/1‑149 (MQ=255)
 accaccCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCttc                                                                                                                                                 <  2:9325009/150‑2 (MQ=255)
  ccaccCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCGGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCttcg                                                                                                                                                <  1:3861117/150‑3 (MQ=255)
   caccCGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCCACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGt                                                                                                                                               <  2:3031215/150‑1 (MQ=255)
      ccGCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAATGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCACTGACCTTCGTCAc                                                                                                                                            >  1:10099089/1‑150 (MQ=255)
        gCAAGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccg                                                                                                                                          >  2:9622912/1‑150 (MQ=255)
          aaGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccgcc                                                                                                                                        <  1:7631893/150‑1 (MQ=255)
          aaGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccgcc                                                                                                                                        <  1:8342580/150‑1 (MQ=255)
           aGGAGAGTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCAccgccg                                                                                                                                       >  1:6406499/1‑150 (MQ=255)
              cgagTCATCCCATGCAACACAGCCACGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCgccgccc                                                                                                                                    <  2:2765366/149‑2 (MQ=255)
              agagTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCgccgccc                                                                                                                                    <  1:4612406/150‑2 (MQ=255)
               gagTCATCCCATGCAACACAGCCACGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCgccgccct                                                                                                                                   <  1:746463/150‑3 (MQ=255)
                 gTCATCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTgc                                                                                                                                 <  1:68341/150‑2 (MQ=255)
                     tCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCgg                                                                                                                             >  2:9204413/1‑150 (MQ=255)
                     tCCCATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCgg                                                                                                                             <  1:8563821/150‑1 (MQ=255)
                       ccATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGcc                                                                                                                           >  1:4159679/1‑150 (MQ=255)
                       ccATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGcc                                                                                                                           >  2:8885291/1‑150 (MQ=255)
                        cATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCa                                                                                                                          >  2:4381164/1‑150 (MQ=255)
                        cATGCAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCa                                                                                                                          <  1:1479417/150‑1 (MQ=255)
                         aTGCAACACAGCCACGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAg                                                                                                                         >  2:9700447/1‑150 (MQ=255)
                            cAACCCCGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCAt                                                                                                                      <  2:1341556/150‑1 (MQ=255)
                            cAACACAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCAt                                                                                                                      >  2:8417390/1‑150 (MQ=255)
                              acacAGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGc                                                                                                                    >  1:2504037/1‑150 (MQ=255)
                               cacaGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGcc                                                                                                                   <  2:8385014/150‑1 (MQ=255)
                                 caGCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCt                                                                                                                 <  2:2306969/150‑1 (MQ=255)
                                   gCCAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGTCCGTGGGCCTTAGCAAGGCCCGCCTGTGCCTGGCCTGGCCGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTGCGTCACCGCCGCCCTGCTCGGCCGGCTTGCCCGGc                                                                                                               >  2:5008645/1‑150 (MQ=255)
                                    ccAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGcc                                                                                                              >  2:9552131/1‑150 (MQ=255)
                                     cAAGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCa                                                                                                             <  1:3647745/150‑1 (MQ=255)
                                      aaGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCag                                                                                                            <  1:4337408/150‑1 (MQ=255)
                                      aaGTGTCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGGCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCGTCGGCACCGCCGCCCTGCTCGGCCAGCATGCCCTGcccg                                                                                                            >  2:10234759/1‑148 (MQ=255)
                                        gtgtCCCCCCACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCagag                                                                                                          <  2:9794595/150‑1 (MQ=255)
                                          gtCCCCCAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggc                                                                                                        <  1:11142968/150‑1 (MQ=255)
                                            cccccAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcgg                                                                                                      >  1:10008705/1‑150 (MQ=255)
                                              cccAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACACGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcg                                                                                                    <  2:6337586/150‑1 (MQ=255)
                                               ccAACGCGCGCCAGCACCTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcgg                                                                                                   <  2:5111231/150‑1 (MQ=255)
                                               ccAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcgg                                                                                                   >  2:3975924/1‑150 (MQ=255)
                                               ccAACGCGCGCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAggcggcgg                                                                                                   >  2:4071724/1‑150 (MQ=255)
                                                   cgcgcgCCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGgcgcg                                                                                               >  2:5738930/1‑150 (MQ=255)
                                                    gcgcgcCAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGgcgcgc                                                                                              >  1:3357341/1‑150 (MQ=255)
                                                         ccAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGt                                                                                         <  2:4519726/150‑1 (MQ=255)
                                                          cAGCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTc                                                                                        >  1:8062627/1‑150 (MQ=255)
                                                            gCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgc                                                                                      >  1:555390/1‑150 (MQ=255)
                                                            gCAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgc                                                                                      <  1:8836472/150‑1 (MQ=255)
                                                             cAACTGGCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTcgcc                                                                                     >  1:5165528/1‑150 (MQ=255)
                                                                  ggCCGAGACCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCa                                                                                >  2:3245693/1‑150 (MQ=255)
                                                                      gagaCCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTc                                                                            <  2:7362456/150‑1 (MQ=255)
                                                                       agaCCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCg                                                                           <  2:3443657/150‑1 (MQ=255)
                                                                        gaCCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCgg                                                                          >  1:10297494/1‑150 (MQ=255)
                                                                        gaCCGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCgg                                                                          <  1:8460883/150‑1 (MQ=255)
                                                                           cGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGcct                                                                       <  1:8414736/150‑3 (MQ=255)
                                                                           cGTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGcct                                                                       >  1:7046954/1‑148 (MQ=255)
                                                                            gTGGTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTg                                                                      >  1:6278639/1‑150 (MQ=255)
                                                                              ggTCCTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGAACTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGga                                                                    <  2:1252179/150‑1 (MQ=255)
                                                                                 ccTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCCCCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacga                                                                 <  2:4193689/150‑1 (MQ=255)
                                                                                 ccTGAACAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTATCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCCACTCGGCCTGgacga                                                                 >  2:3679125/1‑150 (MQ=255)
                                                                                     aaCAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGgacgacgac                                                             <  2:5251782/150‑1 (MQ=255)
                                                                                       cAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGAcgc                                                           <  1:279000/150‑1 (MQ=255)
                                                                                       cAAGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGAcgc                                                           >  2:7017277/1‑150 (MQ=255)
                                                                                        aaGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGAcgcg                                                          <  2:6722101/150‑1 (MQ=255)
                                                                                         aGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCgg                                                         <  1:7375071/150‑1 (MQ=255)
                                                                                         aGGCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCgg                                                         <  1:1135051/150‑1 (MQ=255)
                                                                                          ggCCCGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGt                                                        <  1:9694317/150‑1 (MQ=255)
                                                                                            cccGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGc                                                      <  1:3931726/150‑1 (MQ=255)
                                                                                             ccGCCTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCt                                                     <  1:10778696/150‑1 (MQ=255)
                                                                                                 cTGGGCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctg                                                 >  2:10150723/1‑150 (MQ=255)
                                                                                                   gggCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgct                                               >  2:5628027/1‑150 (MQ=255)
                                                                                                    ggCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGGCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCctgctg                                              >  2:5769262/1‑150 (MQ=255)
                                                                                                     gCCTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCtgctgc                                             <  2:9406249/150‑1 (MQ=255)
                                                                                                       cTGTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCaa                                           >  1:5805814/1‑150 (MQ=255)
                                                                                                         gTCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAg                                         <  1:980065/150‑1 (MQ=255)
                                                                                                          tCCTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGc                                        <  1:200818/150‑1 (MQ=255)
                                                                                                           ccTGGCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCa                                       <  2:1695713/150‑1 (MQ=255)
                                                                                                               gCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATccc                                   <  2:1175685/150‑1 (MQ=255)
                                                                                                               gCAGGATCTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATccc                                   <  1:4975258/150‑1 (MQ=255)
                                                                                                                      cTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTgcgc                            <  2:6637468/150‑1 (MQ=255)
                                                                                                                      cTGGCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTgcgc                            <  1:3120002/150‑1 (MQ=255)
                                                                                                                        ggCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCgg                          >  1:6081815/1‑150 (MQ=255)
                                                                                                                         gCCGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgc                         <  1:9563091/150‑1 (MQ=255)
                                                                                                                          ccGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgca                        <  1:9345824/150‑1 (MQ=255)
                                                                                                                          ccGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgca                        <  1:6508804/150‑1 (MQ=255)
                                                                                                                           cGACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcag                       <  1:7531753/150‑1 (MQ=255)
                                                                                                                            gACGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcagc                      <  1:1749142/150‑1 (MQ=255)
                                                                                                                             aCGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcagca                     >  1:8937534/1‑150 (MQ=255)
                                                                                                                             aCGGCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGgcagca                     <  2:7686418/150‑1 (MQ=255)
                                                                                                                               ggCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCAtt                   >  1:2301997/1‑150 (MQ=255)
                                                                                                                                gCACCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTc                  <  2:7389179/150‑1 (MQ=255)
                                                                                                                                  aCCGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTccc                <  1:5663699/150‑1 (MQ=255)
                                                                                                                                   ccGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCcg               >  2:7064215/1‑150 (MQ=255)
                                                                                                                                    cGGTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCcgg              >  1:10575048/1‑150 (MQ=255)
                                                                                                                                     ggTCTCGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCcggc             >  2:8551707/1‑150 (MQ=255)
                                                                                                                                          cGCCCTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCAt        <  2:5897419/150‑1 (MQ=255)
                                                                                                                                            cccTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCACGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTc      >  2:8467245/1‑150 (MQ=255)
                                                                                                                                            cccTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTc      >  2:4841199/1‑150 (MQ=255)
                                                                                                                                            cccTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTc      >  2:1101898/1‑150 (MQ=255)
                                                                                                                                            cccTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTc      >  1:9051543/1‑150 (MQ=255)
                                                                                                                                            cccTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTc      <  2:6107164/150‑1 (MQ=255)
                                                                                                                                              cTGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTccc    <  1:3475076/150‑1 (MQ=255)
                                                                                                                                               tGACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCa   <  1:1273438/150‑1 (MQ=255)
                                                                                                                                                gACCTTCGTCACCGCCGCCCTGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAg  <  2:10430657/150‑1 (MQ=255)
                                                                                                                                                gACCTTCGTCACCGCCGCCCGGCTCGGCCAGCATGCCCTGCCAGAGGCGGCGGCGCGCAAGGTCGCCGCCCAACTCGGCCTGGACGACGACGCGGTGCTCCTGCTGCAAAGCATCCCGCTGCGCGGCAGCATTCCCGGCGGCATTCCCAg  <  2:10578468/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: